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EXECUTIVE  SUMMARY 


High  frequency  (HF)  links  (2  to  30  MHz)  are  an  alternative  to  the  cost  and  tactical  fragility  of 
commercial  satellite  communications  (SATCOM)  for  beyond  line-of-sight  (LOS)  links.  However,  standard 
HF  systems  operating  over  a  3 -kHz  bandwidth  do  not  provide  sufficient  throughput  for  network 
applications.  Simply  increasing  bandwidth  does  increase  throughput — provided  the  radios  are  properly 
designed  for  the  time-varying  wideband  HF  channel.  The  wideband  HF  channel  is  governed  by  the 
ionospheric  variations  caused  by  latitude,  the  night  and  day  cycle,  the  direction  and  length  of  the  haul,  the 
seasons,  and  ionosphere  disturbances.  This  report  is  the  first  of  a  series  seeking  a  better  understanding  of 
each  of  these  channel  variations  to  aid  in  the  design  of  wideband  HF  systems  and  the  proper  comparison  of 
different  HF  radios  and  waveforms. 

This  report  uses  wideband  mid-latitude  HF  channel  soundings  and  three-dimensional  (3-D)  ray-tracing 
simulations  to  develop  a  statistical  model  of  a  particular  nearly  vertical  incidence  skywave  (NVIS)  channel. 
This  model  is  expressed  in  a  form  that  is  programmable  in  a  channel  emulator  so  that  different  radios  can 
be  tested  under  controlled  conditions.  The  3-D  ray-tracing  code  is  validated  by  measurements  with  the 
intention  of  running  the  code  over  different  times  and  geographic  locations  to  generate  channels  that  have 
not  been  measured  over  the  air.  Then,  those  simulated  channels  will  be  used  to  compute  statistical  models 
for  particular  geographic,  temporal,  and  environmental  conditions. 

The  immediate  application  of  this  particular  channel  is  for  small  unit  mobile  ground  forces  with 
vehicle-mounted  antennas.  Near-term  extensions  support  mid-latitude  short-haul  (200  to  1,500  km)  and 
long-haul  ship-to-ship  and  ship-to-shore  wideband  HF  channels.  Long-term  payoffs  for  the  warfighter  will 
be  a  better  HF  radio  design,  leading  to  increased  data  rates,  more  stable  networks,  reduced  power 
requirements,  and  improved  low  probability  of  detection  (LPD). 
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1.  REPORT  OUTLINE 


This  report  analyzes  a  mid-latitude  wideband  high  frequency  nearly  vertical  incidence  skywave  (HF- 
NVIS)  channel.  The  immediate  objective  is  a  wideband  HF-NVIS  channel  model  suitable  for  software 
simulation  and  hardware  emulation.  The  long-term  objective  is  to  produce  a  three-dimensional  (3-D)  ray¬ 
tracing  code  suitable  for  assessing  HF  waveforms,  HF  antenna  designs  including  polarizations,  automatic 
link  establishment  (ALE),  and  stealthy  HF  communications,  and  modeling  credible  HF  channels  in  regions 
of  tactical  interest.  This  report  is  the  first  of  a  series  to  populate  the  wideband  HF-NVIS  channel  portfolio 
discussed  in  the  Executive  Summary.  Figure  1  shows  a  high-level  view  of  this  process  proceeding  from 
channel  validation  to  statistical  modeling  to  wideband  HF  channel  emulation. 


Simulated  channels 

Ray-tracing  code  using  real  ionograms 

with  physical  Validated  with  HF  over  many  locations, 

ionosphere  model  measurements  times  of  day,  and 


times  of  year 


Figure  1.  HF-NVIS  roadmap:  theory,  simulations,  and  measurements  delivering  channel  models  suitable 
for  emulation. 

Section  2  provides  a  context  for  this  roadmap  and  the  follow-on  reports.  The  initial  Office  of  Naval 
Research  (ONR)  requirements  for  mobile  HF-NVIS  links — a  mobile  satellite  communications  (SATCOM) 
“lite” — are  quantified  using  a  basic  link  budget.  This  link  budget  shows  that  the  vehicle-sized  antenna  is 
the  link’s  bottleneck.  Accordingly,  recent  antenna  designs  are  reviewed  to  bound  the  state  of  the  art.  The 
upshot  is  that  further  gains  must  come  from  HF-NVIS  channel  knowledge. 

Section  3  reviews  HF-NVIS  channels  using  ray  paths.  Figures  2  and  3  show  mid-latitude  HF-NVIS  ray 
paths  linking  a  transmitter  to  its  receiver.  These  ray  paths  are  computed  by  the  3-D  ray-tracing  code 
detailed  in  Appendix  D.  These  ray  paths  are  determined  by  the  electron  density,  the  magnetic  field,  and  the 
ions  and  neutral  particles  in  a  static  and  spherically  homogeneous  ionosphere  determined  from 
neighboring  ionograms.  The  O-  and  X-Modes  are  produced  by  the  index  of  refraction  and  the  Earth’s 
magnetic  field  (O-Mode  bends  north;  X-Mode  bends  south).  The  novelty  of  this  section  is  the  linking  of 
polarization  evolution  along  the  ray  path  and  the  consequences  for  multiple  hops.  This  multi-hop  analysis 
is  applied  to  the  measured  channel  soundings  to  determine  their  delay  structure.  Once  the  ray  paths  are 
known,  the  propagation  losses — spreading,  mode- shedding,  and  absorption — may  be  computed  under  the 
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assumption  these  losses  have  not  distorted  the  ray  paths.  Together,  with  the  antenna  coupling  into  these  ray 
paths,  the  3-D  ray  tracing  permits  the  coherent  summation  of  the  electric  field  to  capture  Faraday  fading 
effects,  O-  and  X-Mode  processing,  dual-polarization  processing,  and  the  antenna  coupling  for  single  input 
multiple  output  (SIMO)  and  multiple  input  multiple  output  (MIMO)  links. 


Figure  2.  Single  HF-NVIS  hop.  Figure  3.  FIF-NVIS  2-Flop. 

Section  4  generalizes  the  discussion  to  an  ionosphere  that  is  time-varying  and  spherically 
inhomogeneous.  In  this  case,  the  single  and  multiple  hops  of  a  nominal  static  ionosphere  are  subject  to  fast 
and  slow  fading.  The  section  starts  by  reviewing  time- varying  channels  to  make  explicit  the  scattering 
function,  the  delay  spread,  the  Doppler  spread,  the  Doppler  shift,  and  the  fading  processes.  HF-NVIS 
channels  cluster  in  two  classes  in  the  Doppler-delay  plane: 

•  Mid-latitude  channels  show  small  Doppler  spread  and  large  delay  spread 

•  High-latitude  channels  show  large  Doppler  spread  and  small  delay  spread 

Section  4.4  lists  the  narrowband  mid-latitude  HF-NVIS  channels  suitable  for  programming  an  HF 
channel  emulator.  Section  4.5  reviews  the  Doppler-delay  spreads  for  high  latitudes.  These  surveys  show 
relatively  few  wideband  HF-NVIS  channel  models  exist  in  a  form  suitable  for  HF  channel  emulators. 

Sections  5,  6,  and  7  pertain  to  the  soundings  and  analysis  of  a  mid-latitude  channel  on  10  February 
2016  at  5.82  MHz  over  the  25-kHz  band.  Section  5  reviews  the  HF-NVIS  channel  sounding  process: 
antennas,  hardware,  waveform,  and  channel  estimation.  Section  6  documents  the  signal,  noise,  and 
interference  along  the  route.  Section  7  is  the  channel  analysis.  Selected  channels  along  the  route  are 
examined  in  detail  using  the  3-D  ray-tracing  code  to  untangle  closely  spaced  O-  and  X-Modes  on  the  first 
hop,  identify  multihops,  and  identify  the  ray  paths  that  are  not  multiples  of  the  first  hop.  The  channel 
statistics  suitable  for  a  channel  emulator  are  extracted.  This  section  concludes  by  processing  all  the 
channels  window  to  deliver  an  “average”  wideband  HF-NVIS  channel  ensemble  suitable  for  a  channel 
emulator. 

Section  8  lays  out  selected  HF  system  analyses  made  possible  by  this  3-D  ray-tracing  code: 

•  Augmenting  HF  Portfolios 

•  O-  and  X-Mode  diversity  gain 

•  Over-the-water  multipath 

•  Multiple  channels  for  HF  Internet  Protocol  (HF-IP) 
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The  close  agreement  between  the  measured  and  simulated  HF  channels  in  this  report  supports  the 
claim  that  the  costly,  time-consuming,  and  non-repeatable  over-the-air  testing  may  be  augmented  by 
simulation  to  populate  HF  channel  ensembles  supporting  HF  system  assessment  both  in  simulation  and 
emulation. 

1.1  CURRENT  RAY-TRACING  LIMITATIONS 

The  3-D  ray-tracing  code  of  this  report  includes  the  O-  and  X-Modes,  multiple  hops,  antennas,  and 
polarization  along  the  ray  path  to  coherently  sum  the  complex-valued  responses  of  the  receive  antenna  to 
the  arriving  electric  field.  Ray-tracing  codes  that  track  power  only  cannot  model  multipath  HF  channel 
responses  nor  encompass  the  multi-antenna  designs.  Nevertheless,  all  ray-tracing  codes  require  the 
electron  density  profile  and  magnetic  field.  Ideally,  ray-tracing  simulations  use  the  time-varying  electron 
density  and  magnetic  fields  at  each  spatial  location  along  the  ray  path.  Currently,  this  3-D  ray-tracing  code 
extracts  electron  densities  from  ionograms  at  selected  locations  where  the  soundings  are  taken  in  time 
increments  ranging  from  5  to  30  minutes.  Consequently,  this  electron  density  estimate  is  static  and  cannot 
model  the  time-varying  Doppler  or  fading  along  ray  paths.  Similarly,  the  magnetic  field  is  obtained  from 
the  Earth’s  magnetic  field  and  is  also  a  static  model.  Including  these  time-varying  ionospheric  processes, 
especially  for  mission  planning,  is  a  substantial  challenge.  Therefore,  this  preface  closes  with  the 
illustration  of  the  time-varying  ionospheric  processes. 

1.2  IONOSPHERIC  PROCESSES 

Figure  4  is  a  superb  graphic  relating  multiple  ionospheric  processes.  Solar  activity  produces  layers  of 
varying  electron  densities  on  the  sun-facing  side  of  the  ionosphere.  The  night-facing  side  of  the  ionosphere 
has  a  reduced  density  and  fewer  layers.  In  this  report,  the  mid-latitude  region  is  relatively  stable  with  only 
one  traveling  ionospheric  disturbance  (TID)  observed  in  these  soundings;  mid-latitude  irregularities — the 
2.5  ray  paths  observed  in  the  soundings — are  likely  caused  by  sporadic  E-layers.  The  high-latitudes  show 
complex  magnetic  activity:  Doppler  spreads  up  to  16  Hz  are  reported  in  the  literature  [73];  HF-NVIS 
propagation  can  break  down  in  the  polar  regions — the  HF  links  propagate  via  the  polar  walls  rather  than 
overhead  [73];  and  Doppler  effects  require  both  time  variation  and  the  frequency  variation  of  the  electron 
density  [36].  Nevertheless,  the  Hamiltonian  encoded  in  the  3-D  ray  tracing  supports  time-varying  effects. 


Figure  4.  Ionospheric  and  thermospheric  processes  (This  file  is  in  the  public  domain  in  the  United  States 
because  it  was  solely  created  by  NASA). 
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In  theory,  simulation  of  HF  system  performance  is  possible — provided  physics-based  time- varying 
models  of  the  ionospheric  processes  are  available.  Maslin’s  observation  [46,  Section  8.2.3]  states 
“predictions  can  only  be  as  good  as  the  model,”  which  requires  the  FIF  ensembles  encompass  the  modeling 
extrema.  Consequently,  Section  9  details  approaches  to  bring  these  time-varying  models  into  this  HF 
portfolio. 

Tables  1 ,  2,  and  3  define  the  symbols  used  for  ionospheric  parameters,  the  plasma  frequencies  and 


layer  heights,  and  HF  channel  parameters,  respectively. 

Table  1 .  Ionospheric  parameters. 

Bo 

Magnitude  of  the  Earth’s  magnetic  induction  field  (Tesla) 

eo 

Permittivity  of  free  space  (F/m) 

fp 

Plasma  frequency  (Hz) 

Id 

Doppler  shift  (Hz) 

A  / 

Doppler  spread  (Hz) 

me 

Electron  mass  (kg) 

Ne 

Electron  density  (number  of  electrons  per  cubic  meter) 

Nn 

Neutral  particle  density  (number  of  neutral  particles  per  cubic  meter) 

Qe 

Charge  of  an  electron  (C) 

T 

- L  e 

Electron  temperature  (K) 

Ulp 

Plasma  frequency  (rad/sec) 

UJr 

Gyro-frequency  (rad/sec) 

uje 

Electron  collision  frequency  (rad/sec) 

me, i 

Electron-ion  collision  frequency  (rad/sec) 

UE.N 

Electron-neutral  collision  frequency  (rad/sec) 

Table  2.  Frequencies  and  heights  characterizing  the  electron  density  profile. 

foF2  F2  layer  plasma  frequency  (MHz) 
foFI  FI  layer  plasma  frequency  (MHz) 
foE  E  layer  plasma  frequency  (MHz) 
hmF2  F2  layer  Ne  maximum  height  (km) 
hmFI  FI  layer  Ne  maximum  height  (km) 
hmE  E  layer  Ne  maximum  height  (km) 


Table  3. 

HF-NVIS  channel  parameters. 

fc 

Center  frequency  (MHz) 

fp 

Bandwidth  (kHz) 

Id 

Doppler  shift  (Hz) 

A  / 

Doppler  spread  (Hz) 

T 

Delay  (sec) 
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ACRONYMS 


ALE  Automatic  Link  Establishment 
BPSK  Binary  Phase-Shift  Keying 
CCIR  Consultative  Committee  on  International  Radio 
DSTG  Defense  Science  and  Technology  Group 
EFV  Expeditionary  Fighting  Vehicle 
EM  Electromagnetic 
GPS  Global  Positioning  System 
GPSDO  GPS-disciplined  Oscillator 
HF  High  Frequency 
HF-IP  HF  Internet  Protocol 
HMMWV  High-Mobility  Multipurpose  Wheeled  Vehicle 
IEEE  Institute  of  Electrical  and  Electronics  Engineers 
LCP  Left  Circularly  Polarized 
LEP  Left  Elliptical  Polarized 
LNA  Low-noise  Amplifier 
LOS  Line-of-sight 
LPD  Low  Probability  of  Detection 
M-ATV  Mine-Resistant  Ambush  Protected  All-Terrain  Vehicle 
MIMO  Multiple  Input  Multiple  Output 
NB  Narrowband 

NVIS  Near  Vertical  Incidence  Skywave 
O-Mode  Ordinary  Mode 

ODE  Ordinary  Differential  Equation 
OFDM  Orthogonal  Frequency  Division  Multiplexing 
ONR  Office  of  Naval  Research 
PA  Power  Amplifier 

PMFS  Phase  Modulation  Fading  Simulator 
PN  Pseudorandom  noise 
PST  Pacific  Standard  Time 
QMFS  Quadrature  Modulation  Fading  Simulator 
RCP  Right  Circularly  Polarized 
REP  Right  Elliptical  Polarized 
RIFF  Resource  Interchange  File  Format 
RF  Radio  Frequency 
SATCOM  Satellite  Communications 

SBIR  Small  Business  Innovation  Research 
SDR  Software-defined  Radio 
SIMO  Single  Input  Multiple  Output 
SISO  Single  Input  Single  Output 
SNR  Signal-to-Noise  Ratio 
SPAWAR  Space  and  Naval  Warfare  Center 
TID  Traveling  Ionospheric  Disturbance 
USRP  Universal  Software  Radio  Peripheral 
VOACAP  Voice  of  America  Coverage  Analysis  Program 
VSWR  Voltage  Standing  Wave  Ratio 
WB  Wideband 
X-Mode  Extraordinary  Mode 
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2.  MOBILE  WIDEBAND  HF-NVIS  COMMUNICATION  DESIGN 


Sklar  observes  that  radio  design  intertwines  requirements  and  the  channel  [72]: 

The  design  of  any  digital  communication  systems  begins  with  a  description  of  the  channel 
(receiver  power,  available  bandwidth,  noise  statistics, . . .  )  and  a  definition  of  the  system 
requirements  (data  rate  and  error  performance). 

This  report  adapts  Sklar’s  system-level  view  for  HF-NVIS  operating  on  a  mobile  platform.  Section  2.1 
sets  out  the  mobile  wideband  HF-NVIS  requirements.  Section  2.2  reviews  HF-NVIS  propagation  and  the 
antenna  patterns  necessary  to  send  and  receive  on  the  HF-NVIS  links.  Section  2.3  analyzes  the  tradeoffs  to 
operate  HF-NVIS  from  a  ground  vehicle.  Section  2.4  reviews  progress  in  wideband  HF  relevant  to  mobile 
wideband  HF-NVIS.  The  vehicle-sized  antennas  are  the  bottleneck  of  these  HF-NVIS  systems.  Therefore, 
understanding  the  HF-NVIS  channel  is  necessary  to  optimize  the  link  from  these  power-  and  bandwidth- 
limited  platforms. 

2.1  MOBILE  WIDEBAND  HF-NVIS  REQUIREMENTS 

ONR  observed  communication  shortfalls  in  expeditionary  and  irregular  warfare  [49]: 

•  SATCOM  is  easy  to  jam 

•  Mountainous  terrain  blocks  line-of-sight  (LOS)  communications 

•  An  HF  skip  zone  exists  that  causes  certain  distances  to  be  unreachable  by  HF  waves  when  using 
certain  antenna  types 

One  solution  exploits  the  nearly  vertical  skywaves  of  HF-NVIS  to  lift  data  “up  and  over”  terrain,  cover 
the  HF  skip  zone,  and  offer  a  SATCOM  “lite”  link.  Therefore,  ONR’s  Improved  HF  Communications 
Small  Business  Innovation  Research  (SBIR)  sought  a  link  suitable  for  Small  Unit  Communications 
Platforms  using  the  following  requirements: 

•  Power:  ground  vehicles  (150  Watts);  manpack  (20  Watts) 

•  Throughput:  100  kbps 

•  Range:  0  to  100  miles 

•  Frequency:  2  to  12  MHz 

•  Bandwidth:  3  to  24  kHz 

Implicit  in  “mobile”  is  the  requirement  that  the  wideband  HF-NVIS  antenna  fits  on  an  operational 
vehicle  and  does  not  interfere  with  mission  requirements.  Section  2.4  illustrates  current  and  proposed 
vehicle-mounted  HF-NVIS  antennas.  Conventional  HF  bandwidth  is  3  kHz.  Wider  bands  typically  step  in 
multiples  of  3-kHz  up  to  24-kHz  (MIL-STD- 188-1 10C  Appendix  D  [26]).  This  report  uses  the  24-kHz 
bandwidth. 
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2.2  NARROWBAND  HF-NVIS 


The  NVIS  Army  Field  Manual  24-18  provides  an  excellent  summary  of  conventional  HF  shortfalls 

[20]: 


Sky  waves,  generated  by  standard  antennas  (for  example,  doublets)  which  efficiently  launch 
the  sky  wave,  will  not  return  to  earth  at  a  range  of  less  than  161  kilometers  (100  mi).  This  can 
leave  a  skip  zone  of  at  least  80  to  113  kilometers  (50  to  70  mi)  where  HF  communications  will 
not  function.  This  means  that  units  such  as  long-range  patrols,  armored  cavalry  deployed  as 
advance  or  covering  forces,  air  defense  early  warning  teams,  and  many  division-corps, 
division-brigade,  division-DISCOM  and  division-DIVARTY  stations  are  in  the  skip  zone  and 
thus  unreachable  by  HF  radio  even  though  HF  is  a  primary  means  of  communication  to  these 
units. 

Figure  5  shows  a  standard  HF  antenna  pattern  that  launches  most  of  the  HF  energy  at  approximately 
20°  elevation  angle.  Figure  6  illustrates  the  skip  zones  produced  by  standard  HF  antennas. 


Take-Off  Angle 


15  10  5  0  -5  -10  -5  0  5  10  15 

dBi 


LEGEND 

- 3  Mhz  -  9  Mhz  - 18Mhz 

dBi  gain  in  decibels  Mhz  Megahertz 


Figure  C-1.  32-foot  vertical  whip,  vertical  antenna  pattern 

Figure  5.  Standard  HF  antenna  pattern  [19]. 


Figure  6.  Sky  wave  propagation  showing  the  skip 
zone  [19]. 


The  Army  Field  Manual  24-18  describes  how  HF-NVIS  propagation  overcomes  the  skip  zone  [20]: 

Energy  radiated  in  a  near- vertical  incidence  direction  is  not  reflected  down  to  a  pinpoint  on  the 
Earth’s  surface.  If  it  is  radiated  on  too  high  a  frequency,  the  energy  penetrates  the  ionosphere 
and  continues  on  out  into  space.  Energy  radiated  on  a  low  enough  frequency  is  reflected  back 
to  earth  at  all  angles  (including  the  zenith),  resulting  in  the  energy  striking  the  earth  in  an 
omnidirectional  pattern  without  dead  spots  (that  is,  without  a  skip  zone).  Such  a  mode  is 
called  a  near- vertical  incidence  sky  wave  (NVIS). 

Figure  7  shows  an  antenna  pattern  radiating  most  of  its  power  in  the  upward  direction.  Figure  8 
illustrates  the  nearly  vertical  skywaves  refracting  off  the  ionosphere  and  illuminating  the  terrain  from  zero 
to  100  km.  Horizontal  dipoles,  rhombic  antennas,  V  antennas,  and  loops  are  common  geometries  for 
HF-NVIS  antennas.  Figures  9  and  10  illustrate  that  these  standard  HF-NVIS  antenna  scale  with 
wavelength  (e.g.,  the  horizontal  dipole  is  98  ft  long  at  5  MHz).  Thus,  the  challenge  for  mobile  WBHF  is  to 
design  a  vehicle-sized  antenna. 
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Figure  M-4.  Typical  elevation  plane  pattern. 


Figure  M-1.  Near-vertical  incidence  sky-wave  propagation 
concept. 


Figure  7.  Atypical  upward-pointing  HF-NVIS 
antenna  pattern  [1 9]. 


Figure  8.  HF-NVIS  refracting  off  of  the  ionosphere 
[19]. 


Figure  3- 1 2.  Center -fed  Hertz  antenna  with  two  upright  supports 

Figure  9.  Horizontal  dipole  [19].  Figure  10.  Rhombic  antenna  [19]. 

2.3  A  MOBILE  WIDEBAND  HF-NVIS  LINK  BUDGET 

This  section  steps  through  a  basic  link  budget  to  quantify  the  antenna  bottleneck  of  the  mobile 
HF-NVIS  link.  The  antenna  dominates  the  link.  Consequently,  only  a  simplified  HF-NVIS  link  budget  is 
necessary  [77].  All  antennas  scale  to  wavelength.  The  wavelengths  in  HF  range  from  10  to  100  m  as 
frequency  sweeps  down  from  30  to  3  MHz.  The  vertical  whip  cut  to  a  quarter  wavelength  sizes  from  2.5  to 
25  m  (e.g.,  the  32-foot  vertical  whip  in  Figure  5  operates  at  7.7  MHz,  assuming  a  quarter  wavelength). 
HF-NVIS  operates  in  2  to  12  MHz  as  pointed  by  the  NVIS  Army  Field  Manual  24-18  [20]: 

The  main  difference  between  this  short-range  NVIS  mode  and  the  standard  long-range  sky- 
wave  HF  mode  is  the  lower  frequency  required  to  avoid  penetrating  the  ionosphere  and  the 
angle  of  incident  signal  upon  the  ionosphere.  To  attain  an  NVIS  effect,  the  energy  must  be 
radiated  strong  enough  at  angles  greater  than  about  75  or  80  degrees  from  the  horizontal  on  a 
frequency  that  the  ionosphere  will  reflect  at  that  location  and  time. 

Antennas  do  not  necessarily  have  to  be  exactly  a  quarter  or  half  wavelength  to  operate,  but  antennas 
smaller  than  this  standard  (electrically-small  antennas)  suffer  from  decreased  gain  and  bandwidth. 
Therefore,  the  long  wavelengths  required  for  HF-NVIS  make  antenna  size  problematic  for 
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vehicle-mounted  antennas.  Figure  1 1  shows  a  59-in  diameter  loop  antenna  superimposed  where  it  would 
be  mounted  on  a  High-Mobility  Multipurpose  Wheeled  Vehicle  (HMMWV).  This  loop  is  the  exemplar 
illustrating  trade-offs  required  for  mobile  HF-NVIS  operations.  An  idealization  of  this  loop,  idealized 
matching  circuits,  and  standard  noise  models  constitute  the  HF-NVIS  link  model  analyzed  in  this  section. 


Figure  1 1 .  HMMVW  with  HF-NVIS  loop  superimposed  in  yellow  of  the  same  dimensions  as  the  Harris 
RF-3134-AT003  antenna  (photo  credit:  Claire  Heininger,  U.S.  Army). 

Figure  12  is  a  high-level  block  diagram  of  the  HF-NVIS  link.  For  this  discussion,  the  HF-NVIS 
channel  is  modeled  as  a  lossy  channel  with  identical  antenna  and  matching  circuits  at  the  transmitter  and 
receiver. 


Figure  12.  Block  diagram  of  the  HF-NVIS  link.  Signal  power:  pT  is  power  from  the  transmitting  amplifier, 
Pa  is  power  delivered  to  the  antenna,  prad  is  power  radiated,  px  is  signal  power  captured  by  receive 
antenna,  and  pR  is  signal  power  remaining  after  matching  circuit.  Gains:  GT  is  transducer  power  gain  of 
the  matching  circuit  and  Ga  is  absolute  gain  of  the  antenna.  Noise:  pN) RF  is  external  radio  frequency  (RF) 
noise  and  pjv,lo  internal  receiver  noise. 
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The  HF-NVIS  link  must  have  a  sufficiently  large  signal-to-noise  ratio  (SNR).  The  SNR  is  measured  at 
the  input  to  the  receiver’s  detector  is 


SNR 


_ PR _ 

Pn,  lo  +  GtG  aPn,  rf 


The  local  noise  pn, lo  is  assumed  to  be  negligible  compared  to  the  received  noise  power  GtG  aPn, rf- 
Consequently,  the  SNR  admits  the  approximation 


SNR 


PR 

GtG  a  Pn, rf 


GtG  a  Prf  _  Prf  _  GufGtGa  Pt 
GtG  a  Pn, rf  Pn,  rf  Pn, rf 


(1) 


The  antenna  gain,  matching  circuits,  and  the  noise  levels  are  reviewed  to  determine  WB HF-NVIS  trade¬ 
offs. 

Antenna  Gain:  The  trade-off  discussion  must  start  with  a  particular  antenna.  Figure  13  shows  a  simulated 
loop  antenna  and  its  gain  pattern  at  12  MHz  over  a  perfect  ground  plane.  There  is  no  vehicle  in  this 
simulation.  Thus,  the  antenna’s  nulls  precisely  align  along  the  x-axis.  Figure  14  compares  gain  patterns 
across  the  HF  band.  Referring  to  Figure  12,  the  antenna  gain  G a  is  the  ratio  of  the  radiated  power  prad  to 
the  power  available  pa  to  the  antenna: 

q  _  Praf> 

PA 

Figure  14  shows  that  antenna  gain  rapidly  falls  off  as  frequency  decreases. 


z 

X 


Figure  13.  Simulated  antenna  pattern  of  a  lossy 
vertical  loop  with  a  59-in  diameter  mounted  four  ft 
above  a  perfect  ground  plane. 


Figure  14.  Simulated  gain  of  the  lossy  vertical  loop 
above  perfect  ground  sliced  along  the  y-z  plane; 
Ga  =  0  dB  at  5  MHz. 
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Figure  15  plots  the  simulated  reflectance  obtained  from  the  simulated  vertical  loop  antenna.  The  Smith 
chart  shows  this  antenna  behaves  as  a  “one-coil”  inductor:  The  loop’s  reflectance  starts  as  an  inductor  at  2 
MHz  and  follows  the  induction  circle  counter-clockwise  around  the  Smith  chart  to  the  open  circuit  with 
increasing  frequency.  An  electrically  small  antenna  requires  a  matching  circuit.  Increasing  the  bandwidth 
severely  challenges  the  matching  circuit.  Figure  16  reports  the  matching  at  5  MHz  over  the  25-kHz  band 
over  all  possible  three-stage  LC  ladders.  The  circuit  decreases  the  maximum  voltage  standing  wave  ratio 
(VSWR)  over  the  25-kHz  band  from  5,177  down  to  22.4.  Equivalently,  this  lossless  matching  sets  the 
matched  transducer  power  gain  Gt  —  —  7.9  dB. 


Figure  15.  Reflectance  of  the  simulated  59-in  Figure  16.  Matching  the  simulated  59-in  diameter 

diameter  vertical  loop  of  Figure  1 4;  The  circle  marks  vertical  loop  at  5  MHz  over  25  kHz  using  a  third- 
the  reflectance  at  2  MHz;  the  end  of  the  blue  line  order  ladder, 
marks  the  reflectance  at  13  MHz. 

Table  4  puts  this  single  matching  circuit  in  context.  The  table  reports  on  this  matching  for  12-,  25-,  and 
50-kHz  bandwidths.  The  line  labeled  “Ladder”  lists  the  transducer  power  gain  over  all  possible  LC  ladders 
containing  at  most  three  inductors  or  capacitors.  The  line  labeled  “State-space”  is  the  best  possible  match 
over  all  lossless  two-ports  that  contain  at  most  three  reactive  elements  [67].  Because  the  state-space 
matching  class  contains  the  LC  ladders,  the  state-space  gain  exceeds  the  ladder  gain.  Looking  at  the  table 
from  top  to  bottom,  the  state-space  matching  shows  there  exist  two-port  topologies  that  deliver  at  least  2  dB 
more  gain  over  the  LC  ladders.  Looking  at  the  table  from  right  to  left,  doubling  the  bandwidth  costs  3  dB 
regardless  of  the  matching  topology.  Because  the  state-space  matching  class  contains  the  LC  ladders,  the 
state-space  gain  exceeds  the  ladder  gain.  Looking  at  the  table  from  top  to  bottom,  the  state-space  matching 
shows  there  exist  two-port  topologies  delivering  at  least  2  dB  more  gain  over  the  LC  ladders.  Looking  at  the 
table  from  right  to  left,  doubling  the  bandwidth  costs  3  dB  regardless  of  the  matching  topology.  The 
forthcoming  link  budget  for  the  25-kHz  bandwidth  uses  the  general  matching  Gt=  —  5.4  dB. 

Table  4.  Transducer  power  gain  GT  (dB)  for  degree  3  matching  at  5  MHz  of  the  simulated  59-in  diameter 
vertical  loop  antenna.  GT  is  the  worst  case  gain  in  the  bandwidth. 


Matching  12  kHz  25  kHz  50  kHz 
Ladder  ^5T  -10.7 

State-space  -2.8  -5.4  -8.3 
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RF  Noise:  The  external  RF  noise  is  modeled  as  zero-mean,  complex- valued  Gaussian  with  flat  spectrum 
across  the  frequency  band.  As  such,  only  the  variance  or,  equivalently,  in-band  noise  power  pn:rf,  is 
required  to  model  this  Gaussian  noise.  This  noise  power  is  modeled  as  [40,  page  76] 

PN,  RF  =  KBTeqfB ,  (2) 

where  Boltzmann’s  constant  is  kb  —  1-38  x  10-23  (W  •  K-1  •  Hz-1),  Teq  is  the  effective  operating 
temperature  (K),  and  fs  is  the  bandwidth  (Hz).  Figure  17  plots  Teq  as  a  function  of  frequency  for  selected 
RF  noise  scenarios.  The  “Urban  Man-Made  Noise”  sets  the  noise  level  for  these  WBHF-NVIS  links. 

Figure  18  plots  the  urban  noise  power  modeled  by  Equation  (2)  for  three  frequency  bands. 


RF  Noise:  Urban  Man  Made 


Figure  17.  Sky  temperatures  Teq  (K)  generated  by  Figure  18.  Urban  wideband  HF  noise  power  from 
background  noise  sources  [40,  page  74].  Equation  (2). 


HF  Gain:  This  report  bounds  the  gain  along  an  HF  ray  path  as  the  product  of  the  spreading,  absorption, 
and  mode- shedding  gains: 

Ghf  =  Gspread  x  G  absorb  x  Gshed-  (3) 

The  spreading  or  free-space  gain  sets  a  convenient  upper  bound  on  the  HF  gain  [87,  Eq.  3.23]: 

A 

47T5 

where  s  is  the  length  of  the  ray  path.  A  single  ray  path  with  a  500-km  length  has  a  spreading  gain  of 
approximately  Gspread  ~  —100  dB  at  5  MHz.  Equation  4  assumes  standard  spherical  spreading 
approximating  the  more  realistic  astigmatic  spreading  caused  by  refraction  through  the  ionosphere  [55]. 
Absorption  and  mode-shedding  gains  are  discussed  in  later  sections.  Only  this  standard  spreading  bounds 
the  link  budget. 


Ghf  <  Gspread  = 
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HF-NVIS  Link:  Transmit  power  px  and  SNR  tradeoffs  are  governed  by  Equation  (1): 

SNR  <  Gt  +  Ga  +  Ghf  +  Pt  —  Pn, rf  [dB].  (5) 

Equation  (5)  directs  attention  to  those  parts  of  the  HF-NVIS  link  that  can  improve  these  Power-SNR 
tradeoffs.  The  transmit  power  px,  the  RF  noise  pn, rf,  and  the  HF  link  gain  Grf  are  fixed  by  the  platform 
and  the  operating  environment.  Examples  1  and  2  assess  the  HF-NVIS  link  budget  using  bandwidth  and 
noise. 

Example  1  (HF-NVIS  and  Bandwidth)  Figure  19  shows  Power-SNR  curves  at  5  MHz  for  three 
bandwidths.  The  link  budget  is  set  by  Figure  12  using  the  59-in  diameter  vertical  loop  antenna  to 
transmit.  Figure  14  sets  the  antenna  gain  as  Ga  =  0  dB  at  5  MHz.  Table  4  sets  the  transducer  power  gains 
Gt  for  each  bandwidth  assuming  best  possible  state-space  matching  limited  to  three  reactive  elements. 
Figure  18  sets  the  urban  noise  Pn,rf  for  each  bandwidth.  The  HF  gain  is  bounded  as  Ghf  =  _100  dB.  If 
the  transmitter  is  limited  to  pt=  100  Watts,  the  50-kHz  bandwidth  is  not  feasible.  The  2  5 -kHz  bandwidth 
forces  the  link  to  operate  with  an  SNR  less  than  5  dB.  Only  the  12-kHz  bandwidth  delivers  an  SNR 
exceeding  10  dB. 


Figure  19.  Power-SNR  tradeoff  for  the  5-MHz  HF-NVIS  link  using  the  59-in  diameter  loop  antenna 
operating  in  urban  noise. 
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Example  2  (HF-NVIS  Link  and  Noise)  Figure  20  shows  power-SNR  curves  at  5  MHz  for  three  noise 
scenarios.  As  in  Example  1,  the  link  budget  is  set  by  Figure  12  using  the  5 9 -in  diameter  vertical  loop 
antenna  to  transmit.  Figure  14  sets  the  antenna  gain  as  G  a  =  0  dB  at  5  MHz .  Table  4  sets  the  transducer 
power  gains  Gt  =  —5.5  dB  for  the  24-kHz  bandwidth  assuming  best  possible  state- space  matching 
limited  to  three  reactive  elements  ( 24-kHz  is  used  for  comparison ).  The  HF  gain  is  bounded  as 
Ghf  =  —100  dB.  Figure  18  sets  the  noise  powers  Pn,rf  for  the  24-kHz  band.  At  transmit  power 
Pt  —  100  Watts,  the  urban  noise  limits  the  SNR  to  5  dB.  However,  suburban  noise  has  an  SNR  exceeding 
15  dB. 


Figure  20.  Power-SNR  tradeoffs  for  the  5-MHz  HF-NVIS  link  using  the  59-in  diameter  loop  antenna 
operating  in  different  noise  models.  The  “X”  is  extrapolated  from  [85],  in  which  the  authors  measure  over 
30-dB  SNR  over  a  3-kHz  bandwidth  using  a  dipole  antenna  on  transmit  that  translates  to  21 -dB  SNR  over 
a  24-kHz  bandwidth. 

Examples  1  and  2  show  that  the  HF-NVIS  antenna  is  the  fundamental  design  parameter  for  this  power- 
only  equation.  Recent  antenna  designs  boosting  antenna  and  matching  gain  are  reviewed  in  Section  2.4. 
Equation  (5)  models  power  only — signal  processing  gain  obtained  by  exploiting  the  HF-NVIS  channel  is 
not  included  in  the  end-to-end  throughput.  For  example,  Section  3  discusses  the  multiple  hops  common  to 
most  HF  channels,  regardless  of  whether  the  channel  is  narrowband  or  wideband.  Receiver  processing 
exploiting  these  multiple  hops  inserts  the  associated  diversity  gain  in  Equation  (5).  Likewise,  wideband  HF 
channels  typically  allow  the  O-  and  X-Mode  paths  to  be  distinguished.  Receiver  processing  exploiting 
these  different  modes  also  inserts  this  diversity  gain  in  Equation  (5). 

The  HF  channels  also  supports  propagation  over  polarization  modes.  If  the  transmitter  is 
dual-polarization  and  the  receiver  is  dual-polarization,  then  2x2  MIMO  gains  are  possible.  The  payoff 
from  this  approach  is  considerable.  Essentially,  there  can  be  up  to  two  channels  rather  than  a  single 
channel.  Thus,  the  12-kHz  channel’s  matching  can  be  exploited  to  deliver  bit  rates  commensurate  with  a 
24-kHz  bandwidth.  However,  antenna  size  on  a  vehicle  limits  this  2x2  MIMO  gain.  Therefore,  Section  2.4 
reviews  progress  in  wideband  HF  communications  emphasizing  mobile  communications  to  guide  the 
analysis  of  the  experimental  measurements. 
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2.4  EVOLUTION  OF  WIDEBAND  HF  AND  HF-NVIS 


This  section  reviews  selected  highlights  in  wideband  HF  communications  emphasizing  wideband  HF 
measurements,  insights  for  wideband  channel  communications,  and  concludes  with  a  review  of  mobile 
HF-NVIS  systems.  The  wideband  HF  literature  goes  back  at  least  five  decades  and  the  progress  by  MITRE 
nicely  brackets  the  topic.  In  1983,  Perry  made  explicit  the  payoff  for  wideband  HF  [60]: 

Wideband  HF  will  offer  an  alternative  to  satellite  and  tropospheric-scatter  communications. 

It  will  permit  the  use  of  low-data-rate  spread  spectrum  signaling,  with  the  attendant  virtues  of 
anti-jam  protection  and  covertness. 

This  paper  cites  the  wideband  HF  measurements  made  in  the  1960s  over  2,000-km  paths  [7]  and  develops 
over-the-air  equalizers  operating  across  1-MHz  bandwidths. 

In  1986,  Saluous  and  Shearman  sought  to  [65]  “determine  the  widest  bandwidth  that  can  be  used  for 
HF  ionospheric  propagation  without  channel  equalization.”  Their  HF-NVIS  measurements  were  collected 
over  a  234-km  link  in  the  United  Kingdom  using  frequency  sweeps  up  to  5  MHz.  These  measurements 
show  coherent  bandwidth  is  limited  by  the  following  effects: 

•  Polarization  interference  between  the  two  magnetoionic  waves  produces  deep  fades,  and 

•  Phase  nonlinearity  with  frequency  increases  the  duration  of  the  pulse  response  and  thereby  limits  the 
rate  of  digital  transmission 

In  1989,  Wagner-Goldstein-Meyers-Bello  measured  HF  channels  using  1-MHz  bandwidths  in  the  mid- 
and  high-latitudes  over  short-  and  long-haul  paths  and  obtained  the  following  fundamental  insight 
regarding  wideband  HF  channel  communications  [81]: 

For  a  simple  specular  channel,  rapid  fading  occurs  as  a  consequence  of  signals  from  more  than 
one  mode  having  different  Doppler  shifts  and  being  received  simultaneously.  The  ability  to 
isolate  and  to  respond  exclusively  to  a  single  specular  return  offers  the  opportunity  of 
eliminating  rapid  fading  for  non-disturbed  channels. 

In  1990,  Katal  summarized  performances  for  the  standard  ground-mounted  HF-NVIS  antennas  [39]. 
The  horizontal  quad  loop  antenna  (quarter  wavelength  per  side;  quarter  wavelength  above  ground)  has  the 
greatest  gain  at  high  take-off  angles.  Katal  also  examined  the  long-wire  antennas  and  reports  that  the 
actual  radiation  pattern  differs  from  the  theoretical  calculation  because  of  ohmic  loss.  In  1992,  Nissen  and 
Bello  measured  HF  channels  at  several  frequencies  using  1-MHz  bands  over  a  3100-km  link  between 
Greenland  and  Massachusetts  [54],  shown  in  Figure  21.  Although  the  channels  linked  the  high  latitudes  to 
the  mid-latitudes,  the  joint  distribution  of  the  measured  delay  spread  and  Doppler  spread  matches  other 
measurements  of  high-latitude  channels  (See  Section  4.5). 

Starting  in  the  early  2000s,  several  countries  began  assessing  HF-MIMO.  All  systems  employ  large 
antennas  (relative  to  a  ground  vehicle)  and  static  ground-based  arrays. 

China:  In  2004,  Xu,  Zhang,  Yang,  and  Wang  simulated  HF-MIMO  over  the  Consultative  Committee  on 
International  Radio  (CCIR)  poor  channel  (2-ms  delay;  1  Hz  Doppler  spread)  showing  that  a  2  x  2  MIMO 
system  delivered  better  performance  over  a  comparable  single  input  single  output  (SISO)  system  [86].  In 
2005,  Zhang,  Yang,  Luo,  and  Xu  again  simulated  HF-MIMO  over  the  CCIR  poor  channel  to  report  that  the 
simulated  2x2  MIMO  system  yields  4-15-dB  diversity  gain  over  a  comparable  SISO  system  [91]. 
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RMS  Doppler  Spread  [Hz] 


Figure  21 .  Left:  delay  spread  versus  Doppler  spread  over  long-haul  wideband  HF  links  (Reprinted  with 
permission  of  the  MITRE  Corporation)  [54].  Right:  geographic  location  of  the  HF  link.  The  lines  other  than 
the  HF  link  are  magnetic  field  lines. 


France:  In  2004,  Perrine,  Erhel,  Lemur,  Bertel,  and  Bourdillon  reported  on  one  of  the  first  over-the-air 
wideband  HF  SIMO  experiments  over  a  300-km  link  in  France  between  Poitiers  and  Monterfil  conducted 
at  the  end  of  year  in  2000  [61].  Figure  22  is  a  schematic  of  the  system  using  a  single  transmit  antenna  and 
multiple  co-located  antennas — four  loops  and  four  dipoles — to  exploit  polarization  diversity. 


8  collocated 


transmitting  amplifier  receiving 

system  system 


Figure  22.  HF-SIMO  link  for  over-the-air  transmissions  recorded  in  2000;  20  kbps  delivered  over  a  6-kHz 
band  (figure  courtesy  of  Annals  of  Geophysics)  [61  ]. 

In  2009,  Ndao  and  associates  begin  an  HF-MIMO  effort  starting  with  simulations  leading  to  over-the- 
air  experiments  over  a  280-km  link  between  Monterfil  and  Lucay  in  April  of  201 1  [51,  52].  Figure  23  is  a 
schematic  of  the  HF-MIMO  system.  Their  2013  paper  summarizes  these  efforts  [53]: 

The  simulation  of  the  global  system,  including  a  realistic  model  of  the  ionospheric  channel, 
indicates  a  significant  capacity  gain  if  compared  to  a  SISO  system.  An  experimental 
validation,  involving  a  280-km  long  radio  link,  confirms  this  prediction  and  underlines  the 
possibility  to  reach  a  data  rate  of  24.09  kbps  within  a  4.2-kHz  wide  bandwidth  with  a  good 
quality  of  service.  This  performance  significantly  exceeds  values  encountered  in  current 
standards. 
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Figure  23.  HF-MIMO  system  for  the  over-the-air  transmissions  recorded  in  201 1 ;  24  kbps  delivered  over  a 
4.2-kHz  band  (figure  courtesy  of  EURASIP  Journal  on  Wireless  Communications  and  Networking)  [53]. 


United  Kingdom:  In  2006,  Strangeway s  estimated  HF-MIMO  channel  capacities  [74],  [75].  From 
September  2007  to  September  2008,  the  University  of  Leicester  and  Durham  University  conducted 
HF-MIMO  channel  soundings  over  the  255-km  path  between  these  two  sites  [1],  [29],  [30],  [22].  Both  the 
transmitter  site  and  the  receiver  site  were  populated  with  heterogeneous  antennas  spread  over  tens  of 
meters.  A  typical  transmit  antenna  was  a  7-m  vertical  whip.  Figure  24  shows  a  Giselle-type  antenna  [45] 
used  for  receive.  Although  the  Giselle-type  antennas  can  be  sized  for  vehicles,  realizing  multiple  HF 
transmitters  from  a  mobile  platform  is  problematic  because  of  the  coupling  between  the  antennas. 

Figure  25  presents  an  upper  bound  of  the  capacity  of  these  HF-MIMO  channels  estimated  from  the 
channel  soundings  and  idealized  assumptions.  Although  the  SISO-to-MIMO  capacity  gains  are  small  at 
low  SNR,  the  reduction  in  SNR  even  at  low  capacity  is  significant.  At  5  bps/Hz,  a  4  x  4  MIMO  system 
needs  an  SNR  of  6  dB  in  comparison  to  the  15  dB  required  by  SISO.  This  SNR  reduction  of  9  dB 
corresponds  to  a  transmitter  power  reduction  by  a  factor  of  8.  Problematic  for  mobile  HF-MIMO  are  the 
array  sizes  and  the  size  of  the  individual  elements  of  the  array. 


Figure  24.  Modified  Giselle  antenna  for  receive  Figure  25.  HF-MIMO  capacity  near  5.255 

(figure  courtesy  of  Dr.  Stuart  Feeney)  [1].  MHz  estimated  from  over-the-air  channel 

soundings  recorded  on  3  September  2008 
[30]  (figure  courtesy  of  Radio  Science). 
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Australia:  The  University  of  Adelaide  and  the  Defense  Science  and  Technology  Group  (DSTG)  have  a 
long  history  of  HF  research,  with  an  emphasis  on  over-the-horizon  HF  radar.  In  2003,  Brine,  Lim,  Massie, 
and  Marwood  produced  idealized  HF-MIMO  channel  capacity  estimates  based  on  HF  soundings  from 
Darwin  to  Adelaide  [10].  In  2005,  Yau  conducted  dual-polarization  channel  soundings  over  a  1,200-km 
path  from  Alice  Springs  to  Lake  Bennet  in  Australia.  Figure  26  shows  a  channel’s  scattering  function  for 
vertical  polarization  at  10.677  MHz  estimated  from  a  46-kHz  band  [88].  The  channel  shows  three  modes 
spread  in  delay  by  2  ms  and  Doppler  spreading  from  —0.5  to  0.8  Hz. 


Channel  Scattering  Function  for  vertical  polarisation  30-Mar-2005  18:54:10 


Doppler  (Hz) 


Figure  26.  Vertical  polarization  scattering  function;  10.677  MHz  (figure  courtesy  of  Dr.  Kin  Shing  Bobby 
Yau)  [88]. 


United  States:  In  October  2012, 

Kuma  Signals,  LLC  conducted 
dual-polarization  HF-NVIS  channel 
soundings  over  a  19.7-km  link  in 
Austin,  TX  at  7.08806  MHz  [15]. 

Both  the  transmitter  and  receiver  used 
two  co-located,  orthogonal,  horizontal 
dipoles  cut  to  85  ft  in  length.  Figure  27 
compares  the  the  ergodic  capacity 
between  the  SISO  and  the  2  x  2  MIMO 
using  the  channel  matrix  estimated  from 
the  soundings.  Kuma  Signals  used  these 
channel  estimates  to  propose  a  MIMO 
waveform,  backwards  compatible 
MIL-STD- 1 88- 1 1 0C  Appendix  D, 
promising  “...a  1 16%  improvement 
in  overall  throughput  and  a  15-dB  SNR 
improvement  for  the  highest-rate  modes 


5  10  15  20  25  30 


MIMO  SNR(dB) 

Figure  27.  Dual-polarization  MIMO  capacity  near  7  MHz 
estimated  from  over-the-air  channel  soundings  recorded  on  19 
October  201 2  (figure  courtesy  of  Kuma  Signals)  [1 5]. 

.[16]” 
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In  December  2013,  MITRE  conducted  wideband  dual-polarization  HF-MIMO  experiments  over  a 
50-km  link  from  Bedford,  MA,  to  Worcester,  MA,  at  5-MHz  using  a  49-kHz  bandwidth.  Figure  28  shows 
the  crossed  dipole  antennas  used  at  both  the  transmitter  and  receiver.  The  system  used  approximately  5  W 
per  channel  to  deliver  94  kbps  using  binary  phase-shift  keying  (BPSK)  and  9  kbps  with  sparse  orthogonal 
frequency  division  multiplexing  (OFDM).  Figure  29  shows  the  bit  error  rate  as  a  function  of  SNR  for 
OFDM  (upper  plot)  and  BPSK  (lower  plot). 


Estimated  BER  vs.  SNR 


Estirrutec  BCR  vs.  SNR 


(C)  Estimated  SNR 


Figure  28.  TX/RX  crossed  dipole  antennas  Figure  29.  OFDM  and  BPSK  waveforms 

(reprinted  with  permission  of  the  MITRE  (reprinted  with  permission  of  the  MITRE 

corporation)  [66].  corporation)  [66]. 

In  2014,  TrellisWare  Technologies,  Inc.  conducted  soundings  over  19  days  across  200  HF-NVIS 
channels  using  24-kHz  bandwidths  sweeping  from  2  to  12  MHz.  The  transmitter  was  located  in  Rancho 
Bernardo,  CA.  Two  synchronized  HF-NVIS  receivers  were  located  in  Encinitas,  CA  (ground  range  ^12 
mi)  and  Fos  Angeles,  CA  (ground  range  ^120  mi).  Synchronization  permits  direct  comparison  of  noise, 
interference,  and  channel  states.  Figure  30  compares  the  noise  and  interference  between  the  two  HF-NVIS 
channels  at  4:30  am  (Pacific  Standard  Time)  PST. 

Both  sites  show  common  interference  near  6  MHz,  7.5  MHz,  9.75  MHz,  and  11.75  MHz.  The  Fos 
Angeles  receiver  measured  wideband  impulsive  noise  (horizontal  lines),  likely  caused  by  lightning,  and 
larger  low  frequency  noise.  In  contrast,  the  Encinitas  receiver  shows  a  higher  noise  floor  above  1 1  MHz. 
These  measurements  support  Bark’s  observation  [5]: 

The  HF  radio  channel  is  one  of  the  most  difficult  channels  for  two  main  reasons:  the 
frequency- selective  fading  caused  by  ionospheric  multipath  propagation,  and  the  severe 
interference  from  other  users  ...  in  many  cases,  the  performance  of  the  HF  systems  is  limited 
rather  by  interference  than  by  bad  propagation  conditions  . . . 
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Figure  30.  Simultaneous  measurement  of  HF-NVIS  noise  and  interference  in  the  early  morning  (figure 
courtesy  of  TrellisWare  Technologies,  Inc.)  [24].  Vertical  axis  is  time  in  seconds. 


Consequently,  any  HF-NVIS  system  must  adapt  to  the  location- specific  noise  and  interference. 
Analysis  of  the  nine  terabytes  of  channel  soundings  showed  that  the  standard  narrowband  models  do  not 
capture  the  wideband  channel  features  [79]: 

traditional  doubly-spread  HF  channel  models  are  . . .  ill-suited  to  describe  complex  NVIS 
propagation  phenomena.  In  particular,  it  is  observed  that  the  NVIS  channel  is  often 
non- stationary  over  a  few  minutes;  with  taps  rapidly  shifting  in  delay,  exhibiting  abrupt 
changes  in  phase  characteristics,  and  periodic  variations  in  magnitude. 

The  limitation  of  this  observation  is  that  the  measurements  are  mid-latitude.  High-latitude  HF  channels 
are  known  to  be  more  complex.  Some  of  the  narrowband  phenomena  are  reviewed  in  Section  4.5. 

Assessment  for  mobile  HF-MIMO:  MIMO  only  makes  sense  when  there  are  multiple  paths.  The 
HF-NVIS  channels  typically  admit  at  least  two  paths — the  O-  and  X-Modes.  Therefore,  dual-polarization 
MIMO  is  possible  provided  the  transmitter  and  receiver  are  equipped  with  antennas  that  match  these 
modes.  Section  3  shows  that  the  O-  and  X-Mode  are  left  and  right  elliptically  polarized.  Therefore, 
equipping  the  transmitter  and  receiver  with  two  orthogonal  antennas  is  one  realization  of  this  particular 
2x2  MIMO.  This  setup  also  avoids  the  fading  caused  by  these  modes  canceling  in  a  single  antenna. 
However,  broadcasting  from  multiple  HF  antennas  on  a  ground  vehicle  is  challenging — the  cross  coupling 
between  antennas  essentially  voids  the  MIMO  gain.  Therefore,  a  SIMO  architecture  is  currently  the  most 
feasible  candidate  for  mobile  HF-NVIS.  An  excellent  analysis  of  the  antenna  selection  for  non-mobile 
HF-SIMO  is  found  in  the  2015  paper  by  Erhel,  Lemur,  Oger,  and  Masson  [21]. 
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Mobile  HF-NVIS:  The  antenna  is  the  bottleneck  for  any  mobile  HF-NVIS  communication  system. 
Filipovic’s  observations  are  particularly  insightful  [33]: 


Any  vehicular  mounted  antenna  will  have  a  quality  factor  far  below  the  fundamental  upper 
limit  for  its  electrical  size.  This  is  because  for  the  [platform]  to  remain  mobile  and  functional, 
an  antenna  cannot  fill  the  entire  space  of  the  . . .  platform. 


One  approach  to  the  mobile  HF-NVIS  antenna  simply  adapts  an  existing  antenna  to  the  vehicle. 
Figure  31  shows  the  canonical  tilted  whip  antenna.  Any  proposed  HF-NVIS  antenna  system  must  deliver 
greater  performance  than  this  existing  antenna  system. 


Figure  31 .  Tilting  a  whip  to  produce  HF-NVIS  [1 9]. 


Several  research  teams  provided  designs  using  the  platform  itself  to  boost  antenna  performance.  In 
2002,  Austin  and  Lui  analyzed  tilted  whips  and  loops  on  a  simulated  generic  vehicle  to  find  that  a 
capacitive-loaded  loop  is  the  preferable  configuration  [4].  In  2010,  Deng  and  associates  proposed  spiral 
antennas  [18].  Figures  32  and  33  show  a  recent  HF/HF-NVIS  antenna  developed  for  ground-mobile 
applications  that  uses  inductive  tuning. 


Figure  32.  Compact  HF-NVIS  antenna  [31] 
(figure  courtesy  of  Hi-Q-MilitaryAntennas, 
Inc.). 
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Figure  33.  Mounted  on  an  HMMWV  [31] 
(Figure  courtesy  of  Hi-Q-MilitaryAntennas, 
Inc.). 


In  2016,  Shih  and  Behdad  applied  characteristic  modes  to  design  mobile  HF-NVIS  antennas  using  the 
platform  as  part  of  the  antenna  system  [71].  Figure  34  shows  the  simulated  HF-NVIS  patterns  generated 
by  low-profile  loop  antennas  mounted  on  a  simplified  model  of  an  Expeditionary  Fighting  Vehicle  (EFV). 
The  vehicle  is  excited  by  low-profile  loop  antennas  mounted  lengthwise  in  the  rear.  The  number  of 
antennas  increases  viewing  the  figure  from  left  to  right.  The  gain  patterns  are  plotted  over  each  antenna 
configuration  and  show  the  NVIS  directionality.  The  current  flows  on  the  EFV  are  excited  by  the  antenna 
configurations.  These  simultaneous  visualization  of  the  current  flows  and  associated  patterns  guide  the 
placement  of  antennas. 


Figure  34.  HF-NVIS  patterns  produced  by  selected  loop  antennas  exciting  characteristic  modes  on  a 
simulated  EFV  [71]  (used  with  the  author’s  permission). 

In  2016,  a  substantial  study  of  vehicle-mounted,  wideband  HF-NVIS  was  completed  by  Filipovic’s 
team  at  the  University  of  Colorado  [33].  Several  types  of  antennas  (e.g.,  half-loop,  inverted-L,  and 
top-loaded  monopole)  were  analyzed  on  several  vehicle  types.  Figure  35  shows  one  double  horizontal  loop 
antenna  mounted  on  an  Oshkosh  Defense®  mine-resistant  ambush  protected  all-terrain  vehicle  (M-ATV). 
This  study  encompassed  the  realized  gain — matching  and  upward  radiation — as  affected  by  realistic 
vehicles  operating  over  various  ground  planes. 

Figure  36  shows  the  effect  of  the  ground  on  these  vehicle-mounted  antennas.  The  resulting  realized 
gain  Gr  encompasses  the  transducer  power  gain  Gt  of  matching  circuit  and  the  antenna  gain  Ga  in  the 
upward  direction: 

Gr  —  Gr  +  Ga  [dB] . 

Figure  36  shows  that  ground  conditions  degrade  performance  upwards  of  5  dB.  Therefore,  any  credible 
link  budget  for  a  mobile  HF-NVIS  system  must  also  include  the  operational  ground. 

2.5  SUMMARY 

This  in-depth  assessment  of  mobile  HF-NVIS  shows  that  the  antenna  is  the  bottleneck  of  this  HF  link. 
MIMO  promises  some  boost  in  throughput,  but  size  limitations  preclude  multiple  transmitters  operating 
from  a  vehicle  platform.  Consequently,  only  SIMO  systems  may  deliver  useful  performance.  Specifically, 
the  SIMO  diversity  gain  could  reduce  bandwidth  requirements.  Reduced  bandwidth  enables  better 
matching  and  lower  power.  Regardless  of  the  MIMO  flavor,  all  these  diversity  schemes  require  knowledge 
of  the  HF-NVIS  channel  and  the  O-  and  X-Modes.  Accordingly,  the  next  section  undertakes  a  review  of 
the  HF-NVIS  channel  emphasizing  SIMO  usage  via  the  O-  and  X-Modes. 
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Figure  35.  One  double  horizontal  loop  antenna  mounted  on  an  M-ATV  [33]. 
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Figure  36.  Realized  gain  for  two  horizontal  loop  antennas  mounted  on  a  simulated  AAV  over  24-kHz 
bandwidths  as  a  function  of  frequency  and  ground  [33]. 


23 


3.  HF-NVIS  CHANNELS  AND  THE  IONOSPHERE 


The  HF-NVIS  channels  are  governed  by  the  ionosphere.  Figure  37  is  a  schematic  showing  some  of  the 
multiple  hops  that  the  HF-NVIS  rays  may  take  when  traveling  from  the  transmitter  to  receiver. 


Figure  37.  HF-NVIS  multiple  hops  [20]. 

Each  hop  induces  a  delay,  attenuation,  fading,  and  frequency  shift  in  the  transmitted  signal  sr(t). 

For  example,  if  the  channel  exhibits  the  rightmost  hops  in  Figure  37,  the  received  signal  sR(t)  should 
have  the  structure 

sR(t)  =  af(t)sT(t  -  Tf)  +  a2p{t)sT{t  -  t2f)  +  -  ^f(Es)f)?  (6) 

where  the  1-Hop  has  delay  tf,  the  2-Hop  has  delay  T2f,  and  the  third  hop  has  delay  tF(es)f-  The 
amplitudes  ap(t)  are  time-varying  and  exhibit  more  attenuation  and  fading  as  the  hop  number  increases. 
The  measured  channels  show  1-,  2-,  2.5-,  3-  and  4-Hops.  Therefore,  the  HF-NVIS  channel  appears  to  be 
completely  determined  by  this  model.  However,  Equation  (6)  implicitly  assumes  that  the  antennas  are 
already  folded  into  the  channel  and  that  polarization  is  absent.  More  to  the  point,  the  ordinary  mode 
(O-Mode)  and  the  extraordinary  mode  (X-Mode)  are  absent  from  this  model.  This  section  reviews  the  O- 
and  X-Modes,  their  presence  in  the  ionograms,  and  concludes  with  a  general  model  encompassing  these 
modes  and  the  antennas  in  the  HF-NVIS  link. 

3.1  O-MODE  AND  X-MODE 

A  high-level  description  of  the  O-  and  X-Modes  is  found  in  Gillies  [27]: 

“. . .  when  a  wave  propagates  through  an  ionized  medium  with  an  external  magnetic  field,  two 
different  indices  of  refraction  exist.  These  two  indices  correspond  to  two  modes  of 
propagation.  The  ordinary  mode  (O-Mode)  and  the  extraordinary  mode  (X-Mode)  are  affected 
differently  by  the  external  magnetic  field  (B-field).” 
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Figure  38  illustrates  the  Earth’s  magnetic  field.  The  magnetic  field  runs  approximately  north-south  in 
the  mid-latitudes.  Waves  traveling  north- to- south  are  less  affected  by  the  magnetic  field  than  waves 
traveling  east-to-west.  Figures  39  and  40  show  simulated  ray  paths  for  HF-NVIS  waves  traveling 
west-to-east.  The  magnetic  field  pushes  the  O-Mode  to  the  north  and  the  X-Mode  to  the  south.  The 
O-Mode  also  travels  higher  but  arrives  sooner  because  of  the  differing  index  of  refraction  along  each  of  the 
paths. 


Figure  38.  Illustration  of  the  Earth’s  magnetic  field. 


Figure  39.  1-Flop  O-  and  X-Modes. 


Figure  40.  Top  view  of  1-Hops. 


Travel  times  for  both  HF-NVIS  rays  can  be  approximated  from  the  associated  ionograms  [47,  page 


51]: 


“An  ionogram  is  a  plot  of  the  apparent  length  of  the  propagation  path  as  a  function  of 
frequency.” 

Figure  41  is  an  annotated  ionogram  showing  the  virtual  heights  where  the  upward  traveling  waves  are 
reflected  back  to  Earth.  This  ionogram  was  recorded  at  Point  Arguello,  CA,  on  10  February  2016  at  09:45 
PST.  At  each  frequency,  the  sounding  first  refracts  off  the  ionosphere  and  travels  back  to  the  station 
producing  the  first  “1-Hop”  layer.  This  sounding  then  reflects  off  the  ground,  travels  upward  a  second  time 
to  refract  off  the  ionosphere,  and  then  travels  back  to  the  station  to  produce  the  “2-Hop”  layer.  The 
difference  in  the  refractive  index  for  the  O-  and  X-Modes  becomes  more  apparent  with  each  hop  and 
increasing  frequency. 
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Figure  41.  lonogram  at  Point  Arguello  on  10  February  2016  at  09:45  PST.  Open  source  ionograms 
courtesy  of  the  Global  Ionospheric  Radio  Observatory  at  the  University  of  Massachusetts,  Lowell  [64]. 


Table  5.  Virtual  heights  and  round-trip  travel 
times  (delay)  at  5.82  MHz  estimated  from 
Figure  41. 

1-Hop  height  (km)  delay  (ms) 


Table  5  lists  the  heights  for  the  1-Hop  modes  and 
their  associated  round-trip  travel  times  by  dividing 
the  virtual  round-trip  distance  by  the  speed  of  light. 

The  “virtual”  heights  are  determined  by  multiplying 
the  actual  time-of-flight  and  the  speed  of  light.  This 
“virtual”  height  model  assumes  the  waves  travel  at  the 
speed  of  light — not  the  actual  physics  where  the  speed 
of  propagation  is  determined  by  the  index  of  refraction 

causing  the  upward  traveling  wave  to  slow  and  then  reverse  as  it  refracts  back  to  Earth.  These  delays  show 
that  the  simple  model  of  Equation  (6)  must  be  extended  to  encompass  the  O-  and  X-Modes. 
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3.2  INDEX  OF  REFRACTION  AND  THE  MAGNETIC  FIELD 

The  O-  and  X-Modes  are  determined  by  the  electron  density  and  the  magnetic  field  of  the  Earth. 
Figure  42  is  an  ionogram  recorded  at  Austin,  TX,  on  10  February  2016  at  12:10  PST  that  displays  the 
plasma  frequency  (black  line)  as  a  function  of  height. 
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Figure  42.  Selected  plasma  frequencies  (MHz)  on  the  horizontal  axis  and  associated  heights  (km)  on  the 
vertical  axis.  Open  source  ionograms  courtesy  of  the  Global  Ionospheric  Radio  Observatory  at  the 
University  of  Massachusetts,  Lowell  [64]. 


The  plasma  frequency  fp  links  to  the  electron  density  Ne  as  [58,  Eq.  6.42]: 


Ne(h)  :=  fp{h)2 


47r2eo  me 


Qe 


[N  •  m  3], 


(7) 


where  Ne  is  the  number  of  electrons  per  cubic  meter,  fp  is  measured  in  Hertz,  and  the  remaining  constants 
are 


eo  8.85418782  •  1CT12  Permittivity  of  free  space  (F/m) 

me  9.10938291  •  10-31  Electron  mass  (kg) 

qe  —1.60217657  •  10-19  Charge  of  an  electron  (C) 

The  electron  density  Ne  and  the  Earth’s  magnetic  field  determine  the  phase  refractive  indices  n± 
governing  the  O-  and  X-Mode  propagation  through  the  ionosphere  [69,  Eq.  5],  [17,  Eq.  3.8]: 

4  :=  1 - ; - ^  j  =^=,  (8) 

1  _  ~  2(1  -X-jZ)  ^  y  +  4(1  -X-jZ)2 

where  the  magnetoionic  frequencies  X,  Y,  and  Z  are  described  in  Equations  (9),  (10),  and  (11). 

Plasma  Frequency:  The  electron  density  determines  the  normalized  plasma  frequency  squared  [89,  Eq.  4.1.5]: 

X=^ Up  =  2nfp,uj  =  2irf  (9) 
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where  the  normalization  is  made  with  respect  to  the  operating  HF  frequency  /. 


Gyro  Frequencies:  The  Earth’s  magnetic  field  causes  the  electrons  to  rotate  around  the  magnetic  lines  of 
force.  The  evolution  of  the  polarization  along  the  ray  path  is  governed  by  this  rotation  or  gyro-frequency 
[89,  Eq.  4.4.3] 

ujh  =  ~~  Bo  [rad  •  sec-1], 
me 


where  Bo  is  the  magnitude  of  the  Earth’s  magnetic  induction  field  vector  Bo  and  determines  the 
normalized  gyro-frequency  [89,  Eq.  4.4.3] 

Y=—. 
c u 


(10) 


The  direction  of  the  Earth’s  magnetic  field  with  respect  to  the  ray  path  appears  in  the  refractive  index 
via  the  longitudinal  and  transverse  components  of  the  normalized  gyro-frequency  [89,  Eq.  4.5.8] 
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where  Oh  is  the  angle  between  the  magnetic  field  and  the  wave  vector. 


Electron  Collision  Frequency:  The  collisions  of  the  electrons  with  ions  and  neutral  molecules  determines 
the  local  absorption  along  the  ray  path  and  the  subsequent  total  absorption.  Absorption  appears  in  the 
index  of  refraction  n  via  the  normalized  electron  collision  frequency  u ue- 


LUe 

UJ 


(ii) 


The  electron  collision  frequency  is  the  sum  of  electron-ion  collision  frequency  uje,i  and  the  electron- 
neutral  collision  frequency  oje,n  [90]: 

ujej  =  10-6  x  Ne  x  T-3/2  1 59  +  4.18  log  j  [rad -sec-1], 

c oe,n  —  5.4  •  10-16  x  Nn  x  Tg1/2  [rad  •  sec-1], 

where  Te  is  the  electron  temperature  and  Nn  is  the  neutral  particle  density.  . 

Channel  Modeling:  The  electron  density  profile  and  the  Earth’s  magnetic  field  govern  the  O-  and 
X-Modes  via  the  phase  refractive  index  n±.  Extracting  the  time-of-flight  and  propagation  losses  along 
these  modes  is  the  basic  channel  modeling.  The  complications  are  that  the  antennas  must  couple  into  these 
modes  and  that  the  polarization  evolves  along  these  ray  paths.  The  next  two  sections  examine  these  effects 
for  1-Hop  and  2-Hop  channels. 

3.3  THE  1-HOP  CHANNELS 

Figure  37  shows  only  1-Hop  channels  whereas  the  O-  and  X-Mode  propagation  produce  two  distinct 
channels.  As  remarked  in  the  introduction  to  this  section,  “channels”  implies  that  the  transmitter  antenna 
launches  an  electromagnetic  (EM)  wave  into  each  mode  and  the  receiver  antenna  intercepts  the  EM  wave 
from  each  mode.  Vatberg  and  Lundborg  offer  an  excellent  description  of  this  O-  and  X-Mode  propagation 
[80]: 
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When  the  wave  first  enters  the  ionosphere,  it  will  split  into  two  wave  modes,  ordinary  (O)  and 
extraordinary  (X).  If  the  polarization  of  the  incident  wave  matches  one  of  the  modes  exactly, 
all  the  power  will  go  into  that  wave  and  the  other  mode  will  be  suppressed,  but  generally,  both 
modes  will  be  generated,  and  these  will  propagate  independently. 

Figure  43  illustrates  the  simplest  propagation  along  the  modes,  the  mode  “shedding”  [27],  and  the 
interaction  of  both  antennas  at  the  transmitter  and  receiver  with  these  modes.  The  left  and  right  elliptical 
polarizations  are  oriented  in  the  right-handed  coordinate  system  with  propagation  into  the  page  marked  by 
the  “0”  [58,  Figure  6.3].  The  descriptor  “simplest”  means  the  polarization  orientation  does  not  change 
along  the  modes.  However,  Booker’s  1934  paper  derived  that  orientation  of  the  polarization  can  change 
along  the  modes  [9]. 
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Figure  43.  Block  diagram  of  1-Hop,  O-  and  X-Mode  propagation  without  polarization  switching. 


Figure  44  illustrates  one  possibility  where  the  polarization  along  both  modes  swings  through  the  linear 
polarization  and  reverses.  Indeed,  Booker’s  1934  paper  shows  polarization  switching  right-left-right  on 
one  path  and  left-right  on  the  other  path. 
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Figure  44.  Block  diagram  of  1-Hop,  O-  and  X-Mode  propagation  showing  polarizations  switching. 


Appleton  established  the  nomenclature  of  the  O-  and  X-Modes  and  their  polarizations  [3]: 

The  polarizations  of  the  magnetic  forces  in  the  two  waves  are  similar  ellipses  the  major  axes 
of  which  are  at  right  angles.  For  a  direction  of  propagation  making  an  acute  angle  with  the 
positive  direction  of  the  magnetic  field,  the  polarization  of  the  ordinary  wave  is  left-handed 
and  that  of  the  extraordinary  wave  right-handed. 
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Table  6  links  the  ±  sign  in  the  refractive  index  of 
Equation  (8)  the  O-  and  X-Modes,  and  their 
polarizations  under  Appleton’s  acute  angle  constraint. 

The  right-  and  left-elliptical  polarizations  determined 
by  the  O-  and  X-Modes  are  necessary  to  calculate  the 
voltage  response  at  the  receive  antenna  generated  by 
the  down-coming  waves.  Likewise,  the  polarization 
of  the  transmit  antenna  that  launches  the  up-going  wave  determines  how  the  electric  field  will  “split  when 
the  wave  enters  the  ionosphere.”  Thus,  interaction  of  the  transmit  and  receive  antennas  with  the  O-  and  X- 
Modes  are  critical  to  the  end-to-end  HF-NVIS  channel  model.  Consequently,  Equation  (6)  should  generalize 
as  a  coherent  sum  of  the  O-  and  X-channels  [801,  [70]: 


Table  6.  Index  of  refraction,  modes,  and 
acute-angle  polarizations. 

n±  Mode  Polarization 

+  O-Mode  Left  Elliptical  (LEP) 

-  X-Mode  Right  Elliptical  (REP) 


SR(t)  =  SRt0(t)  +  SRjx(t). 


The  “O”  and  “X”  label  the  O-Mode  or  X-Mode  channels  measured  at  the  receiver.  In  terms  of  the 
experimental  measurements,  the  O-Mode  channel  is  the  LEP  component  of  the  received  electric  field — 
provided  the  polarization  has  not  switched  on  this  1-Hop  channel.  Likewise,  the  X-Mode  channel  is  the 
REP  component  of  the  received  electric  field — provided  the  polarization  has  not  switched  on  this  1-Hop 
channel.  Consequently,  if  the  HF-NVIS  channel  is  modeled  as  the  1-Hop  of  Figure  43,  the  O-  and  X-Mode 
responses  at  the  receiver  are  modeled  as 

SR,o(t)  =  af^o^srit  -  rF,o), 

sR,x(t)  =  Gf,x(£)st(£  “  rF,x)- 

The  “F”  refers  to  the  F-layer  refection  and  the  “O”  and  “X”  label  the  O-Mode  or  X-Mode  propagation. 

The  channel  fading  coefficients  aF,o(^)  and  ax,o  (t)  encompass  the  transmit  and  receive  antennas  feeding 
into  the  O-  and  X-Mode  channels. 

3.4  THE  2-HOP  CHANNELS 

Including  multiple  hops  requires  that  the  channel  model  include  ground  reflections.  Ground  reflections 
typically  change  the  polarization  on  the  multiple  hops.  Figure  45  is  adapted  from  the  cover  of  Balanis’ 
book  [6].  The  figure  shows  an  incident  left  circular  polarized  (LCP)  planewave  reflecting  as  a  right  circular 
polarized  (RCP)  planewave. 
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Figure  46  is  a  block  diagram  showing  the  propagation  for  two  hops.  The  figure  shows  that  the 
reflection  from  the  ground  bounce  causes  a  mixing  of  the  polarizations. 


Ground 

Reflection 


O-Mode 


Ionosphere 


X-Mode 


Ionosphere 


Figure  46.  Block  diagram  of  the  2-Hop  O-  and  X-Mode  channels  showing  ground  scattering  between 
channels. 


In  the  extreme  case  where  the  O-  and  X-Modes  are  circular,  the  plane-wave  components  of  the  incident 
electric  field 
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The  conversion  between  circular  and  rectangular  field  components  is  [50,  Eq.  6.119] 
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If  the  ground  scattering  matrix  in  rectangular  coordinates  is 
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where  E rH  and  Ey  are  the  horizontal  and  vertical  components  of  the  reflected  field,  ElH  and  Elv  are  the 
horizontal  and  vertical  components  of  the  incident  fields,  and  pn  and  py  are  available  from  measurements, 
the  ground  scattering  matrix  in  circular  coordinates  is 
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The  following  examples  show  two  extremes  of  the  ground  scattering  matrix.  The  first  example  shows 
that  sea  water  almost  completely  “flips”  the  mode  with  little  loss.  The  second  example  shows  that  dry 
ground  still  flips  the  modes  but  also  absorbs  power. 
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Example  3  (Scattering  from  Sea  Water  Surface)  Figures  47  and  48  show  the  horizontal  and  vertical 
reflection  coefficients  for  sea  water  as  a  function  of  the  grazing  angle  or  elevation  from  horizon.  The  plots 
show  for  near  vertical  incident  angles  (grazing  angle  near  90°),  pn  ~  —  1  and  py  ~  1.  Consequently,  the 
circular  scattering  matrix  for  sea  water 
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reflects  the  incident  LCP  as  RCP  and  conversely.  Thus,  2-Hop  channels  reflecting  off  sea  water  support 
only  the  O-to-X  andX-to-0  channels,  unless  the  modes  flip  in  the  ionosphere  as  described  in  Section  3.3. 
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Figure  47.  Horizontal  reflection  coefficient  for  sea 
water. 
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Figure  48.  Vertical  reflection  coefficient  for  sea  water. 
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Example  4  (Scattering  from  very  dry  ground)  Figures  49  and  50  show  the  horizontal  and  vertical 
reflection  coefficients  for  “very  dry  ground”  as  a  function  of  the  grazing  angle  or  elevation  from  horizon. 
The  plots  show  for  near  vertical  incident  angles  (grazing  angle  near  90°),  pH  ~  —  1/2  and  py  ~  1/2. 
Consequently,  the  circular  scattering  matrix  for  very  dry  ground 
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reflects  most  of  the  attenuated  incident  LCP  as  RCP  and  conversely.  Thus,  2-Hop  channels  reflecting  off 
very  dry  ground  support  only  the  O-to-X  andX-to-0  channels,  unless  the  modes  flip  in  the  ionosphere  as 
described  in  Section  3.3. 
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Figure  49.  Horizontal  reflection  coefficient  for  “very 
dry  ground.” 
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Figure  50.  Vertical  reflection  coefficient  for  “very  dry 
ground.” 
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Figure  51  tabulates  the  delays  modeled  by  the  2-Hop  channels.  At  launch,  the  first  hop  enters  the 
ionosphere  either  left  or  right  circular  polarization.  Upon  exiting  the  ionosphere  from  either  the  O-Mode  or 
X-Mode  first  hop,  the  polarization  may  be  reversed.  The  ground  reflection  strongly  reverses  this 
polarization.  The  polarization  mismatch  further  filters  the  reflected  wave  to  match  the  mode  for  the  second 
hop.  To  model  the  2-Hop  time  delays,  all  that  matters  is  tracking  the  modes  that  the  rays  use  to  travel 
through  the  ionosphere. 


Ground 

Reflection 


T0 


L 


Ionosphere  L 
O-Mode 


Ionosphere  _ ^  xn+xn=  xnn 

O-Mode  o  o  oo 


> 


Ionosphere 

X-Mode 


R  Ionosphere  T 7 
X-Mode 

R 


Figure  51 .  Delays  caused  by  polarization  switching  and  ground  reflection. 

If  the  HF-NVIS  channel  is  modeled  as  the  1-Hop  of  Figure  43  and  the  2-Hop  of  Figure  46,  the  the 
received  O-  and  X-Mode  responses  are  modeled  as 

SR,0(t)  ~  CLf,0(t)sT(t  ~  Tf:o)  +  CLf^ox(t)sT(t  -  TF, Ox)  +  «F,00 (t)sT(t  -  TF, Oo), 

sR,x(t)  «  CLf,x(t)sT(t  ~  TF,x)  +  ttF5XO (t)sr(t  ~  TF  xo)  +  ,Xx(t)sT(t  ~  ^"F,Xx), 

using  the  convention  of  Figure  5 1 . 

3.5  SUMMARY 

This  report  finds  that  the  polarization  reverses  during  the  hops.  Consequently,  O-  and  X-Mode 
responses  encompassing  both  1-Hop  and  2-Hop — specific  to  this  polarization  switching — are 
approximated  as  follows: 

sr, oOO  ~  <^f,o (t)sr(t  -  tf? o)  +  cLf,oo{t)sT(t  ~  ^,00), 

SR,x(t)  «  CLf ,x(t)sT(t  -  Tf  x)  +  CLf:Xx(t)sT(t  -  Tp^x)- 

These  responses  are  generally  sufficient  to  model  the  measured  HF-NVIS  channels,  although  ground 
waves  and  weak  3-Hop  modes  are  also  observed. 
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4.  HF-NVIS  CHANNEL  MODELS 


This  section  reviews  two  approaches  for  modeling  the  HF-NVIS  channel  functions.  The  first  approach 
uses  specific  channel  models  simulating  transmissions  for  waveform  assessments.  These  models  specify 
channel  delays  and  Doppler  shifts  that  spread  the  transmitted  signal  in  time  and  frequency.  If  sr(t) 
denotes  the  transmitted  signal,  the  received  signal  sr(£)  is  modeled  as 

K 

sn{t)  =  y]  ake+j27TfktsT(t  -  Tfc), 
k= 1 

where  the  delays  r^’s,  Doppler  shifts  /&’ s,  and  amplitudes  s  are  specified  as  either  a  few  deterministic 
values  or  as  statistical  ensembles  modeling  more  complex  multipath  propagation.  The  second  approach  is 
based  on  the  statistics  of  the  delay  {r^}  and  Doppler  {/*.}  ensembles.  Figure  52  shows  how  these 
HF-NVIS  ensembles  jointly  layout  in  the  delay  Doppler  plane.  Measured  HF-NVIS  channels  determine 
two  distinct  clusters.  The  mid-latitude  channels  spread  delays  up  to  8  ms  while  limited  to  1-Hz  spreads  in 
frequency.  The  high-latitude  channels  spread  Doppler  up  to  16  Hz  while  limited  to  3  ms  in  delay  spreads. 


Fieur*  3-1:  Survey  of  reported  HF-NVIS  channel  measurements 


Figure  52.  Narrowband  HF-NVIS  channel  classification — delay  spread  versus  Doppler  spread  (figure 
courtesy  of  TrellisWare  Technologies,  Inc.  [41]). 

This  section  surveys  selected  narrowband  HF-NVIS  channels  using  the  delay  Doppler  spreads  of  Figure 
52.  Section  4.1  reviews  time-varying  channels  necessary  to  describe  the  delay  spread,  the  Doppler  spread, 
and  the  scattering  function.  Channel  simulation  and  emulation  requires  models  of  the  fading  on  the  ray 
paths,  the  associated  delays  and  Doppler.  Section  4.2  reviews  propagation  losses  along  the  ray  paths,  their 
connection  to  the  magneto-ionic  frequencies,  and  the  ionogram.  Section  4.3  reviews  the  Doppler  shifts  and 
spreads  and  points  out  the  fine  sampling  of  the  electron  density  required  for  modeling  Doppler.  Section  4.4 
tabulates  selected  narrowband  mid-latitude  HF-NVIS  channels  and  the  associated  channel  parameters  for 
HF  channel  emulation.  Section  4.5  describes  selected  narrowband  high-latitude  HF-NVIS  channels.  One 
payoff  of  this  report  is  extending  Figure  52  with  the  wideband  mid-latitude  measurements.  Section  4.6 
remarks  on  the  gaps  in  the  measured  HF-NVIS  channels  and  the  HF-NVIS  channel  models. 
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4.1  GENERAL  CHANNEL  MODELS 


The  fundamentals  of  time-varying  channels  were  laid  out  by  Bello  [8]  and  Kailath  [38]  in  the  early 
1960s.  Mathematically,  a  time- varying  channel  h(t ,  r)  maps  the  transmitted  signal  sr(t)  to  the  received 
signal  «§#(£)  as  [8,  Eq.  8] 

/oo 

h(t,T)sT(t  -  r)dr.  (12) 

-OO 


This  delay-spread  function  h{t ,  r)  models 
the  channel  as  a  collection  of  impulse  responses 
indexed  by  time.  Typically,  t  is  the  “slow 
time”  and  the  delay  r  is  the  “fast”  time.  The 
delay-spread  function  /i(£,  r)  is  part  of  Bello’s 
time-varying  channel  framework  illustrated  in 
Figure  53. 

Taking  the  Fourier  transform  in  delay  produces 
the  Time-variant  Transfer  function  [8,  Eq.  18]: 

/oo 

e-j27 ^Th{t^  r)dr 
-oo 

that  models  the  channel  as  a  collection  of  transfer 
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Figure  53.  Bello’s  channel  functions. 


functions  indexed  by  time  [8,  Eq.  19].  Taking  the 

Fourier  transform  in  time  produces  the  Delay  Doppler  function  [8,  Eq.  18]: 


/oo 

e+j27 r)dt. 

-OO 


The  covariance  of  the  delay  Doppler  function — the  scattering  function — has  emerged  as  a  workhorse  for 
stochastic  channels.  The  delay  Doppler  function  has  covariance 


Ruu(f,  r;  /',  r')  =  E[U(f,  r)U(f\  t')} . 


Under  the  wide-sense  stationary ,  uncorrelated  scattering  assumptions, 

Ruuif,  t;  /',  r')  =  P(f,  t)S(t  -  r')8{f  -  /'),  (13) 

where  P(r,  /)  is  called  the  scattering  function.  The  practical  value  of  the  scattering  function  is  succinctly 
described  [81]: 


Channel  characteristics  such  as  the  extent  of  dispersive  and  multipath  delay  spread  and  the 
extent  of  Doppler  shift  and  spread  are  concisely  represented  by  the  channel  scattering 
function.  The  character  of  the  return,  specular  or  diffuse,  and  its  multimodal  or  multipath 
composition  are  also  usually  evident  from  an  examination  of  the  channel  scattering  function. 

Figure  54  shows  an  estimate  of  the  scattering  function  for  a  long-haul  wideband  HF  channel  from 
Greenland  to  Massachusetts  (See  Figure  21).  This  channel  has  multiple  modes  spread  in  between  100  to 
200  /is  in  delay  and  spread  between  5  to  10  Hz  in  Doppler. 

Theoretical  channel  models  show  that  the  scattering  function  measures  the  joint  distribution  delay  and 
Doppler.  Two  common  time-varying  channel  models  are  the  quadrature  modulation  fading  simulator 
(QMFS)  and  the  phase  modulation  fading  simulator  (PMFS). 
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Figure  54.  Scattering  function  estimated  from  3100-km  link  operating  at  15.5  MHz  on  1  April  1992  at 
07:35  UTC  (reprinted  with  permission  of  the  MITRE  Corporation)  [54], 

Example  5  (QMFS)  [56],  [12],  [63,  pages  697-703]:  The  QMFS  is  the  stochastic  channel  model 

K 

h{t,  t)  =  ^2  ak{t)5(T  -  Tk), 
k=  1 

where  the  {a^(t)}}s  are  complex-valued,  jointly  wide-sense  stationary  (JWSS),  uncorrelated  Gaussian 
random  processes  with  mean  and  spectral  shape  determined  by  the  propagation  environment. 

Remarks: 

L  The  QMFS  maps  the  transmitted  signal  sr(t)  to  the  received  signal  sn(t)  as  faded  and  delayed 
copies  of  sr(t) 

/oo  K 

h(t ,  r)sT(t  -  t )dr  =  ^  ak(t)sT(t  -  rfc). 

-°°  k= l 

2.  Physically,  this  model  “leads  to  a  physical  picture  of  the  channel  as  a  continuum  of  non-moving 
scintillating  scatterers”  [8]. 

3.  More  general  models  admit  time- varying  delays  rn{t). 

4.  Most  RF  channel  emulators  model  an  RF  channel  via  QMFS  where  each  path  is  parameterized  by 
delay  and  fading. 

Example  6  (PMFS)  [ 32 ]:  The  PMFS  is  the  stochastic  channel  model 

1  K 

hit,  t)  =  —=  ^2  ake+j2nfkt5iT  -  Tk).  (14) 

^  71=1 

This  model  is  governed  by  the  following  assumptions: 
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•  PMFS-1:  The  (a&,  /&,  Tfe)  s  are  independent  and  identically-distributed. 

•  PMFS-2:  The  a k ’s  have  zero  mean. 

•  PMFS-3:  The  ak ’s  are  independent  from  the  fk ’s  and  the  Tk ’s. 


Remarks: 

1.  The  PMFS  maps  the  transmitted  signal  srif)  to  the  received  signal  «§#(£)  as  shifted  and  delayed 
copies  of  srif)  : 


/OO  1 

h(t,T)sT(t  -  t)cLt  =  —=y2ake+j2nfktsT{t  ~  Tk). 

-°o  VK  k=1 

2.  PMFS-1  makes  explicit  that  there  exists  a  density  pajfT(a,  /,  r).  PMFS-2  “uncorrelates”  each  path. 
PMFS-3  separates  the  density  paifl)  °f  the  a&’s  from  the  Doppler  and  delay: 


Pa,/,r(o,  /,  T)  =  pa(a)pf,T(f ,  t). 

3.  The  Delay-Doppler  function  [/(/,  r)  takes  the  elegant  form  showing  how  the  signal  “spreads”  in 
time  and  frequency: 

1  K 

U(f,r )  =  —y2ak8(f-fn)8(T-rk). 

k= l 

4.  Assumptions  PMFS-1,  PMFS-2,  and  PMFS-3  force  E[U (/,  r)]  =  0  and  covariance: 

Ruu(f,  r;  /',  t')  =  E[U(f,  t)C/(/',t')] 

A'  a: 

=  £  E[afeafc/]£;[5(/  -  fk)8(j  -  Tk)S(f  -  fk>)8(r'  -  t¥)\ 

k= 1  /c'=l 

=  ^53  aa5(/'  -  /M(r'  -  T)pf,r(f,  T ) 

n=l 

=  CrlPf,r(f,T)8{f  -  f')8(T  -V).  (15) 

Thus,  the  PMFS  gives  a  straight-forward  derivation  that  the  scattering  function  is  proportional  to  the 
probability  distribution  of  the  Doppler  and  delays.  The  marginal  distributions  determine  the  delay  and 
Doppler  spreads  of  Figure  52. 

Doppler  Spread:  The  Doppler  spread  has  probability  density 

POO 

Pf(f)=  Pr,f(rJ)dT. 

J  0 

Typically,  the  mean  is  the  Doppler  shift 

/OO 

fpfU)df 

-oo 

and  the  variance 

/oo 

(/  -  fD)2Pf(f)df 

-OO 
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determines  the  Doppler  spread  or  frequency  spread  [34] 

A  f  —  2  x  <jf 

that  is  the  “2  x  cr”  label  in  the  Doppler- spread  axis  of  Figure  52. 

Delay  Spread:  The  delay  spread  has  probability  density 

/oo 

Prj(r,f)df. 

-OO 

The  cumulative  distribution 

Ft(t)  =  /  pT(u)du 
Jo 

determines  the  “99%”  used  in  the  delay-spread  axis  of  Figure  52.  The  following  examples  specialize  these 
general  channel  models  to  the  HF-NVIS  channels. 

Example  7  (WBHF-NVIS  1-Hop)  The  wideband  HF-NVIS  channel  model  based  on  the  O-  and  X-Modes 
from  Section  3  has  the  form 

h(t,  t )  =  h0(t,  t )  +  hx(t,  t ) 

where  the  delay-spread  functions  for  the  O-  and  X-Modes  are 

ho(t,r)  =  aFt0(t)S(T  -  rF,o) 
hx(t,T )  =  aF}x(t)S(T  -  rFjX). 

This  1-Hop  model  describes  most  of  the  HF-NVIS  channels  measured  for  this  report.  There  are  2-Hop 
and  3-Hop  observed  in  the  measured  HF-NVIS  channels.  However,  the  dry  ground  absorbs  3  to  7  dB  of 
the  power  on  reflection.  Consequently,  2-Hop  amplitudes  are  typically  10  to  20  dB  below  the  first  hop.  In 
addition,  the  O-  and  X-Modes  appear  to  switch  polarizations  (See  Figure  51).  Therefore,  these  1-  and  2- 
Hop  channels  are  modeled  as  follows. 

Example  8  (WBHF-NVIS  2-Hop  Channels  With  Polarization-Switching)  The  wideband  HF-NVIS 
channel  model  based  on  the  O-  and  X-Modes  from  Section  3.4  has  the  form 

h(t,r)  =  h0(t,T) +  hx(t,r)  (16) 

where  the  delay -spread  functions  for  the  O-  and  X-Modes  are 

ho(t,r )  =  aF}0(t)8(T  -  tF}0)  +  aFj0o(t)S(T  -  Tf^oo) 
hx(t,  t )  =  aF}x(t)8(T  -  rF;x)  +  aFtxx(t)S(T  -  rF;Xx), 

assuming  the  polarizations  switch  in  the  ionosphere  on  both  modes. 

The  bulk  of  the  HF-NVIS  channel  literature  is  narrowband.  The  standard  3-kHz  bandwidth  limits  path 
resolution  to  0.33  ms.  The  forthcoming  HF-NVIS  measurements  and  simulations  report  the  O-  and 
X-Modes  on  a  given  hop  are  separated  by  at  most  0.15  ms.  Consequently,  narrowband  HF-NVIS  channels 
must  lump  both  modes  into  a  single  path.  This  narrowband  fading  differs  from  the  wideband  fading 
because  the  O-  and  X-Mode  interference  causes  deep  fades  when  the  receiver  cannot  separate  the  two 
modes.  That  is,  the  narrowband  receiver  may  lose  information  in  the  deep  fades  whereas  a  wideband 
system  would  not. 
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(17) 


Example  9  (NBHF-NVIS)  The  narrowband  HF-NVIS  channel  model  has  the  form 

h(t,  t )  =  ai(t)S(T  -  tl)  +  a2(t)8(T  -  r2)  +  a3(i)5(r  -  r3), 

where  the  fading  processes  o  i-  (  t  )  encompass  the  antennas  at  both  the  transmitter  and  receiver  and  lumps 
the  O-  and  X-Modes  into  a  single  path. 

The  fading  processes  for  HF-NVIS  channels  are  typically  parameterized  by  their  mean  amplitude,  the 
delays,  the  fading  model  or  Doppler  spread,  and  the  Doppler  shift. 

4.2  PROPAGATION  LOSSES  IN  HF  AND  HF-NVIS 

Discounting  the  transmit  and  receive  antennas  coupling  into  the  O-  and  X-Mode  ray  paths  and  the  loss 
caused  by  ground  reflections  on  multiple  hops,  there  are  at  least  three  propagation  effects  that  reduce  the 
transmitted  power  along  every  ray  path  [46,  page  25],  [47,  Section  6.8]: 

•  Spherical  Spreading 

•  Mode  Shedding 

•  Absorption 

Equation  (3),  repeated  here  for  convenience,  models  these  effects  as 

Ghf  =  ^SPREAD  X  ^ABSORB  X  Gshed- 

Spreading  is  approximated  by  the  standard  free-space  loss  caused  by  spherical  propagation  [47,  Section 

6.8]: 

^spread  = 

where  s  is  the  length  of  the  ray  path  and  typically  ranges  from  —90  to  —110  dB  for  HF-NVIS.  The  length 
of  the  ray  path  is  a  relatively  stable  computation.  Therefore,  the  spreading  loss  is  considered  a  stable 
computation. 

Mode  shedding  caused  by  the  polarization  evolving  along  each  ray  path  and  is  in  Sections  3.3  and  3.4  . 
The  3-D  ray  trace  code  estimates  mode  shedding  to  range  from  —3  to  —2  dB  for  this  simulation.  This 
polarization  is  governed  by  the  Earth’s  magnetic  field — relatively  stable  and  well-known — so  tracking 
polarization  along  the  ray  path  is  considered  a  stable  computation.  Thus,  mode-shedding  estimate  is 
considered  a  stable  computation. 

Absorption  is  the  loss  caused  by  the  electron  collisions  with  the  heavy  neutral  molecules  and  with  the 
ions  in  the  ionosphere  [25].  Absorption  appears  in  the  refractive  index  n  via  the  normalized 
electron-collision  frequency  X  (Equation  11).  Although  absorption  is  reported  to  range  from  —20  to  0  dB 
for  HF  propagation  in  general  [47,  Section  4.8],  absorption  modeling  shortfalls  are  significant  and 
absorption  research  is  on-going  [59].  If  the  refractive  index  n  has  real  and  imaginary  parts  n  =  /i  +  j\  the 
attenuation  along  the  path  is  modeled  as  [1 1,  Eq.  14.35] 

Gabsorb  =  exp  kco  J  j^dt 

where  k  is  the  wave  number  and  the  integral  over  time  t  is  along  the  ray  path.  Implicit  in  this  integral  is 
that  the  index  of  refraction  takes  on  four  values — the  Booker  quartic — that  requires  knowing  the  upward 
and  downward  directions  of  propagation.  The  Hamiltonians  used  in  the  3-D  ray  tracing  avoid  the  Booker 
quartic  but  have  “computational  difficulties”  [36].  Consequently,  semi-empirical  absorption  models  are 
still  under  development  [59],  [69].  Therefore,  the  3-D  ray-tracing  code  will  estimate  the  absorption  from 
measurements  by  modeling  the  other  losses  and  compare  to  the  absorption  computed  by  the  Hamiltonians 
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4.3  DOPPLER  SHIFT  AND  DOPPLER  SPREAD  IN  HF 


Doppler  spread  and  Doppler  shift  are  caused  by  the  time-varying  ionosphere,  vehicle  motion,  and  the 
variation  of  the  ionosphere  with  respect  to  frequency.  The  access  to  ionograms  limit  Doppler  modeling 
and  comparison  with  measurements.  Therefore,  this  section  reviews  the  Doppler  mechanisms  to  put  the 
Doppler  measurements  in  context.  In  1994,  Malachias  measured  both  the  Doppler  spread  and  Doppler 
shift  in  HF  signals  and  describes  two  sources  of  Doppler  shift  caused  by  the  ionosphere  [44]: 

If  the  layer  descends,  ...  the  virtual  height  descends  a  proportional  amount  and  the  group  path 
. . .  from  transmitter  to  receiver  becomes  shorter.  The  Doppler  shift  is  thus  proportional  to  the 
path  length  change.  Analogous  to  this  situation,  is  when  the  electron  density  beneath  the  layer 
increases,  lowering  the  reflection  height,  with  a  similar  decrease  in  the  group  path.  With 
respect  to  the  propagating  wave,  the  two  actions  are  indistinguishable,  and  a  Doppler  shift 
likewise  occurs  for  a  change  in  electron  density  of  the  layer. 

This  description  makes  explicit  the  Doppler  caused  by  the  changing  geometry  and  changes  in  the  electron 
density.  Recall  that  electron  density  is  proportional  to  the  square  of  the  plasma  frequency  (Equation  7). 
Figure  55  plots  the  plasma  frequency  fp(h ,  t)  as  a  function  of  height  and  time  The  layer  height  h(fc,  t) 
is  determined  by  fc  =  fp(h(fc ,  t),t)  during  the  time  of  the  channel  sounding.  The  15-minute  sampling 
of  the  ionograms  causes  the  discrete  appearance.  Roughly  speaking,  the  layer  height  decreases  during  the 
interval  08:00  to  1 1:00  PST  remains  relatively  constant  to  14:00,  and  then  begins  to  increase.  There  is  an 
exceptional  event  at  1 1:00  PST. 


Point  Arguello:  Plasma  Frequency 
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Figure  55.  Plasma  frequency  fP(t ,  h)  at  Point  Arguello;  missing  sample  at  09:30  PST. 
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Figure  56  shows  one  estimate  of  the  Doppler  shift  caused  by  the  layer-movement  only.  The  associated 
Doppler  shift  is 

=  ac  =  ^, 

where  this  differential  is  approximated  by  computing  the  numerical  derivative  of  a  spline  fit  to  the  height. 
Consequently,  the  15 -minute  sampling  implies  that  this  Doppler  estimate  is  only  a  rough  estimate  of  the 
Doppler  shift. 


This  estimate  omits  the  change  in  the  electron  density  as  a  function  of  time  and  frequency.  Indeed,  the 
infinitesimal  Doppler  shift  along  the  ray  path  given  in  Equation  D-l  requires  both  dNe/dt  and  dNe/du. 
Therefore,  time-varying  changes  in  the  electron  density  will  not  be  modeled  in  this  report — only  the 
measured  Doppler  shift  will  be  reported. 


Moving-Layer  Dopper  shift:  f  = 5.82  MHz 


Figure  56.  Doppler  shift  caused  by  moving  layer  height  for  fP(t,  h)  =  5.82  MHz  at  Point  Arguello. 


Doppler  spread  is  an  ensemble  of  the  Doppler  shifts.  Referring  to  Equation  (14),  with  a  single  delay, 
such  as  the  O-  and  X-Mode  ray  path  observed  in  the  over-the-air  measurements,  the  delay-spread  function 
takes  the  form 

h(t,  t )  =  a0S(t  -  r0)  exp(</>(f)) 

where  the  instantaneous  phase  4>(t)  is  the  integral  of  the  instantaneous  frequency  { J)  (/)]■: 


<f>(t)  =  (f)0  +  2i r 


If  the  instantaneous  frequency  is  narrowband,  the  Doppler  shift  is 


fD  =  E[fi(t)\ 


and  the  power  spectrum  of  the  phase  modulation 

<b(f)  :=  exp (0(f)) 
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is  proportional  to  the  probability  of  the  Doppler  arising  in  Equation  (15) 

0 CP/(/) 

by  Woodward’s  Theorem  [57].  Consequently,  an  estimate  of  the  delay  spread  A /  is  obtained  by  fitting  the 
Gaussian  distribution  to  a  sampling  of  the  instantaneous  frequency. 

4.4  MID-LATITUDE  HF-NVIS  CHANNEL  MODELS 

Table  7  lists  the  current  narrowband  HF-NVIS  channel  models.  The  model  is  a  narrowband  model  of 
Example  9  Equation  (17)  is  repeated  here  for  convenience: 

h(t,  t)  =  ai(t)S(r-Ti)  +  a2(i)5(r-  r2)  +  a3(f)5(r  -  r3). 

The  “Disturbed”  model  comes  from  the  ITU-R  F.1487  [34],  where  only  the  delay  difference  of  7  ms  is 
reported.  For  consistency,  the  first  path  is  assigned  a  delay  of  t\  —  1.7  ms  [2].  This  delay  then  sets  the 
delay  of  the  second  path  as  i~2  =  8.7  ms. 

Table  7.  Mid-latitude  narrowband  HF-NVIS  channel  models;  hk  magnitude  (abs);  rk  delay  (ms);  A fk 
Doppler  spread  (Hz);  Doppler  shift  is  zero. 


Descriptor 

hi 

Tl 

A/i 

hi 

T2 

A/a 

h3 

T3  a  f3 

Reference 

Good 

0.4 

1.7 

0.2 

[2] 

Medium 

0.4 

1.7 

0.5 

0.2 

4.6 

0.5 

[2] 

Poor 

0.4 

1.7 

1.0 

0.2 

4.6 

1.0 

0.15 

8.3  1.0 

[2] 

Disturbed 

1.0 

1.7 

1.0 

1.0 

8.7 

1.0 

[34] 

This  report  assumes  the  simple  fading  on  the  taps  [76,  Section  1.2.2]: 

ak(t)  =  e+j27TfD’kt{sk  +  gk(t)}, 


where  {#&(£)}  is  a  zero-mean,  complex-valued,  circular,  Gaussian  stochastic  process  with  variance  2 crk. 
The  gain  of  the  tap  is  the  sum  of  the  specular  and  diffuse  “power” 

h2k  :=  E[\ak(t)\2}  =  \sk\2  +  2a2k. 


The  ratio  of  the  specular  to  the  diffused  power  sets  the  A  Tactor: 


K  = 


2al ' 


The  Doppler  spread  A  fk  is  determined  by  the  power  spectra  of  the  fading  processes.  A  Gaussian  power 
spectrum  is  a  standard  assumption  for  mid-latitude  medium-haul  HF  [34,  Eq.  2]: 


Pk{f)  := 


af,k 


y/2i t 


exp  - 


2(j2f,k , 


The  multimode  power  spectrum  is  not  needed  because  the  wide  bandwidth  separates  modes  [34].  This 
particular  spectrum  sets  Doppler  spread  as  [34,  page  3]: 


Afk  : —  2(7 f,k- 

This  report  extracts  these  channel  parameters  from  the  HF-NVIS  channel  soundings  and  compares  with  the 
3-D  ray-tracing  code. 
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4.5  HIGH-LATITUDE  HF  AND  HF-NVIS  MODELS 


There  is  long-term  interest  in  high-latitude  HF  channels  and  an  emerging  interest  in  the  high-latitude 
HF-NVIS  channels.  This  section  reviews  selected  high-latitude  HF  channel  phenomena  to  complete  the 
discussion  of  delay  spread  versus  Doppler  spread  (Figure  52)  supporting  future  channel  modeling  efforts. 
At  the  system  level,  Figure  52  reports  high-latitude  channels  have  relatively  small  delay  spread  but  larger 
Doppler  spreads  compared  to  the  mid-latitude  channels.  The  larger  Doppler  spreads  are  typically  attributed 
to  magnetic  activity  [82]. 


Table  8.  Doppler  spreading  and  magnetic  activity  [82,  Table  2]. 


Magnetic  condition 

Time 

Doppler  spread  (Hz) 

Delay  spread  (ms) 

Quiet 

Day 

0.09 

0.080 

Perturbed 

Day 

6.2 

0.812 

Post-perturbation 

Day 

0.6 

0.545 

Quiet 

Night 

1.6 

0.583 

Disturbed 

Night 

7.6 

1.159 

Disturbed 

Night 

12.6 

0.571 

At  the  physical  level,  Feinblum  and  Horan  offer  the  following  insights  on  the  high-latitude  ionosphere 
[23]: 


The  most  striking  feature  of  the  mid-to-high-latitude  ionosphere  is  probably  the  highly 
localized  region  of  very  low  electron  densities  at  the  F2  peak  which  is  observed  at  night  at 
geographic  latitudes  from  about  35°  N  to  about  70°N  . . .  referred  to  as  the  “main  trough” 

...  the  “mid-latitude  trough,”  the  “high-latitude  trough,”  or  the  “Ottawa  trough.” 

Figure  57  shows  typical  poleward  and  equatorward  walls  that  delineate  the  trough.  Although  the 
trough  is  a  night-time  feature  during  the  winter  months,  the  trough  changes  the  nature  of  the  HF-NVIS 
propagation  from  vertical  reflections  off  the  ionosphere  to  wall  reflections  [73]: 

The  reduction  in  electron  density  means  that  the  signal  radiated  by  HF  radio  communication 
systems  operating  under  the  trough  may  no  longer  be  ionospherically  reflected  along  the  great 
circle  path.  However,  the  signal  can  still  be  reflected  from  gradients  in  the  poleward  and 
equatorward  walls  of  the  trough  or  scattered  from  irregularities  embedded  in  the  trough  or  in 
the  auroral  region  which  lies  just  poleward  of  the  trough. 

Figure  57  also  shows  two  mid-latitude  HF  links  marked  in  blue.  These  links  are  within  the  trough. 
There  is  a  high-latitude  channel  above  63°  North  marked  in  red.  Figure  58  reports  the  delay  and  Doppler 
spreads  measured  over  both  these  medium-haul  HF  links  (4.64-18.38  MHz).  The  mid-latitude  channels 
spread  to  4.5  ms  in  delay  while  the  high-latitude  channels  spread  to  14  Hz  in  Doppler.  Thus,  the 
medium-haul  HF  soundings  show  the  same  Delay-Doppler  structure  as  the  HF-NVIS  channels — the 
mid-latitude  channels  spread  in  delay  while  the  high-latitude  channels  spread  in  Doppler. 

The  final  observation  on  these  high-latitude  channels  occurs  near  the  auroral  oval  (shaded  region)  [73]. 

Within  the  polar  cap  ionosphere,  features  of  particular  note  are  convecting  patches  and 
sun-aligned  arcs  of  enhanced  electron  density  that  can  lead  to  signals  arriving  in  directions 
displaced  from  the  great  circle  path  by  up  to  100°. 
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Figure  57.  Typical  trough,  its  bounding  walls,  and  the  auroral  oval  (shaded)  [73]  (figure  courtesy  of  Radio 
Science). 
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Figure  58.  HF  medium-haul  channels  in  the  mid-  and  upper  latitudes  [73]. 


Consequently,  HF  radio  systems  operating  in  the  high-latitudes  must  adapt  waveforms,  frequencies,  and 
antenna  patterns  to  account  for  the  large  Doppler  spreads  and  non- vertical  paths. 
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4.6  SUMMARY 


This  section  reviewed  two  approaches  for  modeling  the  HF-NVIS  channel  functions: 

•  Realizations  of  channel  parameters  for  HF  channel  emulation 

•  The  delay  Doppler  spreads  describing  classes  of  HF  channels 

The  superb  wideband  HF  measurements  of  Wager  and  Goldstein  taken  in  1995  are  an  excellent  example  of 
the  latter  [82].  Figure  58  shows  that  the  medium-haul  HF  channels  cluster  in  the  mid-latitude  or 
high-latitude  behaviors  of  Figure  52.  If  these  measurements  are  still  extant,  a  re-analysis  to  extract  ray  path 
statistics  would  deliver  channel  parameters  for  the  wideband  HF  medium-haul  in  the  high 
latitudes — without  the  substantial  cost  of  high-latitude  field  measurements.  There  still  remains  the  need  for 
high-latitude  HF-NVIS  measurements  to  clarify  antenna  performance,  measure  ray  path  delays  and 
Doppler  shifts  as  affected  by  the  magnetic  field,  estimate  polarization- switching  along  the  ray  paths,  and 
compare  absorption  models  to  the  measurements. 

Table  7  shows  the  limited  HF-NVIS  channel  models  available  for  HF-NVIS  channel  emulation. 
Although  these  limited  narrowband  HF-NVIS  channel  soundings  have  been  boosted  by  the  wideband 
HF-NVIS  soundings  of  TrellisWare  [79],  publication  of  that  massive  database  is  pending  as  of  this  report. 

The  wideband  soundings  described  in  this  report  and  the  follow-on  reports  serve  three  purposes.  The 
first  and  immediate  goal  is  to  provide  wideband  HF-NVIS  channel  models  for  both  channel  simulation  and 
emulation.  Second  and  more  lasting,  is  that  these  measurements  validate  the  3-D  ray-tracing  code. 
Consequently,  credible  mid-latitude  HF-NVIS  channels  may  be  simulated  without  costly  field  tests.  Third, 
and  contributing  to  the  understanding  of  the  HF  wireless  channel,  the  3-D  ray-tracing  code  explicitly 
models  the  antenna’s  coupling  into  the  electric  field.  Consequently,  the  coherent  interaction  of  multiple 
antennas  and  time- varying  polarizations — necessary  for  the  SIMO,  MIMO,  and  dual-polarization 
configurations — is  correctly  modeled  and  linked  to  the  polarization  evaluation  on  the  O-  and  X-Modes.  As 
such,  this  3-D  ray-tracing  code  is  better  suited  for  HF  wireless  engineering  than  the  power-tracking  codes. 
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5.  THE  HF-NVIS  CHANNEL  SOUNDING  SYSTEM 


Figure  59  is  a  block  diagram  of  the  HF-NVIS  channel  sounding  system.  The  transmitter  is  an  85-ft 
horizontal  dipole  antenna  located  at  the  Model  Range  of  Space  and  Naval  Warfare  Systems  Center  Pacific. 
(SSC  Pacific).  This  stationary  transmitter  broadcasts  a  pseudorandom  noise  (PN)  sequence  for  wideband 
channel  sounding.  The  receiver  is  carried  in  a  mobile  van  equipped  with  a  two-port  antenna.  The  receiver 
carried  in  the  mobile  van  phases  the  ports  so  that  the  sounding  sequence  is  simultaneously  received  on  a 
vertical  whip  that  has  the  standard  HF  antenna  pattern  similar  to  Figure  5  and  a  vertical  loop  that  has  an 
HF-NVIS  pattern  similar  to  Figure  7.  Consequently,  this  system  design  allows  simultaneous  side-by-side 
comparisons  between  a  standard  HF  system  and  an  HF-NVIS  system  in  the  same  over-the-air  experiment. 
The  following  sections  detail  key  components  of  this  system.  Section  5.1  documents  the  HF  antennas  and 
their  measured  patterns.  Section  5.2  describes  the  sounding  waveform  and  the  archiving  scheme  of  the 
channel  soundings.  Section  5.3  reviews  the  channel  estimation  scheme  and  synthetic  examples  showing 
the  effect  of  additive  noise.  Section  5.4  makes  explicit  the  relative  calibration  used  in  this  report. 
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Figure  59.  Block  diagram  of  the  HF-NVIS  channel  sounding  system. 

5.1  TRANSMIT  AND  RECEIVE  HF  ANTENNAS 

The  transmit  antenna  is  a  horizontal  half-wave  wire  dipole  approximately  84.5  ft  long.  It  is  hung  about 
30  ft  in  the  air  between  two  wooden  masts.  The  antenna’s  braided  wire  arms  are  connected  to  a  Unadilla® 
W2DU  HF  1:1  balun.  Also  connected  to  the  balun  is  about  50  ft  of  coaxial  cable  that  goes  through  a  hole 
in  the  ground  to  the  transmitter  in  an  underground  tunnel.  The  balun  has  a  ferrite  sleeve  that  chokes  the 
current  on  the  cable  shield.  The  ground  plane  beneath  the  antenna  and  masts  is  a  metal  mesh  that 
approximates  a  perfect  ground  for  HF  frequencies.  The  dipole  and  ground  plane  occupy  the  open  space 
above  a  circular  outdoor  pattern  range  shown  in  Figure  60. 

Figure  61  shows  the  receiving  half-loop  mounted  on  the  top  of  the  mobile  receiver.  This  half-loop  is 
based  on  the  design  in  [83]  and  adapted  for  the  HF  band  in  [13]  and  [14].  The  loop  is  connected  to  a  metal 
ground  plane  slightly  smaller  than  the  minivan  roof.  This  ground  plane  is  secured  to  a  luggage  rack  and 
not  electrically  connected  to  the  van  body. 
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Figure  60.  Transmitter:  84.5-ft  horizontal  dipole  located  approximately  30  ft  above  the  ground  (yellow  line). 


Figure  61 .  Mobile  FIF-NVIS  receiver  carrying  the  dual-function  antenna. 


The  half-loop  is  electrically  small  for  the  HF  band  and  on  a  small  ground  plane,  and  is  inefficient  at 
receiving  radiation.  However,  because  the  external  HF  environment  is  much  noisier  compared  to  higher 
frequencies  due  to  the  favorable  propagation  of  interference  over  long  distances,  the  receive  system  is 
external  noise  limited.  This  means  that  the  internal  noise  in  the  receive  chain  is  insignificant  compared  to 
the  external  noise  received  by  the  loop.  Making  the  loop  larger  would  increase  the  received  signal  strength 
but  also  the  received  noise,  and  therefore  would  not  significantly  improve  the  SNR.  The  ends  of  the  loops 
float  about  an  inch  over  the  ground  plane.  These  ends  are  connected  to  coaxial  cables  with  4: 1  Guanella 
unun  transformers  [28].  The  coaxial  cables  from  each  port  then  connect  to  the  two  input  ports  of  a  180° 
hybrid  coupler.  The  sum  port  adds  both  ports  of  the  loop  and  generates  a  dipole  radiation  pattern.  The 
difference  port  adds  both  loop  ports  out  of  phase  and  generates  a  loop  radiation  pattern. 

The  hemispherical  radiation  pattern  of  the  half-loop  on  its  ground  plane  was  measured  at  the  outdoor 
pattern  range  (Model  Range).  Figures  62,  63,  and  64  show  the  directivities  in  three  cut  planes  measured  at 
5.8  MHz.  Figure  62  shows  the  approximate  azimuth  cut  plane  around  the  antenna.  This  is  only 
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approximate  because  the  probe  dipole  antenna  was  5°  above  the  horizon,  which  was  as  low  to  the  horizon 
as  the  dipole  could  get  without  touching  the  ground  plane.  The  dominant  polarization  is  Eq. 

In  Figure  62,  the  dipole  (whip)  mode  of  operation  produces  the  expected  omni-directional  Eq  radiation 
pattern  in  azimuth.  The  loop  mode  Eq  pattern  has  two  nulls  at  loop  broadside  as  expected.  In  Figure  63, 
half  of  the  elevation  cut  plane  that  is  broadside  to  the  loop  is  shown.  The  dominant  polarizations  for  both 
antenna  modes  is  Eq.  The  whip  mode  has  a  null  at  zenith  and  a  maximum  at  the  horizon.  The  loop  mode 
has  a  null  at  the  horizon  and  a  maximum  at  zenith.  In  Figure  64,  the  whip  mode  again  has  a  maximum  at 
the  horizon  and  a  null  at  zenith  in  the  dominant  Eq  polarization.  The  loop  mode  has  a  high  Eq  level  both 
at  the  zenith  and  close  to  the  horizon  because  the  cut  plane  is  in  the  plane  of  the  loop.  At  the  horizon,  the 
pattern  still  goes  to  zero  because  the  antenna  is  directly  above  a  ground  plane.  Overall,  the  measured  half¬ 
loop  pattern  has  the  main  characteristics  of  an  ideal  half-loop  above  a  ground  plane. 


Azimuth  Cut  at  9= 5°;  5.8  MHz 
0° 


Figure  62.  Radiation  pattern  of  dual-port  loop  at  5.8  MHz  in  the  azimuthal  plane. 

Elevation  Cut:  Broadside  to  Loop;  5.8  MHz 


Figure  63.  Radiation  pattern  of  dual-port  loop  at  5.8  MHz  in  the  elevation  plane  orthogonal  to  the  plane  of 
the  loop. 

5.2  SOUNDING  WAVEFORM 

Figure  65  shows  the  channel-sounding  waveforms  organized  into  a  frame  bundle.  Each  frame  bundle 
consists  of  58  1-second  frames.  Each  frame  consists  of  22  PN  bursts  followed  by  dead  time.  Recording 
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Elevation  Cut:  Plane  of  the  Loop;  5.8  MHz 
0° 


Figure  64.  Radiation  pattern  of  dual-port  loop  at  5.8  MHz  in  the  elevation  plane  of  the  loop. 

a  frame  bundle  starts  on  each  minute.  The  final  2  seconds  of  each  minute  are  reserved  for  saving  of  the 
recordings  of  each  software-defined  radio  to  disk. 
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Figure  65.  Frame  bundle. 
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Each  PN  burst  spn(^)  is  a  waveform  obtained  from  the  BPSK  sequence  with  a  square-root  raised 
cosine  pulse  shaping: 


NP 

5Pn(^)  —  ^  ^  ^n<S\/Rc(^  _  ^-^sym )  5 
n=l 


(18) 


where  Np  is  the  length  of  the  PN10  sequence:  Np  =  1023,  the  6n’s  are  the  PN10  symbols:  bn  =  ±1,  and 
5\/rc(£)  denotes  the  square-root  raised  cosine  pulse  with  symbol  period  Tsym.  For  the  channel  sounding 
shown  in  Figure  65,  the  symbol  rate  is  Rs ym  =  25  kbps  so  that  symbol  duration  is 


Tsym  —  RSym 


40  [fj,s\ . 


Each  frame  sp(t )  consists  of  Np  =  22  PN  bursts 


nb 

SF(t)  =  ^  <§PN(f  -  n^Tp), 

nb= 1 


where  each  burst  has  duration 

Tp  =  NpTsym  =  40.092  [ms]. 

The  Np  —  22  PN  bursts  produce  a  transmission  “on”  time  of  duration 

Tp  =  NBTT  =  900.240  [ms]. 

The  remainder  of  the  1 -second  frame  consists  is  the  dead  time  lasting 

Tp  =  99.760  [ms]. 

Each  frame  bundle  is  recorded  in  a  .WAV  file  corresponding  to  the  whip  or  NVIS  antenna  and  stamped 
with  the  date  and  UTC  stamp.  On  10  February  2016,  370  files  were  recorded  from  the  whip  and  another 
370  files  were  recorded  from  the  loop  as  the  van  (Figure  61)  traveled  from  Point  Foma,  San  Diego,  CA,  to 
El  Centro,  CA,  and  back  (Figures  73  and  74)  starting  at  8:30  a.m.  PST  and  ending  at  3:08  p.m.  PST.  Table 
9  shows  the  labeling  and  interpretation  for  two  sample  files. 

Table  9.  Labeling  scheme  for  the  .WAV  files. 

File  Antenna  Date  PST 

NVis_20l602l0_l62800.wav  Loop  10-Feb-2016  8:28  a.m. 

GND_20160210_230800.wav  Whip  10-Feb-2016  3:08  p.m. 


Conversion  from  UTC  time  to  local  PST,  on  this  date,  requires  that  8  hours  be  subtracted  from  the 
UTC  hours: 


162800  ^  (16  -  8):28:00  =  08:28  PST 

230800  ^  (23  -  8):08:00  =  15:08  PST 
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5.3  CHANNEL  ESTIMATION 


The  HF-NVIS  channel  is  sounded  by  the  PN  bursts  in  the  frame  bundle  of  Figure  65.  The  HF-NVIS 
channel  is  estimated  by  standard  correlation  that  assumes  the  channel  is  stationary  over  all  the  PN  bursts  in 
a  frame  as  a  sum  K  delays  where  K  is  unknown 


K 

h(t )  =  -  Tk). 

k= 1 

The  received  signal  «§#(£)  is  modeled  as  noise-corrupted  multipath, 

/oo 

h(t')sPN(t  -  t')dt'  +  g(t), 

-OO 

where  {g(t)}  is  the  zero-mean  noise.  Correlation  processing  simply  filters  the  received  signal  against  the 
“extraction”  signal  sx(t )  designed  to  register  a  strong  unit  peak  with  each  PN  sequence.  A  noisy  version 
of  a  stationary  channel  is  then  extracted  as  follows: 

sx*SR,(t)  =  h*sx  *sT(t)  +  sx  *g(t) 

K 

=  ^hk(sx  *  sT)(t  -  Tk)  +  sx  *s(£) 
k= 1 
K 

~  ^2  hkS(t  -  rk)  +  sx  *  g{t). 

k= 1 

The  extraction  signal  is  scaled  to  register  a  unit  response  for  the  channel.  Consequently,  the  noise  level  in 
the  channel  estimate  is  suppressed  by  the  magnitude  of  this  filter.  The  next  example  computes  the 
suppression  for  white  noise. 


Example  10  (White  Noise)  Although  zero-mean  continuous-time  white  noise  exists  only  as  a  generalized 
random  process  with  covariance  [57,  Eq.  9-40] 

Rgg(T)  =  °2gS(T) 

this  limiting  case  provides  insight  into  the  filtered  noise: 

/oo  poo  _ 

dtl  /  dt2  Sx{tl)SX{t2)Rgg{t2  -  h) 

-OO  J  — OO 

/oo  poo  _ 

dti  /  dt2  sx{ti)sx{t2)5(t2  -  ti) 

-OO  J  — OO 

/OO 

\sx{t)\2dt. 

-OO 

Channel  estimation  is  performed  in  the  digital  domain.  The  extraction  signal  is  the  time-reversed  copy 
of  the  digital  PN  signal  extended  by  a  few  extra  symbols.  The  time-limited  square-root  raised  cosine 
function  as  the  square-root  raised  cosine  function  observed  over  Nsym  symbols  centered  at  t  —  0  is  defined 
as: 

5VRC (NSym',t)  ~  s^/wc(t)  x  1[-Td/2,Td](£);  =  NsymTsyin. 


E[\sx*g(t)\2}  = 
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The  time-limited  PN  signal  as  the  time-limited  version  of  Equation  (18)is  defined  as: 


NP 

5pN(^Vsym?  f)  =  ^  ^  ^ri^v/RC  (^symj  t  ^Xsym )  •  (19) 

n=  1 


The  digital,  time-limited  square-root  raised  cosine  function  by  sample  using  Nsps  samples  per  symbol  is 
defined  as: 

k 

£PN  (^Vsymj  ^Ysps?  ^k)  ~  ^PN  (^sym5  ^k)  5  ^ k  ~  ~T~  i  f s  ~  -Rsym-^sps* 

Is 


The  extraction  signal  is 


sx(tk)  •= 


5pN(^Vsym?  -Ngpsj  ^ k ) 


1 2  • 
1 2 


||5pN(^Vsym?  -^sps) 

Table  10  lists  the  parameters  for  the  extraction  signal  specific  for  this  channel  sounding  and  sets 

||spN(^VSym,^VSps)||2  =  31  dB. 


Table  10.  Parameters  for  the  digital  PN  signal. 


Variable 

Value 

Description 

fs 

100 

Digital  sampling  rate  (kHz) 

T  p 

0.35 

roll-off  factor 

N sym 

11 

Number  of  symbols 

Nsps 

4 

Number  of  samples  per  symbol 

-f^sym 

25 

Symbol  rate  (kHz) 

Example  11  (Simulated  Channel)  Figure  66  reports  on  the  performance  of  this  channel  estimation 
scheme  using  the  PN  signal  parameterized  by  Table  10.  The  upper  panel  shows  a  simulated  channel  h(t ) 
with  taps  spread  over  10  ms.  The  received  signal  sn(t)  is  the  time-limited  PN  signal  spn(^)  of  Equation 
(19)  filtered  through  this  channel  and  corrupted  by  white  Gaussian  noise: 

Sr (t)  =  h  *  spnW  +  gif)- 


The  specified  SNR  sets  the  variance  of  the  noise  reported  as  the  “wideband  input  noise”: 

SNR  = 


II h  *  spnIH 


af, 


The  lower  panel  reports  on  the  channel  estimate: 

^est(^)  =  sx  *  SR(t)  +  sx  *  g(t)  =  h(t)  +  sx  *  #(*)• 

The  peaks  of  the  channel  register  while  the  output  noise  level  is  reduced  by  approximately  31  dB. 
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Channel 


Correlation  Recovery:  pT= 0  dBW 


milli-sec:  SNR=10  dB;  Wideband  Input  Noise=  1  dB;  Output  Noise=  -29  dB 

Figure  66.  Simulated  channel  (top  panel);  Estimated  channel  (lower  panel). 


Example  1 1  uses  a  bandwidth  of  25  kHz  that  is  oversampled  by  a  factor  of  four:  fs  =100  kHz.  The 
digital  noise  in  this  example  has  flat  spectrum  spread  over  the  entire  100-kHz  band.  The  next  example 
assesses  channel  recovery  when  the  in-band  noise  is  confined  to  the  25-kHz  band. 


Example  12  (In-Band  Noise)  Figure  67  shows  the  low -pass  filter  that  determines  the  in-band  noise 
spectrum.  Figure  68  shows  the  performance  of  this  channel  estimation  scheme  against  a  realization  of  this 
in-band  noise.  The  upper  panel  shows  a  simulated  channel  h(t )  with  taps  spread  over  10  ms.  The  received 
signal  is  the  spn(^)  of  Equation  (19)  filtered  through  this  channel  and  corrupted  by  the  filtered  Gaussian 
noise: 

SR{t )  =  h*  spn{t)  +  g{t). 

The  specified  SNR  sets  the  variance  of  the  noise  reported  as  the  u  In-Band  Input  Noise”: 


SNR  = 


||  h  *  SpN  | 


(7ft 


2 

2 


The  lower  panel  of  Figure  68  reports  on  the  channel  estimate: 


^ESt(^)  =  *  Sr(T)  +  Sx  *  g(t)  =  h(t)  +  Sx  *  g{p)- 


The  peaks  of  the  channel  still  register,  the  output  noise  is  reduced  by  approximately  25  dB,  and  the  noise 
floor  is  smoother. 
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Magnitude  Response  (dB) 


Figure  67.  Low-pass  filter  to  generate  in-band  noise. 


Channel 


Correlation  Recovery:  />  = 0  dBW 


Figure  68.  Simulated  channel  (top  panel);  Estimated  channel  (lower  panel). 
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5.4  CALIBRATION 


The  preceding  section  reported  on  channel  estimation  in  an  ideal  digital  simulation  using  the  implicit 
assumptions: 

•  Perfect  waveform  synthesis 

•  Perfect  synchronization  between  transmitter  and  receiver 

•  No  delay  between  transmitter  and  receiver 

•  Perfect  amplitude  control 

•  Zero  system  noise 

In  contrast,  all  real-world  channel  sounding  systems  do  not  generate  ideal  waveforms,  suffer  from  jitter 
in  synchronization,  and  exhibit  end-to-end  delays  caused  by  processing.  The  end-to-end  amplitude  can 
only  be  known  by  measurement.  Therefore,  timing  and  gain  calibrations  are  necessary  to  compensate  these 
effects. 

On  2  March  2016,  the  transmitter  and  receiver  shown  in  Figure  59  were  connected  with  cables  for 
back-to-back  testing.  This  section  reports  on  the  testing  and  determines  the  necessary  timing  and 
amplitude  compensation — assuming  linear  operation  of  the  amplifiers  and  demodulators — so  that  the 
frame-bundle  processing  reports  no  delay  and  unit  amplitude  for  this  back-to-back  channel.  Therefore, 
applying  this  calibration  to  the  over-the-air  measurements  should  only  register  the  HF  channels,  the 
transmitter’s  power,  and  the  antenna  effects. 

Figure  69  shows  the  received  signal  reshaped  into  the  format  of  the  ideal  frame-bundle  of  Figure  65. 
The  crisp  edges  and  straight  lines  demonstrates  that  the  channel  sounding  system  is  maintaining  tight 
synchronization 

Figure  70  records  the  first  PN  sequence  of  the  first  frame  of  the  frame  bundle.  This  waveform  shows 
that  the  transmitter  is  producing  high-quality  analog  waveforms  and  digitizing  the  waveform  in  the 
receiver. 

Figure  7 1  is  a  close-up  view  comparing  the  real  part  of  this  measured  waveform  to  its  ideal  digital 
waveform  delayed  to  compensate  for  the  system  delays.  The  scaling  was  selected  to  show  that  the  received 
wave  aligns  with  the  ideal  waveform.  That  is,  any  processing  done  in  the  receiver  can  assume  the  received 
waveform  is  well-modeled  by  its  associated  ideal  waveform. 

Figure  72  correlates  the  real-world  received  waveform  sr(£)  against  the  idealized  extraction  signal 
sx(t)  using  the  compensating  system  delays.  The  single  peak  shows  that  the  system  delays  have  been 
compensated.  The  peak  value  of  0.025136  is  amplitude  scaling  necessary  to  produce  the  unit  height.  A 
slightly  improved  value  is  the  “gCai  —  —15.9976”  (dB)  reported  on  the  forthcoming  channel  estimation 
plots. 
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Figure  69.  Received  signal  shaped  into  the  frame-bundle  format. 
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Figure  70.  First  PN  sequence  of  the  first  frame  of  the  frame  bundle. 
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Figure  71 .  Comparing  the  real  part  of  the  measured  received  waveform  to  an  ideal  digital  waveform.  The 
scaling  was  selected  to  show  the  timing  alignment. 
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Figure  72.  Extraction  of  the  back-to-back  channel  to  compute  amplitude  scaling. 
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6.  SIGNAL,  NOISE,  AND  INTERFERENCE  ALONG  THE  MOBILE  ROUTE 


Figure  73  shows  the  route  of  the  mobile  HF-NVIS  receiver  on  10  February  2016  marked  in  red.  The 
starting  point  is  Point  Loma,  CA,  at  8:30  a.m.  PST.  The  outbound  route  traveled  on  Interstate  8  East 
arriving  near  El  Centro,  CA,  at  10:30  a.m.  PST.  The  return  loop  traveled  north  to  the  Salton  Sea,  crossed 
Anza-Borrego  Desert  Park,  and  returned  through  the  Cuyamaca  Mountains.  Figure  74  isolates  the  route  for 
clarity.  Figure  75  reports  the  ground  range  as  a  function  of  the  time  along  the  route. 


Figure  73.  Route  of  the  mobile  HF-NVIS  vehicle. 


Time-to-Distance:  10-Fcb-2016 


Figure  74.  Mobile  HF-NVIS  route. 


Figure  75.  Ground  range  as  a  function  of 
time. 


6.1  SIGNAL  POWER  ALONG  THE  ROUTE 

The  transmit  power  delivered  to  the  horizontal  dipole  was  approximately  90  Watts  at  5.82  MHz  over  a 
25-kHz  band.  The  transmit  signal  consists  of  a  repeating  sequence  of  frame  bundles  as  shown  in  Figure  65. 
Each  frame  bundle  starts  at  the  beginning  of  each  minute  and  lasts  58  seconds.  These  frame  bundles  are 
simultaneously  received  through  the  whip  and  loop  antennas  enabling  side-by-side  comparisons  between 
the  standard  HF  and  the  HF-NVIS  channels. 
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Figures  76  and  77  compare  the  “raw”  power  received  by  the  whip  and  loop  antennas.  The  quotes 
around  “raw”  mean  that  the  power  is  computed  by  averaging  the  entire  received  digital  signal  recorded  in 
each  frame  bundle.  All  58  frames  are  included  in  this  average,  including  the  dead  time  at  the  end  of  each 
frame.  The  dB  reported  is  computed  directly  from  the  digital  audio  signal  in  the  .WAV  files — no 
conversion  to  dBW  is  applied.  Comparing  the  plots  shows  that  the  noise  floor  of  the  loop  antenna  is 
approximately  25  dB  lower  than  the  whip,  excepting  the  interference  “spikes”  and  the  broad  region 
between  10:00  and  1 1:00  PST.  This  comparison  opens  the  conjecture  that  the  up-ward  pointing  loop 
rejects  much  of  the  HF  noise  arriving  from  the  side  that  is  swamping  the  whip. 


Figure  76.  Raw  whip  power.  Figure  77.  Raw  loop  power. 

6.2  NOISE  POWER  ALONG  THE  ROUTE 

Figure  78  and  79  support  the  conjecture  that  the  HF-NVIS  loop  antenna  is  rejecting  the  local  HF  noise 
by  comparing  the  “noise  only”  portion  at  the  end  of  each  frame.  The  dB  reported  is  computed  directly 
from  the  digital  audio  signal  in  the  .WAV  files — no  conversion  to  dBW  is  applied.  Comparison  of  these 
plots  show  that  the  loop  noise  is  still  25  dB  lower  than  the  whip  noise,  excepting  the  interference  “spikes” 
and  the  broad  region  between  10:00  and  1 1:00  PST.  Comparison  with  the  preceding  raw  power  plots  also 
shows  that  the  signal  received  by  the  whip  is  essentially  lost  once  the  ground  wave  is  blocked — there  is 
negligible  shift  in  power  levels  between  the  “on”  time  in  the  frame  bundle  in  comparison  to  the  dead  time. 


Figure  78.  Raw  whip  noise.  Figure  79.  Raw  loop  noise. 

A  complicating  factor  between  the  difference  in  unwanted  received  power  between  whip  and  loop 
antennas  is  the  presence  of  interfering  signals.  The  next  section  examines  possible  sources  of  these 
interfering  signals:  cell  towers  are  associated  with  the  “spikes”  and  the  city  of  El  Centro,  CA  lines  up  with 
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the  broad  swath  between  10:00  and  1 1:00  PST.  One  solution  for  mitigating  cell  tower  effects  is  to  move 
the  receiver  a  few  hundred  feet.  However,  the  swath  of  urban  interference  is  persistent  as  the  vehicle 
travels  through  El  Centro,  CA.  This  RF  interference  opens  a  design  objective  for  placing  an  HF-NVIS 
antenna  on  a  platform:  these  platform-mounted  antennas  should  adapt  to  their  platform  by  maximizing  a 
null  at  the  horizon  over  all  azimuth  angles. 

6.3  INTERFERENCE  ALONG  THE  ROUTE 

Figure  80  reports  on  possible  sources  of  the  spikes  observed  in  the  HF-NVIS  channel.  The  assumption 
is  that  the  loop  antenna  has  an  upward  pointing  pattern.  Therefore,  overhead  sources  such  as  high-voltage 
power  lines  passing  over  the  freeway  and  cell  towers  immediately  adjacent  to  the  freeway  couple  into  this 
antenna. 


9  iO  II  12  13  14  15 

PST:  1 0-Feb-20 1 6 


Figure  80.  Average  power  received  by  the  whip  antenna. 


Identifying  the  source  of  all  spikes  is  problematic.  For  example,  Figures  81  and  82  locate  a  cell  tower 
near  Desert  View  point  and  high-voltage  lines  parallel  to  the  highway.  However,  two  of  these  cell  towers 
are  in  relatively  isolated  areas.  Figures  83  and  84  locate  an  isolated  radio  tower  at  10:03  PST.  Figures  85 
and  86  locate  isolated  radio  tower  at  12:14  PST.  The  isolation  of  these  towers  makes  associating  the  spikes 
with  these  sources  more  plausible  than  the  cell  tower  near  high-voltage  lines. 
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Track:  10-Feb-2016 


Figure  82.  Radio  tower  near  Desert  View  Point  at  09:54  PST. 
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Track:  10-Feb-2016 


Exit  Street  View 


Figure  84.  Radio  tower  after  passing  Ocotillo  at  10:03  PST  but  prior  to  Coyote  Wells. 
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Track:  10-Feb-2016 


Figure  86.  Cell  tower  enroute  at  12:14  PST  approaching  Ocotillo  Wells. 
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7.  HF-NVIS  CHANNEL  ESTIMATES  ALONG  THE  MOBILE  ROUTE 


Figure  87  displays  all  the  HF-NVIS  channel  soundings  by  plotting  the  magnitude  of  the  delay-spread 
function  h(t ,  r)  described  in  Section  4.1.  Each  vertical  slice  at  time  t  is  an  estimate  of  the  impulse 
response  of  the  channel  with  delay  r  measured  in  milliseconds.  The  vertical  striations  are  the  interference 
sources  reported  in  Section  6.3.  The  strong  responses  at  the  beginning  and  end  of  the  day  are  caused  by  the 
mobile  vehicle  departing  and  arriving  back  at  the  transmitter’s  site.  The  dB  scale  was  selected  to  display 
the  HF-NVIS  channels  for  most  of  the  day. 


NVIS  Channel  Estimate: /c=5.82 
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Figure  87.  HF-NVIS  channel  estimate  extracted  from  the  loop  antenna  over  the  entire  experiment.  The 
ground  wave  and  the  3-Hop  response  are  not  labeled. 


On  this  day  and  at  this  frequency,  this  mobile  wideband  HF-NVIS  channel  is  well-modeled  by 
Equation  (16)  and  repeated  here  for  convenience: 


h(t,  r )  =  ho(t,  t )  +  hx(t,  t), 
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where  the  channel  functions  for  the  O-  and  X-Modes  admit  the  approximations 

h0(t,  t )  «  aF)0(t)5(T  -  rFi0)  +  aFtOX(t)S(T  -  rF> ox) 

hx(t,r )  «  aF,x(0^(r  ~  rF,x)  +aF,xo(i)<5(>  “  tf,xo)- 

The  ground  wave,  the  2.5-Hop,  and  the  3-Hop  channels  are  not  included  in  this  model.  The  2.5-Hop  is 
tentatively  identified  similar  to  the  “FEsF”  ray  paths  of  Figure  37  produced  by  the  sporadic  E-layer  (See 
[17,  Figure  6.27(f)]).  From  the  perspective  of  HF  radio  testing,  this  2-Hop  model  captures  the  dominant 
1-Hop  ray  paths  and  models  self-interference  with  the  2-Hop  ray  paths.  However,  the  1-Hop 
approximation  is  still  a  credible  model  for  HF  radio  testing  because  most  of  the  2-Hop  ray  paths  arrive 
with  10  dB  or  more  attenuation  relative  to  the  1-Hops: 

h(t,  t )  =  aFt0(t)5(T  -  rF,o)  +  aF:X(t)S(T  -  rFjx)-  (20) 

Figure  87  shows  that  both  the  1-Hop  ray  paths  dominate  across  most  of  the  day.  The  O-Mode  arrives  first 
with  typically  greater  power.  The  arrival  of  the  modes  is  slightly  bowed  in  the  middle  of  the  plot 
corresponding  to  the  distance  of  the  mobile  from  the  transmitters.  The  arrival  time  of  both  modes  also 
modulated  starting  an  1 1:00  PST  and  lasting  to  13:00  PST.  The  current  conjecture  is  that  this  modulation  is 
a  TID  in  the  f2  layer. 

This  section  examines  the  channel  at  selected  times  along  the  route.  Each  channel  assessed  against  the 
ionogram  and  the  synthetic  channel  model  predicted  by  the  3-D  ray-tracing  code.  The  fading  processes, 
the  delay  spread,  the  Doppler  shifts,  and  the  Doppler  spreads  are  discussed  relative  to  each  channel.  This 
section  concludes  by  gathering  the  preceding  observations  to  model  the  1-Hop  model  of  Equation  (20)  in  a 
form  suitable  to  program  a  channel  emulator. 
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7.1  HF-NVIS  CHANNEL  AT  09:07  PST 


Figure  88  locates  the  mobile  receiver  approximately  40  km  north-east  of  the  transmitter  as  it  enters 
the  foothills  of  the  Cuyamaca  Mountains  near  Flynn  Springs,  CA.  This  is  shown  on  the  map  in  Figure  89. 
There  appears  to  be  little  interference  from  cell  towers  or  other  sites  along  the  route.  The  vehicle  is  pointed 
at  83.5°  and  traveling  60-80  mph. 


Track:  10-Feb-2016 


Figure  89.  Receiver  location  marked  by  yellow  “thumbtack”  at  09:07  PST. 


67 


Figure  90  shows  the  magnitude  of  estimated  channel  averaged  over  the  sounding.  The  channel  contains 
multiple  ray  paths:  the  ground  wave,  the  1-Hops,  2-Hop,  2.5-Hop,  3-Hop,  and  4-Hop.  Each  hop,  excepting 
the  ground  wave,  contains  the  O-  and  X-Mode  ray  paths. 


HF-NVIS  Channel  Estimate:  10-Feb-2016  9:07:00  PST; /c=5.82 


Figure  90.  Channel  estimate  at  09:07  PST. 


The  distance  between  the  transmitter  and  receiver  is  approximately 

do  ~  40.1223  [km] 

and  sets  the  delay  in  the  ground  wave  as 

tG  =  —  m  0.1337  [ms], 
c0 

The  accuracy  of  this  ground  wave  delay  is  one  verification  that  the  system  delays  are  correctly  calibrated  in 
Section  5.4.  The  discussion  now  focuses  on  these  measured  hop  times  and  the  hop  times  computed  by  the 
3-D  ray-tracing  code. 
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Figure  91  shows  the  nearest  ionogram  measured  at  Point  Arguello,  CA,  at  09: 15  PST.  The  ionogram 
shows  four  vertical  bounces  where  the  O-Mode  arrives  before  the  X-Mode.  The  black  curve  on  this 
ionogram  is  the  estimated  plasma  frequency  that  determines  the  electron  density;  see  Equation  (7).  The 
ray-tracing  code  uses  the  hybrid  ionogram  obtained  by  spatially  averaging  the  ionograms  at  Point 
Arguello,  CA,  Boulder,  CO,  and  Austin,  TX.  Figure  92  shows  this  hybrid  plasma  frequency  profile  used  in 
the  3-D  ray- tracing  code. 
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Figure  91.  Ionogram  at  09:15  PST.  Open  source  ionograms  courtesy  of  the  Global  Ionospheric  Radio 
Observatory  at  the  University  of  Massachusetts,  Lowell  [64], 


Plasma  Frequency  (MHz) 


Figure  92.  Estimated  plasma  frequency  profile  at  09:07  PST. 

Figures  93  and  94  show  the  1-Hop  O-  and  X-Modes  traveling  northeast  from  the  transmitter  to  the 
mobile  receiver.  The  Earth’s  magnetic  field  bends  the  O-Mode  north  and  the  X-Mode  south.  Figures  95 
and  96  show  the  polarization  along  these  paths.  The  O-Mode  path  starts  with  LCP  and  ends  with  RCP 
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passing  through  vertical  polarization.  The  X-Mode  path  starts  with  RCP  and  ends  with  LCP  passing 
through  horizontal  polarization.  This  polarization  switching  and  the  ground  reflection  determine  that  the 
2-Hop  ray  paths  having  the  OO-  and  XX-Mode  structure  discussed  in  Example  8. 


Figure  93.  1-Hops  at  09:07  PST.  Figure  94.  1-Hops  at  09:07  PST  (top  view). 

LCP  LCP 


Figure  95.  Polarization  along  the  1-Hop  Figure  96.  Polarization  along  the  1-Hop 

O-Mode  ray  path  at  09:07  PST.  X-Mode  ray  path  at  09:07  PST. 

Table  1 1  compares  the  hop  times  between  this  sounding,  the  ray  paths  computed  by  the  3-D  ray-tracing 
code,  and  estimates  obtained  from  the  ionogram.  The  table  reports  that  the  hop  times  produced  by  the 
ray-tracing  code  have  a  consistent  positive  bias  of  0.1  ms.  The  bottom  layer  of  the  plasma  frequency 
profile  is  critical  for  computing  accurate  ray  paths — raising  the  layer  causes  hop  times  to  increase  because 
of  the  longer  path  in  the  ionosphere.  Comparison  between  the  measured  and  simulated  plasma  profile 
shows  the  simulated  profile  starts  at  a  slightly  lower  altitude.  Although  tweaking  the  ionogram  layers  can 
force  the  simulated  hop  times  to  line  up  with  the  measured  hop  times,  this  “fitting  the  model”  after  the 
experiment  is  not  an  approach  that  supports  HF  mission  planning.  Rather,  improved  ionogram  modeling 
boosts  performance  of  any  ray-tracing  codes. 
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Table  1 1 .  Comparing  hop  times  between  the  sounding,  simulation,  and  the  nearest  ionogram  (Point 
Arguello,  CA)  at  09:07  PST. 


09:07  PST 

Measured 

Simulation 

Ionogram 

1-Hop  O-Mode 

1 .66  ms 

1.71  ms 

1 .67  ms 

1-Hop  X-Mode 

1.74  ms 

1.86  ms 

1 .85  ms 

2-Hop  O-Mode 

3.34  ms 

3.42  ms 

3.34  ms 

2-Hop  X-Mode 

3.48  ms 

3.70  ms 

3.36  ms 

2.5-Hop  O-Mode  (?) 

4.38  ms 

— 

— 

2.5-Hop  X-Mode  (?) 

4.46  ms 

— 

3-Hop  O-Mode 

5.00  ms 

5.12ms 

5.07  ms 

3-Hop  X-Mode 

5.22  ms 

5.54  ms 

— 

4-Hop  O-Mode 

6.06  ms 

6.84  ms 

6.80  ms 

This  section  closes  with  an  analysis  of  the  Doppler  shifts,  the  Doppler  spreading,  and  the  fading. 
Figure  97  shows  an  estimate  of  the  delay-spread  function  /i(£,  r)  of  Equation  (12)  where  “f  ’  is 
determined  by  the  timing  of  the  PN  sequences.  Averaging  the  magnitude  of  this  delay-spread  function 
produces  the  estimated  HF-NVIS  channel  estimate  of  Figure  90. 


NVIS  Channel  Estimate:  10-Fcb-2016  9:07:00  PST;/C=5.82 


200  400  600  800  1000 

PN  index 

Figure  97.  Delay  spread  estimate  at  09:07  PST. 

The  closely  spaced  O-  and  X-Mode  channels  show  slow  modulations  with  the  noise  floor  switching  to 
lower  level  midway  through  the  sounding.  The  ground  wave  is  just  visible  at  0.14  ms.  The  1-Hop  ray  paths 
are  the  strongest  just  before  2  ms.  The  2-Hop  at  3.4  ms  show  the  strongest  modulation.  The  3-Hop  are  at  5  ms 
while  the  unidentified  2.5-Hop  are  between  the  2-  and  3-Hop. 
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Figure  98  shows  an  estimate  of  the  scattering  function  described  in  Equation  (13).  This  estimate  is 
computed  by  taking  the  discrete  Fourier  transform  of  the  delay-spread  function  h(t,  r)  over  the  time 
variable  t  where  uniform  time  samples  are  obtained  by  resampling  over  the  PN  bursts  to  handle  the  “dead” 
time  shown  in  Figure  65.  The  expectation  operator  in  Equation  (13)  is  approximated  by  Wiener  filtering. 
This  estimate  of  the  scattering  function  shows  the  delay  and  Doppler  of  the  1-,  2-,  2.5-,  3-,  and  4-Hop. 

This  estimate  even  registers  the  ground  wave  and  shows  the  ground  wave  with  Doppler  shift 
fo, grnd  =  —0.5833  Hz  producing  a  credible  vehicle  speed  of  [87,  Eq.  4.9] 

v  =  ^c/d,G  RND  =  -67  [mph], 

where  the  minus  sign  has  the  vehicle  moving  away  from  the  transmitter. 

HF-NVIS  Scattering  Function:  10-Feb-2016  9:07:00  PST;/C=5.82 
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Figure  98.  Scattering  function  estimated  at  09:07  PST. 

The  limitation  of  the  scattering  function  estimate  is  that  the  relatively  small  Doppler  shifts  and  Doppler 
spreads  are  not  well-resolved.  Therefore,  examination  of  the  individual  ray  paths  is  necessary  to  extract 
more  reliable  Doppler  and  fading  estimates. 

Figure  99  displays  the  amplitude  and  phase  of  the  1-Hop  O-Mode  channel  obtained  from  the  delay- 
spread  function  at  the  O-Mode  delay: 

ho(t)  =  \ho(t)\Zho(t)  =  h(t,ro). 

The  amplitude  in  the  upper  panel  shows  slow  and  fast  fading.  The  Rayleigh  fading  o  ~  0.159  Hz 
measured  around  the  median  line  is  a  one-parameter  compromise  to  this  non- stationary  fading.  However, 
the  lower  panel  shows  the  unwrapped  phase  and  an  excellent  fit  of  the  quadratic  phase  model 

Zh0(t)  «  <t>o{t)  :=  0o  +  0i£  +  02 12.  (21) 


The  linear  term  estimates  the  Doppler  shift  reported  on  this  plot: 


0  10  20  30  40  50 

sec:  fNQ=  0.159  Hz  (Rayleigh  fading) 


Figure  99.  1-Hop  O-Mode  amplitude  and  phase  at  09:07  PST. 


Figure  100  reports  the  same  computations  for  the  1-Hop  X-Mode  channel’s  amplitude  and  phase: 

hx(t)  =  \hx(t)\Zhx(t)  =  h(t,TX). 

The  amplitude  in  the  upper  panel  shows  fast  and  slow  fading  with  a  relatively  stationary  region  between  10 
and  30  seconds.  The  Rayleigh  fading  /^o  ~  0.33  Hz  is  a  forced  “fit”  to  this  envelope  as  measured  by  the 
median  line.  The  lower  panel  shows  the  phase  and  the  excellent  fit  of  a  quadratic  phase  model: 

Ahx(t)  «  0x(i)  :=  0O  +  01 1  +  02i2  (22) 


The  Doppler  shift  reported  on  this  panel  is  determined  from  the  linear  term: 


fx  = 


1  d0xm 
2tt  dt  1  j 


t= 0 


Table  12  compares  the  Doppler  shifts  extracted  from  the  O-  and  X-Mode  over-the-air  measurements 
with  estimates  obtained  from  the  vehicle  speed  extracted  from  the  scattering  function  and  simulated 
elevation  angles  computed  by  the  3-D  ray-tracing  code.  At  best,  the  comparison  shows  that  all  the  Doppler 
shifts  are  small. 


Table  12.  Comparing  measured  and  estimated  Doppler  shifts  in  Hz  based  on  the  estimated  mobile  speed 
v  =  -67  mph. 
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10-Feb-2016  9:07:00  PST:  X-Modc 


Figure  100.  1-Hop  X-Mode  amplitude  and  phase  at  09:07  PST. 


Doppler  structure  is  revealed  by  the  time-varying  transfer  function  repeated  here  for  convenience: 


/oo 

e-j2n fTh(t,  r)d,T. 

-OO 


Figure  101  shows  the  magnitude  of  the  time-varying  transfer  function  for  this  channel.  The  channel’s 
frequency  response  as  a  function  of  time  is  a  “spectrogram.” 


-50  0  50 

kHz:  /  =5.82  MHz 

J  c 

Figure  101.  Spectrogram  at  09:07  PST. 
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Nulls  caused  by  different  Doppler  shifts  on  the  1-Hop  O-  and  X-Modes  appear  to  wrap  over  the 
25-kHz  band.  Appendix  C  shows  how  the  1-Hop  linear  phase  model 

h(t,  t)  =  hoe+i27rfot5(t  —  tq)  +  hxe+^2n^xtS{t  —  rx) 


produces  the  time-varying  transfer  function 

T(i,  /)  =  h0e+j2*fote-j2*fT°  +  hxe+j27rfxte~j27rf ^ 

with  constant  phase  along  straight  lines  in  the  time-frequency  plane.  However,  the  wrapping  in  Figure 
101  is  better  modeled  by  the  quadratic  phase  model  of  Equations  (21)  and  (22).  The  corresponding  time- 
varying  transfer  function 

T(t,  f)  =  hoe+jMt)e~j2nfTO  +  hxe+j^e-j2 nfr* 

has  constant  phase  along  quadratic  curves  in  the  time-frequency  plane.  Figure  102  overlays  the  lines  of 
constant  quadratic  phase  on  the  spectrogram. 


10-Feb-2016  9:07:00  PST 
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Figure  102.  Spectrogram  at  09:07  PST  and  the  1-Hop  model  with  quadratic  phase. 

This  relatively  good  fit  of  the  quadratic  phase  model  argues  that  the  1-Hop  delay-spread  function 
should  have  the  form 


h(t,  r)  =  hoe+^°^5(t  —  to)  +  hxe+^x^5(t  —  rx). 

Most  channel  emulators  handle  the  linear  phase  term  as  a  Doppler  shift.  The  significance  of  this  quadratic 
phase  model  for  channel  emulation  is  an  open  question. 

Turning  to  the  problem  of  identifying  the  fading  processes,  the  best  approach  is  based  on  the  ray  path 
physics.  Rician  fading  is  claimed  to  be  a  plausible  fading  model  on  the  ray  paths  [17,  Section  7.7.6]. 
However,  fading  identification  is  confounded  by  fast  and  slow  fading  [87,  Section  4.10.2].  Figure  100 
shows  the  X-Mode  fading  is  relatively  stationary  over  the  first  30  seconds.  Figure  99  shows  the  O-Mode 
fading  has  both  fast  and  slow  fading.  Fading  identification  is  restricted  to  the  X-Mode  over  the  first 
30  seconds. 
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Figure  103  shows  the  “goodness-of-fit”  ranking  estimated  over  this  window.  Although  Rician  fading 
has  a  significance  level  over  80%  for  this  particular  estimate,  these  “goodness-of-fit”  tests  are  extremely 
sensitive  to  the  correlation  between  samples.  This  particular  plot  was  produced  using  every  10th  sample. 
Figure  104  fits  selected  distributions  to  the  sample  histogram  of  the  X-Mode  fading  restricted  to  the  first  30 
seconds.  Although  no  distribution  shows  a  compelling  fit,  the  physical  basis  for  the  Rician  fading  makes  it 
a  natural  candidate  for  the  fading  process. 
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Figure  1 03.  “Goodness-of-fit”  for  the  1  -Hop  X-Mode  Figure  1 04.  Fitting  probability  density  functions  to  the 
amplitude  at  09:07  PST.  1  -Hop  X-Mode  amplitude  at  09:07  PST. 


Table  13  encapsulates  this  channel  in  a  form  suitable  for  a  channel  emulator  as  discussed  in 
Section  4.4.  The  model  employs  only  the  1-Hop  O-  and  X-Modes.  The  total  path  loss  is  scaled  to  the 
1-Hop  O-Mode  and  encompasses  both  the  transmit  and  receive  antenna,  including  the  latter’s  orientation 
and  velocity  on  the  mobile. 

Table  13.  Wideband  HF-NVIS  channel  model;  5.82  MHz;  10  February  2016  09:07  PST. 


09:07  PST 

Path 

Total  path 

Doppler 

Fade 

Rician  K 

delay 

loss 

shift 

rate 

factor 

1-Hop  O-Mode 

1 .66  ms 

OdB 

-0.12  Hz 

0.16Hz 

8.7  dB 

1-Hop  X-Mode 

1.83  ms 

-5.9  dB 

-0.04  Hz 

0.33  Hz 

10.0  dB 
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7.2  HF-NVIS  CHANNEL  AT  09:41  PST 


Figure  105  locates  the  mobile  receiver  approximately  90  km  due  east  of  the  transmitter.  Figure  106 
shows  the  mobile  receiver  traveling  through  the  rolling  foothills  near  Boulevard,  CA.  There  appears  to 
be  little  interference  from  cell  towers  or  other  sites  along  the  route.  The  vehicle  is  pointed  at  116.3°  and 
traveling  60-80  mph. 


Track:  10-Feb-2016 


Figure  106.  Receiver  location  marked  with  yellow  “thumbtack”  at  09:41  PST. 
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Figure  107  shows  the  estimated  channel  contains  multiple  ray  paths.  Table  14  compares  the  hop  times 
between  this  sounding,  the  ray  paths  computed  by  the  3-D  ray-tracing  code,  and  estimates  obtained  from 
the  nearest  ionogram. 


HF-NVIS  Channel  Estimate:  10-Feb-2016  9:41:00  PST;/C=5.82 


Figure  107.  Channel  estimate  at  09:41  PST. 


Table  14.  Comparing  hop  times  between  the  sounding,  simulation,  and  ionogram  at  09:41  PST. 


09:41  PST 

Measured 

Simulation 

Ionogram 

1-Hop  O-Mode 

1.71  ms 

1 .69  ms 

1 .60  ms 

1-Hop  X-Mode 

1.83  ms 

1 .85  ms 

1 .73  ms 

2-Hop  O-Mode 

3.37  ms 

3.34  ms 

3.23  ms 

2-Hop  X-Mode 

3.64  ms 

3.65  ms 

3.33  ms 

2.5-Hop  O-Mode  (?) 

4.43  ms 

— 

— 

2.5-Hop  X-Mode  (?) 

4.55  ms 

— 

— 

3-Hop  O-Mode 

5.08  ms 

5.00  ms 

— 

3-Hop  X-Mode 

— 

— 

— 

4-Hop  O-Mode 

6.71  ms 

6.66  ms 

— 

Figure  108  is  the  nearest  ionogram.  This  ionogram  shows  the  1-Hop  and  2-Hop  paths  arrive  with  the 
O-Mode  followed  by  the  X-Mode. 
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Figure  108.  Point  Arguello,  CA  ionogram  at  09:45  PST.  Open  source  ionograms  courtesy  of  the  Global 
Ionospheric  Radio  Observatory  at  the  University  of  Massachusetts,  Lowell  [64]. 


The  ray-tracing  code  uses  the  hybrid  ionogram  obtained  by  spatially  averaging  the  ionograms  at  Point 
Arguello,  CA,  Boulder,  CO,  and  Austin,  TX.  Figure  109  shows  the  plasma  frequency  fp{h )  obtained  from 
this  hybrid  ionogram.  Figure  110  and  111  shows  the  simulated  1-Hop  ray  paths  with  the  O-Mode  bending 
north  and  the  X-Mode  bending  south.  Figures  112  and  113  show  that  the  polarization  switches  on  both  ray 
paths. 
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This  polarization  switching  and  the  reflection  determine  that  the  2-Hop  ray  paths  will  be  OO-  and 
XX-Modes  discussed  in  Example  8.  Table  14  reports  that  the  ray-tracing  code  produces  very  accurate  hop 
times  using  these  1-,  2-,  3-,  and  4-Hop  models. 

Figures  110  and  111  display  the  O-  and  X-Mode  1-Hops  traveling  from  the  transmitter  to  the  mobile 
receiver  due  east.  This  orientation  reveals  the  maximum  bending  caused  by  the  Earth’s  magnetic  field. 
Figures  112  and  113  show  the  polarization  along  these  paths.  The  O-Mode  path  starts  with  LCP  and  ends 
with  RCP  passing  through  vertical  polarization.  The  X-Mode  path  starts  with  RCP  and  ends  with  LCP 
passing  through  horizontal  polarization.  This  polarization  switching  and  the  ground  reflection  determine 
that  the  2-Hop  ray  paths  reported  in  Table  14  are  the  OO-  and  XX-Modes  discussed  in  Example  8. 


Figure  110.1  -Hops  at  09:41  PST.  Figure  111.1  -Hops  at  09:41  PST  (top  view). 

LCP  LCP 


Figure  1 1 2.  Polarization  along  the  1  -Hop  Figure  1 1 3.  Polarization  along  the  1  -Hop 
O-Mode  ray  path  at  09:41  PST.  X-Mode  ray  path  at  09:41  PST. 
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Figure  1 14  estimates  the  delay-spread  function  h(t ,  r)  where  is  determined  by  the  timing  of  the  PN 
sequences.  The  closely  spaced  O-  and  X-Mode  channels  show  slow  modulations  against  a  relatively  flat  noise 
floor.  The  1-Hop  ray  paths  are  the  strongest  just  before  2  ms.  The  2-Hop  at  3.4  ms  show  the  strongest 
modulation.  The  3-Hop  are  at  5  ms  while  the  unidentified  2.5-Hop  are  between  the  2-  and  3-Hops. 

NVIS  Channel  Estimate:  10-Feb-2016  9:41:00  PST;/C=5.82 
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Figure  114.  Delay  spread  estimate  at  09:41  PST. 

Figure  115  estimates  the  scattering  function  for  this  channel  This  estimate  is  computed  by  taking  the 
discrete  Fourier  transform  of  the  delay-spread  function  h(t ,  r)  over  the  time  variable  t  without  resampling 
over  the  PN  bursts  to  handle  the  “dead”  time  shown  in  Figure  65.  The  expectation  operator  in  Equation  (13) 
is  approximated  by  Wiener  filtering.  This  estimate  of  the  scattering  function  shows  the  delay  and  Doppler  of 
the  1-,  2-,  2.5-,  3-,  and  4-Hop.  This  estimate  shows  the  increasing  delays  between  the  O-  and  X-Modes  as 
the  hops  increase  and  the  slight  Doppler  differences  between  the  hops. 


Hz:  Wiener  Filter — no  gap  correction 


Figure  115.  Scattering  function  estimated  at  09:41  PST. 
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The  limitation  of  the  scattering  function  estimate  is  that  the  relatively  small  Doppler  shifts  and  Doppler 
spreads  are  not  well-resolved.  Therefore,  examination  of  the  individual  ray  paths  is  necessary  to  extract 
more  reliable  Doppler  and  fading  estimates. 

Figure  116  displays  the  amplitude  and  phase  of  the  1-Hop  O-Mode  channel  obtained  from  the  delay- 
spread  function  at  the  O-Mode  delay 

ho  (t)  =  \ho(t)\Zho(t)  =  h(t,ro). 

The  amplitude  in  the  upper  panel  appears  to  show  multiple  fading  processes.  However,  the  mean  level  is 
-44  dB  while  the  variation  is  only  ±3  dB.  That  is,  the  fading  is  approximately  Rician.  Consequently,  the 
Rayleigh  fading  estimate  /jv,o  ~  0.477  Hz  measured  around  the  median  line  is  a  reasonable  fit  of  a  one- 
parameter  model  to  this  nearly  Rician  fading.  The  lower  panel  shows  the  phase  and  a  quadratic  fit 

ZhQ(t)  «  (po (t)  :=  <po  +  <f>i t  +  02 t2 


where  the  linear  term  estimates  the  Doppler  shift  as 


f  1  d(^°  (+\ 

,0  -df(t) 


«  -0.247  [Hz]. 

t= o 


0  10  20  30  40  50 

sec:  /  =  0.477  Hz  (Rayleigh  fading) 


Figure  116.  1-Hop  O-Mode  amplitude  and  phase  at  09:41  PST. 
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Figure  117  reports  the  same  computations  for  the  1-Hop  X-Mode  channel.  The  amplitude  in  the  upper 
panel  shows  relatively  stationary  fast  fading  so  that  the  Rayleigh  fading  fx,x  ~  0.416  Hz  is  a  reasonable 
estimate  from  this  measurement. 


10-Feb-2016  9:41:00  PST:  X-Modc 


Figure  1 1 7.  1  -Hop  X-Mode  amplitude  and  phase  at  09:41  PST. 
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Figure  118  shows  the  spectrogram  for  this  channel  sounding  and  the  quadratic  phase  lines.  This 
spectrogram  is  the  magnitude  of  the  time-varying  transfer  function.  Similar  to  the  spectrogram  at 
09:07  PST,  the  estimated  Doppler  extracted  from  the  O-  and  X-Modes  provides  an  excellent  fit  to  the 
nulls  wrapping  over  the  25-kHz  band.  Consequently,  these  Doppler  estimates,  albeit  small,  are  reported  in 
the  fading  description.  Table  15  encapsulates  this  channel  in  a  form  suitable  for  a  channel  emulator.  The 
model  employs  only  the  1-Hop  O-  and  X-Modes.  The  total  path  loss  is  scaled  to  the  1-Hop  O-Mode  and 
encompasses  both  the  transmit  and  receive  antennas,  including  the  latter’s  orientation  and  velocity  on  the 
mobile  receiver. 


10-Feb-2016  9:41:00  PST 
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Figure  118.  Spectrogram  at  09:41  PST. 

Table  15  encapsulates  this  channel  in  a  form  suitable  for  channel  emulator.  The  model  employs  only 
the  1-Hop  O-  and  X-Modes.  The  total  path  loss  is  scaled  to  the  1-Hop  O-Mode  and  encompasses  both  the 
transmit  and  receive  antenna,  including  the  latter’s  orientation  and  velocity  on  the  mobile. 

Table  15.  Wideband  HF-NVIS  channel  model;  5.82  MHz;  10  February  2016  09:41  PST. 


09:41  PST 

Path 

Total  path 

Doppler 

Fade 

Rician  K 

delay 

loss 

shift 

rate 

factor 

1-Hop  O-Mode 

1.71  ms 

OdB 

-0.02  Hz 

0.48  Hz 

12.6  dB 

1-Hop  X-Mode 

1.83  ms 

-6.3  dB 

-0.08  Hz 

0.41  Hz 

13.4  dB 
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7.3  HF-NVIS  CHANNEL  AT  10:23  PST 


Figure  119  locates  the  mobile  at  the 
southeastern  part  of  the  loop  almost  160  km  from 
the  transmitter.  The  terrain  is  flat  desert.  There 
is  no  interference  from  cell  towers  or  El  Centro, 

CA.  The  vehicle  is  pointed  at  89.4°  or  almost 
due  east  and  traveling  between  60  to  80  mph. 

Therefore,  the  maximum  Doppler  shift  is  0.7  Hz. 

Figure  120  shows  the  estimated  channel  consists  of 
two  strong  ray  paths  followed  by  some  smaller 
ray  paths.  Figure  121  is  the  nearest  ionogram 
showing  the  ray  paths  at  5.82  MHz  arrive  in 
O-Mode  followed  by  X-Mode  order.  Table  16 
compares  the  hop  times  between  this  sounding,  the 
ray  paths  computed  by  the  3-D  ray-tracing  code, 
and  estimates  obtained  from  the  ionogram.  The 
ray-tracing  code  predicts  O-Mode  ray  paths  also 
arrive  before  their  associated  X-Mode  ray  paths.  RiJ“re  1 1 9-  Receiver  location  on  the  track  at  1 0:23 
The  ray-tracing  code  uses  a  hybrid  ionogram 
obtained  by  spatially  averaging  the  ionograms  at 

Point  Arguello,  CA,  Boulder,  CO,  and  Austin,  TX.  The  ray-tracing  code  produces  accurate  hop  times  for 
the  O-Mode  hops  but  overestimates  the  X-Mode  hops  time  by  appropriately  0.2  ms.  The  table  also  lists 
two  paths  observed  between  the  2-Hop  and  3-Hop.  These  are  tentatively  identified  “FEsF”  ray  paths  of 
Figure  37  that  might  be  produced  by  a  sporadic  E-layer  [17,  Figure  6.27(f)]. 

HF-NVIS  Channel  Estimate:  10-Feb-2016  10:23:00  PST;/ =5.82 
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Figure  120.  Channel  estimate  at  10:23  PST. 


- 1 

1 - 

-45.9 

1.71ms 

1 - j 

i - 1 

-49.6 

1.83ms 

-66.7 

-67.3 

3.29ms  4  44ms 
.  -67.7 

ft  |i  4.55ms 

0  2  4  6  8  10 

absolute  delay  (ms):  A  =11;  g  =-15.9976 

J  v  '  sym  °cal 


Track:  10-Feb-2016 
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Figure  121.  Point  Arguello,  CA  ionogram  at  10:30  PST.  Open  source  ionograms  courtesy  of  the  Global 
Ionospheric  Radio  Observatory  at  the  University  of  Massachusetts,  Lowell  [64]. 


Table  16.  Comparing  O-  and  X-Mode  hop  times  between  the  sounding,  simulation,  and  the  nearest 
ionogram  (Point  Arguello,  CA)  at  10:23  PST. 


10:23  PST 

Measured 

Simulation 

Ionogram 

1-Hop  O-Mode 

1.71  ms 

1.73  ms 

1.67  ms 

1-Hop  X-Mode 

1.83  ms 

1.91  ms 

1.80  ms 

2-Hop  O-Mode 

3.29  ms 

3.33  ms 

3.27  ms 

2-Hop  X-Mode 

— 

3.67  ms 

3.53  ms 

2.5-Hop  O-Mode  (?) 

4.44  ms 

— 

4.30  ms 

2.5-Hop  X-Mode  (?) 

4.55  ms 

— 

4.60  ms 

3-Hop  O-Mode 

— 

4.96  ms 

4.93  ms 

3-Hop  X-Mode 

— 

5.47  ms 

— 
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Figure  122  shows  the  plasma  frequency  fp{h )  estimated  by  spatially  averaging  the  nearest  three 
ionograms.  This  estimated  plasma  frequency  determines  the  electron  density  profile  Ne{h)  and  the 
associated  index  of  refraction  used  by  the  3-D  ray-tracing  code  to  compute  the  ray  paths.  Comparison  of 
this  hybrid  plasma  frequency  with  the  plasma  frequency  at  Figure  121  shows  that  the  gross  shapes  are 
similar.  However,  the  profile’s  lower  edge  is  critical  for  any  ray-tracing  coding.  Comparison  shows  that  the 
hybrid  plasma  frequency  may  start  too  low  and  introduce  a  bias  for  longer  times  in  the  ionosphere. 


Plasma  Frequency  Profile 


Plasma  Frequency  (MHz) 

Figure  122.  Plasma  frequency  at  10:23  PST  estimated  from  the  ionograms  at  Point  Arguello,  CA,  Boulder, 
CO,  and  Austin,  TX. 

Figures  123  and  124  show  the  O-  and  X-Mode  1-Hop.  The  rays  travel  west  to  east  so  are 
approximately  transverse  to  the  Earth’s  magnetic  field.  The  magnetic  field  bends  the  O-Mode  to  the  north 
and  the  X-Mode  to  the  south.  Figures  125  and  126  report  the  polarization  along  these  paths.  The  O-Mode 
path  starts  with  LCP  and  ends  with  RCP  while  passing  through  vertical  polarization.  The  X-Mode  path 
starts  with  RCP  and  ends  with  LCP  while  passing  through  horizontal  polarization.  This  polarization 
switching  and  the  ground  reflection  determine  that  the  2-Hop  ray  paths  reported  in  Table  16  are  the  OO- 
and  XX-Modes  discussed  in  Example  8. 


Ray  Altitude  Ray  Altitude 


Figure  1 23.  1  -Hops  at  1 0:23  PST.  Figure  1 24.  1  -Hops  at  1 0:23  PST  (top  view). 
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LCP  LCP 


Figure  125.  Polarization  along  the  1-Hop  Figure  126.  Polarization  along  the  1-Hop 
O-Mode  ray  path  at  1 0:23  PST.  X-Mode  ray  path  at  1 0:23  PST. 


Rather  than  present  the  scattering  function,  this  section  extracts  the  fading  and  the  Doppler  shifts  from 
the  1-Hop  O-  and  X-Modes  in  the  delay-spread  function.  Figure  127  shows  an  estimate  of  the  delay-spread 
function.  The  average  of  this  function  across  the  PN  sequence  is  the  channel  estimate  of  Figure  120.  The 
closely  spaced  O-  and  X-Mode  channels  show  slow  modulations.  The  noise  floor  rises  toward  the  end  of 
the  sounding. 


NVIS  Channel  Estimate:  10-Feb-2016  10:23:00  PST; /c=5.82 
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Figure  127.  Delay-spread  estimate  at  10:23  PST. 


Figure  128  reports  the  amplitude  \ho  \  and  phase  ZHq  of  the  1-Hop  O-Mode  channel.  The  upper  panel 
shows  that  the  amplitude  is  relatively  stationary  over  the  first  40  seconds.  Rayleigh  fading  is  estimated  to 
be  /at, o  ~  0.5  Hz  measured  around  the  median  line.  The  lower  panel  shows  the  phase  and  a  quadratic  fit 
where  the  linear  term  estimates  the  Doppler  shift.  The  deep  fades  at  9,  25,  50  seconds  cause  the  jumps  in 
the  phase-unwrapping.  The  quadratic  fit  provides  a  more  robust  estimate  of  the  Doppler  shift  in  the 
presence  of  these  phase  discontinuities  rather  than  averaging  the  instantaneous  frequency.  Figure  129 
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reports  the  same  computations  for  the  1-Hop  X-Mode  channel.  The  upper  panel  shows  multiple  fading 
regions  with  the  noisy  burst  near  40  seconds  causing  the  perturbation  in  the  phase. 
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Figure  129.  1-Hop  X-Mode  amplitude  and  phase  at  10:23  PST. 

The  best  approach  to  fading  identification  is  based  on  the  physics  of  the  problem.  Although  Rician 
fading  has  been  presented  plausible  fading  model  on  the  ray  paths  [17,  Section  7.7.6],  fading  identification 
is  compounded  by  the  presence  of  multiple  fading  processes.  Figure  128  shows  the  O-Mode  fading  is 
relatively  stationary  compared  to  the  X-Mode  fading  in  Figure  129.  Therefore,  fading  identification  is 
restricted  to  the  O-Mode.  The  increasing  background  noise  limits  the  fading  analysis  window  to  the  first 
48  seconds.  Figure  130  shows  one  “goodness-of-fit”  ranking  estimated  over  this  window.  These  standard 
“goodness-of-fit”  tests  are  very  sensitive  to  outliers,  correlation,  and  the  number  of  samples.  No 
distribution  delivers  a  significance  level  over  60%.  Figure  131  compares  the  “fit”  of  selected  distributions 
to  the  sample  histogram  of  the  O-Mode  fading.  All  distributions  provide  a  reasonable  fit.  The  physical 
basis  for  the  Rician  fading  makes  it  a  natural  candidate  for  the  fading  process. 
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Figure  1 30.  “Goodness-of-fit”  for  the  1  -Hop  O-Mode  Figure  1 31 .  Fitting  probability  density  functions  to  the 
amplitude  at  1 0:23  PST.  1  -Hop  O-Mode  amplitude  at  1 0:23  PST. 


Table  17  encapsulates  this  channel  in  a  form  suitable  for  channel  emulator.  The  model  employs  only 
the  1-Hop  O-  and  X-Modes.  The  total  path  loss  is  scaled  to  the  1-Hop  O-Mode  and  encompasses  both  the 
transmit  and  receive  antenna,  including  the  latter’s  orientation  and  velocity  on  the  mobile. 

Table  17.  Wideband  HF-NVIS  channel  model;  5.82  MHz;  10  February  2016  10:23  PST. 


10:23  PST 

Path 

Total  path 

Doppler 

Fade 

Rician  K 

delay 

loss 

shift 

rate 

factor 

1-Hop  O-Mode 

1.71  ms 

OdB 

0.15Hz 

0.54  Hz 

5.2  dB 

1-Hop  X-Mode 

1.83  ms 

-3.7  dB 

-0.03  Hz 

0.33  Hz 

3.6  dB 

Remarks:  The  3-D  ray  tracing  estimates  the  O-  and  X-Mode  hops  arrive  at  the  vehicle  with  elevation 
angles  approximately  73°  and  72°,  respectively.  The  maximum  Doppler  shift  is  estimated  to  be  0.7  Hz. 
These  angles  put  an  upper  bound  on  the  Doppler  shift  induced  by  vehicular  motion  of  approximately  0.2 
Hz  (see  Appendix  B.l  for  the  calculation).  Although  the  Doppler  shifts  extracted  from  the  fading 
processes  are  within  this  bound,  the  difference  in  signs  between  O-  and  X-Mode  shifts  is  not  understood. 
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7.4  HF-NVIS  CHANNEL  AT  13:00  PST 


Figure  132  tracks  the  mobile  approximately  100  km  northeast  of  the  transmitter  as  it  exits  the  Banner 
Grade  in  the  Cuyamaca  Mountains  north  of  Julian,  CA.  The  vehicle  is  pointed  at  291°  and  traveling  40  to 
60  mph.  Figure  133  is  the  topographic  view  showing  that  the  deep  canyon  containing  the  Banner  Grade 
and  the  terrain  of  the  Cuyamaca  Mountains. 


Track:  10-Fcb-2016 


Figure  133.  Receiver  location  at  13:00  PST  marked  by  the  yellow  “thumbtack.” 

Figure  134  shows  the  estimated  HF-NVIS  channel  is  essentially  the  1-Hop  O-  and  X-Modes.  Table  18 
compares  the  hop  times  between  this  sounding,  the  ray  paths  computed  by  the  3-D  ray-tracing  code,  and 
estimates  obtained  from  the  ionogram.  The  ordering  of  the  ray  paths  is  consistent  with  the  O-Modes 
always  arriving  before  the  X-Modes.  The  3-D  ray  tracing  has  a  consistent  bias  of  +0.1  ms  on  the  1-Hop 
O-Mode  that  produces  a  +0.15  ms  bias  in  the  2-Hop  O-Mode  and  a  +0.3  ms  bias  on  the  3-Hop  O-Mode. 
This  bias  is  larger  for  the  X-Mode  with  a  +0.32  ms  bias  for  the  1-Hop  X-Mode  that  runs  into  a  +0.65  ms 
bias  in  the  2-Hop  X-Mode. 
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HF-NVIS  Channel  Estimate:  02-0ct-2016  13:00:00  PST; /(  =5.82 


Figure  134.  Channel  estimate  from  the  loop  antenna  at  13:00  PST. 


Table  18.  Comparing  O-  and  X-Mode  hop  times  between  the  sounding,  simulation,  and  the  nearest 
ionogram  (Point  Arguello,  CA)  at  13:00  PST. 


13:00  PST 

Measured 

Simulation 

Ionogram 

1-Hop  O-Mode 

1.74  ms 

1.83  ms 

1.73  ms 

1-Hop  X-Mode 

1.84  ms 

2.16  ms 

1.87  ms 

2-Hop  O-Mode 

3.47  ms 

3.63  ms 

3.47  ms 

2-Hop  X-Mode 

3.63  ms 

4.28  ms 

3.67  ms 

2.5-Hop  O-Mode  (?) 

4.33  ms 

— 

— 

2.5-Hop  X-Mode  (?) 

4.43  ms 

— 

— 

3  Hop  O-Mode 

5.13ms 

5.43  ms 

— 
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The  ray-tracing  code  uses  the  hybrid  ionogram  obtained  by  spatially  averaging  the  ionograms  at  Point 
Arguello,  CA,  Boulder,  CO,  and  Austin,  TX.  Figure  135  shows  the  plasma  frequency  obtained  from  this 
hybrid  ionogram.  Figure  136  is  the  nearest  ionogram  measured  at  Point  Arguello,  CA.  Comparison  of  the 
lower  edge  of  the  plasma  frequencies  shows  the  hybrid  is  slightly  higher  than  the  measured  ionogram.  The 
resulting  bias  shows  the  limitations  of  any  ray-tracing  code — the  ray  paths  are  only  as  accurate  as  the 
estimated  plasma  frequency  or,  equivalently,  the  electron  density. 


Plasma  Frequency  (MHz) 


Figure  135.  Hybrid  plasma  frequency  at  13:00  PST. 
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Figure  136.  Ionogram  at  Point  Arguello,  CA  at  13:00  PST.  Open  source  ionograms  courtesy  of  the  Global 
Ionospheric  Radio  Observatory  at  the  University  of  Massachusetts,  Lowell  [64], 
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Figure  137  shows  an  estimate  of  the  delay-spread  function  h(t,  r)  where  is  determined  by  the 
timing  of  the  PN  sequences.  The  average  of  this  function  across  the  PN  sequence  is  the  channel  estimate  of 
Figure  134.  The  closely  spaced  O-  and  X-Mode  channels  show  slow  modulations  over  a  relatively  flat 
noise  floor.  The  1-Hop  ray  paths  are  the  strongest  just  before  2  ms.  The  2-Hop  and  3-Hop  are  negligible. 


NVIS  Channel  Estimate:  10-Feb-2016  13:00:00  PST; /c=5.82 
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Figure  137.  Delay  spread  estimate  at  13:00  PST. 


This  section  extracts  the  fading  and  the  Doppler  shifts  from  the  1-Hop  O-  and  X-Modes  in  the 
delay- spread  function  rather  than  the  resolution-limited  scattering  function.  The  1-Hop  O-Mode  fading  has 
amplitude  and  phase 

h(t,ro )  =:  h0(t)  =  \ho(t)\Zho(t). 

Likewise,  the  1-Hop  X-Mode  fading  has  amplitude  and  phase 

h(t,TX)  -■  hx(t)  =  \hx(t)\Zhx(t). 


Figure  138  reports  1-Hop  O-Mode  fading.  The  upper  panel  shows  that  the  amplitude  is  relatively 
stationary:  the  envelope  varies  ±6  dB  around  a  mean  level  of  approximately  —49  dB.  The  Rician  Kq 
factor  is  the  ratio  of  the  specular  to  the  diffuse  power  [76,  page  52].  Under  the  Rician  assumptions,  the 
O-Mode  has  a  strong  specular  component.  The  lower  panel  shows  the  phase  and  a  quadratic  fit 

Zh0(t)  «  0o  +  0i t  +  02t2 

where  the  linear  term  estimates  the  Doppler  shift. 

The  Rician  Kq  factor  is  the  ratio  of  the  specular  to  the  diffuse  power  [76,  page  52].  Under  the  Rician 
assumptions,  the  O-Mode  has  a  strong  specular  component.  The  lower  panel  shows  the  phase  and  a 
quadratic  fit 

Zh0(t)  «  0o  +  0i t  +  02f 2 1 

where  the  linear  term  estimates  the  Doppler  shift. 
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Figure  139  presents  the  corresponding  plots  for  the  X-Mode.  The  deep  fades  are  reflected  in  the  small 
Rician  Kx  factor.  Under  the  Rician  assumptions,  the  X-Mode  is  tending  toward  Rayleigh  fading.  The  deep 
fades  in  the  envelope  cause  multiple  jumps  in  the  unwrapped  phase.  Therefore,  comparison  with  other 
Doppler  estimators  concludes  this  section. 

The  Doppler  spread  is  estimated  from  the  O-  and  X-Mode  phases  using  Woodward’s  approach  [57]. 
Figures  140  and  141  show  the  estimated  Doppler  spectra  of  the  O-  and  X-Modes.  The  Doppler  spectra  are 
proportional  to  the  PDF  of  the  instantaneous  frequencies  obtained  from  O-  and  X-Mode  phases: 

O-Mode/*  (t )  =  T <^-  (t) ;  X-Mode/*  (t)  =  T  ^  (t) . 

Both  figures  compare  the  histograms  of  these  instantaneous  frequencies  against  the  Gaussian  and  the 
Stable  PDF’s.  The  large  tails  of  instantaneous  frequency  bias  the  mean  and  standard  deviation  of  the 
Gaussian  fit.  Nevertheless,  the  standard  deviation  of  the  Gaussian  determines  the  reported  Doppler  spread 
A/ =  2  x  a  as  discussed  in  Section  4. 1 . 


10-Fcb-2016  13:00:00  PST:  O-Mode 


Figure  138.  1-Hop  O-Mode  amplitude  and  phase  at  13:00  PST. 
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10-Fcb-2016  13:00:00  PST:  X-Mode 


Figure  139.  1-Hop  X-Mode  amplitude  and  phase  at  13:00  PST. 


13:00  PST:  0-Modc/.(I)  13:00  PST:  X-Modc/|.(t) 


Figure  140.  O-Mode  Doppler  spectrum 
estimate  at  13:00  PST. 


Figure  141.  X-Mode  Doppler  spectrum 
estimate  at  13:00  PST. 
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Table  19  reports  these  channel  estimates  in  a  form  suitable  for  channel  emulator.  The  model  employs 
only  the  1-Hop  O-  and  X-Modes.  The  total  path  loss  is  scaled  to  the  1-Hop  O-Mode  and  encompasses  both 
the  transmit  and  receive  antennas,  including  the  latter’s  orientation  and  velocity  on  the  mobile.  The 
positive  Doppler  shifts  are  consistent  with  the  receiver  traveling  toward  the  transmitter. 


Table  19.  Wideband  HF-NVIS  channel  model;  5.82  MHz;  10  February  2016  13:00  PST. 


13:00  PST 

Path 

Total  path 

Doppler 

Doppler 

Fade 

Rician  K 

delay 

loss 

shift 

spread 

rate 

factor 

1-Hop  O-Mode 

1.74  ms 

OdB 

+0.12  Hz 

0.23  Hz 

0.15Hz 

7.2  dB 

1-Hop  X-Mode 

1.84  ms 

-9.6  dB 

+0.14  Hz 

0.67  Hz 

0.32  Hz 

-27.4  dB 
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7.5  HF-NVIS  CHANNEL  AT  13:41  PST 


Figure  142  locates  the  mobile  receiver  approximately  71  km  northeast  of  the  transmitter  and  on  the 
plateau  heading  down  to  Santa  Ysabel,  CA.  The  vehicle  is  pointed  at  268.3°  and  traveling  between  40  to 
60  mph.  Figure  143  is  the  topographic  view  showing  the  escarpment  and  the  terrain  of  the  Cuyamaca 
Mountains. 


Track:  10-Feb-2016 


Figure  143.  Receiver  location  marked  with  yellow  “thumbtack”  at  13:41  PST. 
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Figure  144  shows  the  estimated  HF-NVIS  channel.  The  first  spike  is  the  ground  wave,  the  large  spikes 
are  the  1-Hop  O-  and  X-Modes,  followed  by  a  weak  2-Hop  O-Mode,  and  the  2.5  Hop.  As  a 
communications  channel,  this  HF-NVIS  channel  is  well-modeled  by  the  1-Hop  O-  and  X-Modes. 

HF-NVIS  Channel  Estimate:  10-Feb-2016  13:41:00  PST;/ =5.82 
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Figure  144.  Channel  estimate  from  the  loop  antenna  at  13:41  PST. 

There  is  good  agreement  between  this  measured  and  ideal  accuracy  of  this  ground  wave  delay.  The 
distance  between  the  transmitter  and  receiver  is  approximately 

71.1241  [km]. 

The  ground  wave’s  ideal  delay 

tG=d^-K,  0.2371  [ms] 
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aligns  with  the  measured  delay.  Thus,  the  system  delays  seem  to  be  correctly  calibrated  in  Section  5.4. 
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Table  20  compares  the  hop  times  between  this  sounding,  the  ray  paths  computed  by  the  3-D  ray-tracing 
code,  and  estimates  obtained  from  the  ionogram.  The  ordering  of  the  ray  paths  is  consistent  with  the 
O-Modes  always  arriving  before  the  X-Modes.  The  3-D  ray  tracing  has  a  bias  of  +0.06  ms  on  the  1-Hop 
O-Mode.  This  bias  produces  a  +0.1 1  ms  bias  in  the  2-Hop  O-Mode.  This  bias  is  larger  for  the  X-Mode. 

The  1-Hop  X-Mode  bias  is  +  0.13  ms  that  runs  into  a  +0.34  ms  bias  in  the  2-Hop  X-Mode. 

Table  20.  Comparing  O-  and  X-Mode  hop  times  between  the  sounding,  simulation,  and  the  nearest 
ionogram  (Point  Arguello,  CA)  at  13:41  PST. 


13:41  PST 

Measured 

Simulation 

Ionogram 

1-Hop  O-Mode 

1 .67  ms 

1.73  ms 

1.63  ms 

1-Hop  X-Mode 

1.79  ms 

1 .92  ms 

1.83  ms 

2-Hop  O-Mode 

3.32  ms 

3.43  ms 

3.67  ms 

2-Hop  X-Mode 

3.46  ms 

3.80  ms 

3.53  ms 

2.5-Hop  O-Mode  (?) 

4.39  ms 

— 

— 

2.5-Hop  X-Mode  (?) 

4.51  ms 

— 

— 

Comparison  with  the  bias  at  13:00  PST  shows  that  these  errors  are  decreasing.  Indeed,  only  a 
41 -minute  difference  and  a  30-km  difference  reduces  the  O-Mode  bias  from  0.1  ms  at  13:00  PST  to  +0.06 
ms  at  13:41  PST.  As  expected,  the  hybrid  ionogram  used  by  the  3-D  ray-tracing  code  has  a  slightly  higher 
plasma  frequency  profile  than  the  ionogram  measured  at  Point  Arguello,  CA.  Figure  145  shows  the 
ionogram  at  13:45  PST.  The  black  line  is  the  plasma  frequency  profile.  The  3-D  ray-tracing  code  uses  the 
hybrid  ionogram  obtained  by  spatially  averaging  the  ionograms  at  Point  Arguello,  CA,  Boulder,  CO,  and 
Austin,  TX.  Figure  146  shows  the  plasma  frequency  profile  obtained  from  this  hybrid  ionogram. 
Comparison  of  the  lower  edge  of  the  plasma  frequencies  shows  the  hybrid  plasma  frequency  profile  is 
slightly  higher  than  the  measured  ionogram.  The  resulting  bias  shows  the  limitations  of  any  ray-tracing 
code — the  ray  paths  are  only  as  accurate  as  the  estimated  plasma  frequency  or,  equivalently,  the  electron 
density. 
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Figure  145.  Point  Arguello,  CA,  ionogram  at  13:45  PST.  Open  source  ionograms  courtesy  of  the 
Global  Ionospheric  Radio  Observatory  at  the  University  of  Massachusetts,  Lowell  [64], 


Plasma  Frequency  Profile 
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Figure  146.  Hybrid  plasma  frequency  estimate  at  13:41  PST. 


101 


Figure  147  shows  the  delay-spread  function  h(t,  r)  where  is  determined  by  the  timing  of  the  PN 
sequences.  The  average  of  this  function  across  the  PN  sequence  is  the  channel  estimate  of  Figure  144.  The 
closely  spaced  O-  and  X-Mode  channels  show  slow  modulations  over  a  relatively  flat  noise  floor.  The  1- 
Hop  ray  paths  are  the  strongest  just  before  2  ms.  The  ground  wave,  the  2-,  and  the  2.5-Hop  are  negligible. 
Figures  148  and  149  show  the  respective  O-  and  X-Mode  1-Hop  fading  at  13:41  PST. 


NVIS  Channel  Estimate:  10-Feb-2016  13:41:00  PST;/C=5.82 


200  400  600  800  1000 

PN  index 

Figure  147.  Delay  spread  estimate  at  13:41  PST. 


As  in  the  preceding  sections,  Doppler  and  fading  estimates  of  these  modes  are  obtained  by  extracting 
the  amplitude  and  phase  from  the  delay-spread  function.  Figure  148  displays  the  amplitude  and  phase  of 
the  1-Hop  O-Mode  channel  obtained  by  scanning  the  delay-spread  function  at  the  O-Mode  delay: 


ho(t)  =  \ho(t)\Zho(t)  =  h(t,ro). 

The  amplitude  in  the  upper  panel  shows  slow  and  fast  fading  The  Rayleigh  fading  ~  0.184  Hz  is 
measured  around  the  median  line  is  not  the  best  model  for  this  non- stationary  fading.  However,  the  lower 
panel  shows  the  unwrapped  phase  and  an  excellent  fit  of  the  quadratic  phase  model 

Zh0(t)  «  0O (t)  :=  0o  +  01 1  +  02t2. 


The  linear  term  estimates  the  Doppler  shift  reported  on  this  plot: 


fo  = 


1  #0  /  ,x 

27T  dt  [) 


t= 0 


In  contrast,  Figure  149  shows  a  deep  fade  in  the  amplitude  1-Hop  X-Mode  channel  and  an  associated 
phase  jump.  The  deep  fade  degrades  the  associated  fading  statistics.  Therefore,  estimates  of  the  Doppler 
spread  will  omit  this  deep  fade. 
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10-Fcb-2016  13:41:00  PST:  O-Mode 


Figure  148.  1-Hop  O-Mode  amplitude  and  phase  at  13:41  PST. 


10-Fcb-2016  13:41:00  PST:  X-Mode 


Figure  149.  1-Hop  X-Mode  amplitude  and  phase  at  13:41  PST. 
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The  last  computation  of  this  section  estimates  the  Doppler  spread.  The  Doppler  spread  is  extracted 
from  O-  and  X-Mode  phases  using  Woodward’s  approach  [57].  Figures  150  and  151  estimate  O-  and  X- 
Mode  Doppler  spectra  from  the  PDF  of  the  instantaneous  frequencies: 

O-Mode/;  (f )  =  T  d^9-  (*) ;  X-Mode/;  (t)  =  L  ^  (f ) . 

Both  figures  compare  the  histograms  of  these  instantaneous  frequencies  against  the  Gaussian  PDF  and  the 
t  Location- Stable  PDF.  The  phase  jumps  in  the  X-Mode  are  handled  by  omitting  the  instantaneous 
frequency  samples  when  | hx(t)  \  <  —  65  dB.  This  omission  is  marked  by  the  in  Table  21.  The  large 
tails  of  instantaneous  frequency  bias  the  estimation  of  mean  and  standard  deviation  to  the  Gaussian  fit. 
Nevertheless,  the  standard  deviation  of  the  Gaussian  determines  the  reported  Doppler  spread  A /  =  2  x  a 
as  discussed  in  Section  4.1. 


Figure  150.  O-Mode  Doppler  spectrum  Figure  151.  X-Mode  Doppler  spectrum 

estimate  at  13:41  PST.  estimate  at  13:41  PST. 


Table  21  reports  these  channel  estimates  in  a  form  suitable  for  a  channel  emulator.  The  model  employs 
only  the  1-Hop  O-  and  X-Modes.  The  total  path  loss  is  scaled  to  the  1-Hop  O-Mode  and  encompasses  both 
the  transmit  and  receive  antennas,  including  the  latter’s  orientation  and  velocity  on  the  mobile.  The 
positive  Doppler  shifts  are  consistent  with  the  mobile  receiver  traveling  toward  the  transmitter. 


Table  21.  Wideband  HF-NVIS  channel  model;  5.82  MHz;  10  February  2016  13:41  PST. 


13:41  PST 

Path 

Total  path 

Doppler 

Doppler 

Fade 

Rician  K 

delay 

loss 

shift 

spread 

rate 

factor 

1-Hop  O-Mode 

1.67  ms 

0  dB 

-0.02  Hz 

0.38  Hz 

0.18Hz 

2  dB 

1-Hop  X-Mode 

1.79  ms 

-3.9  dB 

-0.07*  Hz 

0.33  Hz 

0.36  Hz 

-0.2  dB 
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7.6  MOBILE  WIDEBAND  HF-NVIS  CHANNEL  MODEL:  NON-TID  REGION 


The  preceding  sections  extracted  HF-NVIS  channel  models  at  specific  locations.  This  section  extracts 
over  multiple  locations  to  produce  representative  wideband  HF-NVIS  models  suitable  for  HF  channel 
simulation  and  emulation.  Referring  to  Figure  87,  the  strong  O-  and  X-Modes  show  that  the  HF-NVIS 
channel  is  well-modeled  using  the  1-Hop  channel  of  Example  7 

h(t,  t )  =  aF:o(t)5(T  -  tF;o)  +  af:x(t)S(T  -  tFjX). 

Figure  152  plots  the  estimated  1-Hop  delays  over  the  entire  measurement  day.  The  oscillations  between 
10:30  to  13:00  PST  are  likely  show  the  presence  of  a  TID.  A  relatively  stationary  non-TID  channel  exists 
between  09:30  to  10:30  PST.  This  section  determines  a  stochastic  ensemble  parameterizing  these 
HF-NVIS  channels  in  this  stationary  and  non-TID  region. 


Figure  152.  1-Hop  O-  and  X-Mode  delays. 


Table  22  summarizes  this  1-Hop  model.  The  average  channel  parameters  are  listed  for  comparison  to 
the  narrowband  HF-NVIS  models  of  Table  7.  The  following  subsections  document  the  probabilistic 
models  used  in  this  model. 


Table  22.  Mid-Latitude  mobile  wideband  HF-NVIS  channel  model  using  only  the  1-Hop  O-  and  X-Modes; 
5.82  MHz  over  the  non-TID  region  from  9:30  to  10:30  PST;  interference  and  mode  correlation  are  omitted 
from  this  model. 


1-Hop:  No  TID 
5.82  MHz 

Delay 

(ms) 

Amplitude 

H=10log10(|/j|2) 

Rician  K  factor 
(9:30  to  10:00  PST) 

Doppler  shift 
(Hz) 

Doppler  spread 
(Hz) 

O-Mode 

1.72  ms 

To-lognormal 
p  =  0.540 
a  =  0.011 

/-/-Gaussian 
p  =  -44  dB 
a  =  2.4  dB 

K0-Gaussian 

13  dB 

3.4  dB 

/b-Gaussian 
p  =  0  Hz 
a  =  0.13  Hz 

A/b~exponential 
p  =  0.44  Hz 

X-Mode 

1 .84  ms 

*x=*o  + 
iD-Weibull 

A  =  0.131 

B  =  9.25 

H-Gaussian 
p  =  -50  dB 
a  =  3.4  dB 

Kx~Gaussian 

12  dB 

2.8  dB 

/x~Gaussian 
p  =  0  Hz 
a  =  0.12  Hz 

Afx~exponential 
p  =  0.38  Hz 
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7.6.1  Delay 


Channel  emulators  always  require  the  delays.  Figure  153  is  a  histogram  of  the  joint  delays  tq  and  tx 
between  9:30  to  10:30  PST.  The  conditional  dependence  is  approximated  by  modeling  the  O-Mode  delay 
to  as  a  lognormal  random  variable  based  on  Figure  154  and  the  X-Mode  delay  rx  conditioned  on  to  based 
on  Figure  155. 

to  ~  lognormal(/x,  cr), 
rx  =  t0  +  td;  to  ~  Weibull(^4,  B ) 


Figure  153.  Joint  distribution  of  1-Hop  delays  from  09:30  to  10:30  PST. 

O-Mode  delay  TQ(t):  9.5  <  t  <  10.5  PST  1-Hop  delay  difference  rD(f):  9.5  <  t  <  10.5  PST 


Figure  1 54.  O-Mode  delay  r0  and  fitted  Figure  1 55.  Excess  delay  rD  and  fitted 
distributions.  distributions. 

7.6.2  Fading 

Channel  emulators  associate  a  mean  power  level  and  a  fading  process  with  each  delay.  Rician  fading  is 
the  assumed  model  for  both  the  O-  and  the  X-Mode.  Rician  estimation  splits  the  measured  envelope  power 
into  specular  and  the  scattered  components  [76,  Section  1.2.2]: 

E[\h0(t)\2}  =  s0  +  2cto, 

E[\hx(t)\2}  =  sx  +  2a2x. 

Figures  156  and  157  plot  the  measured  envelope  power  and  the  scattered  powers,  assuming  Rician 
fading. 


106 


Figure  156.  0-Mode  envelope  and 
scattered  power. 


Figure  157.  X-Mode  envelope  and 
scattered  power. 


The  10-dB  gap  between  envelopes  and  their  scattered  powers  argue  that  Rician  fading  is  a  credible 
model  of  both  modes  between  09:30  and  10:30  PST.  However,  the  scattered  power  closes  up  to  the 
envelope  power  starting  around  10:00  PST.  Therefore,  the  table  reports  the  envelope  power — a  nearly 
direct  measurement — over  the  entire  interval  but  reports  the  Rician  K  factors 

Ko  =  10  x  log10  and  Ax  =  10  x  log10 

only  over  a  portion  of  the  non-TID  region  where  the  K  factors  are  approximately  12  dB.  Small  K  factors 
call  into  question  the  Rician  assumption.  Therefore,  modeling  a  non- stationary  K  factor  over  the  entire 
experiment  is  one  approach  to  assess  this  Rician  fading  assumption  [48],  [78]. 

7.6.3  Doppler 

Channel  emulators  typically  specify  a  Doppler  shift  and  a  Doppler  spread.  Figures  158  and  159 
displays  estimates  of  the  Doppler  shift  and  spread  for  the  O-  and  X-Modes  in  the  no-TID  region.  The  large 
spike  at  9.9  PST  is  likely  caused  by  interference  by  the  radio  tower  near  Desert  View  Point  described  in 
Section  6.3.  The  small  values  of  the  Doppler  shift  agree  with  the  narrowband  model  of  Table  7  and  the 
values  predicted  by  the  3-D  ray-tracing  code.  Likewise,  the  mean  values  of  the  exponential  distribution 
fitting  the  Doppler  spreads  align  with  Table  7. 


Figure  158.  O-Mode  Doppler. 


Figure  159.  X-Mode  Doppler. 
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7.6.4  Mode  Correlation 


Channel  emulators  typically  allow  correlation  between  taps.  Correlated  taps  challenge  radio  design 
because  the  ray  paths  will  experience  deep  fades  at  the  same  time.  Figure  160  is  an  estimate  of  the  cross 
correlation  between  the  two  modes.  If  ho(t ,  r)  and  hx(t,  r)  are  the  delay-spread  functions,  the  cross 
correlation 


p(t,  t ) 


Rox(t,T) 

/RooM)/R  xx(t,  0) 


is  the  normalized  cross  covariance  of  the  centered  random  process: 


Rox(t,r)\  :=  E 


(h0(t,  t  +  t') 


Mo(*))  {hx(t,T')  -  nx(t))  , 


where  /jlq (t)  :=  E[ho(t,  r)]  and  ~  E[hx(t ,  r)].  Both  the  expectation  and  associated  estimators 

assume  jointly  stationary  random  processes.  The  plot  shows  that  the  correlation  between  the  estimated 
channels  is  likely  caused  by  the  interference  reported  in  Section  6.3.  In  the  interference-free  regions,  the 
O-  and  X-Modes  show  low  correlation — on  this  day,  at  this  frequency,  and  along  this  route.  Therefore, 
mode  correlation  is  omitted  from  Table  22. 
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8.  EXTENDING  THE  WIDEBAND  HF-NVIS  CHANNEL  ESTIMATES 

Comparison  between  the  measured  channel  and  the  3-D  ray-tracing  code  verified  that  credible  ray 
paths  could  be  simulated  assuming  only  the  the  hybrid  ionogram  and  the  dry  ground  reflections. 
Consequently,  wideband  channel  estimates  should  be  possible  at  other  frequencies,  other  times  of  day, 
other  seasons,  and  other  locations — assuming  that  the  real-world  plasma  frequency  and  ground  reflections 
are  well-modeled  by  the  hybrid  ionogram  and  the  simulated  ground  reflections.  The  following  examples 
use  the  3-D  ray-tracing  code  to  explore  some  of  these  channels  and  implications  for  HF  system  design. 

8.1  O-  AND  X-MODE  ARRIVAL  TIMES 

These  measured  wideband  HF-NVIS  channels  always  show  the  O-Mode  arriving  before  the  X-Mode. 
Therefore,  a  single-antenna  wideband  HF  receiver  could  separate  these  taps  to  mitigate  the  Faraday  fading 
that  plagues  the  narrowband  HF  radios,  provided  the  delay  difference  is  resolvable  in  the  bandwidth. 
Therefore,  the  3-D  ray-tracing  code  simulates  these  delays  as  a  function  of  frequency.  The  ionogram  at 
09:41  PST  is  selected  to  sweep  the  simulated  frequencies  because  of  the  excellent  agreement  with  the 
measured  delays  at  5.82  MHz.  Table  23  lists  the  measured  delays  at  5.82  MHz  and  the  simulated  delays  at 
5.82,  4,  and  3  MHz.  Table  23  predicts  that  the  O-  and  X-Modes  cannot  be  separated  in  delay  at  some 
frequencies.  At  5.82  MHz,  the  measured  O-Modes  arrived  before  the  X-Modes.  This  ordering  holds  at  4 
MHz,  but  at  3  MHz,  the  ordering  reverses — the  X-Modes  arrive  before  the  O-Modes.  Consequently, 
wideband  HF  signal  processing  seeking  to  mitigate  Faraday  fading  by  resolving  these  modes  using  delay 
only  will  fail. 

Table  23.  Measured  and  simulated  hop  times  as  a  function  of  frequency  using  the  ionogram  at  09:41  PST. 


09:41  PST 

Measured 
5.82  MHz 

Simulation 
5.82  MHz 

Simulation 

4  MHz 

Simulation 

3  MHz 

1-Hop  O-Mode 

1.71  ms 

1.71  ms 

1.36  ms 

0.82  ms 

1-Hop  X-Mode 

1.83  ms 

1.85  ms 

1.71  ms 

0.75  ms 

2-Hop  O-Mode 

3.37  ms 

3.34  ms 

2.70  ms 

1 .66  ms 

2-Hop  X-Mode 

3.64  ms 

3.65  ms 

3.28  ms 

1.43  ms 

These  simulations  highlight  the  payoff  of  equipping  the  receiver  with  dual-polarization  antennas: 
Faraday  fading  will  be  mitigated  regardless  of  mode  crossing — provided  that  O-  and  X-Modes  carry 
different  polarizations.  For  this  reason,  the  analysis  of  the  channel  soundings  displayed  the  evolution  of  the 
polarization  along  each  mode.  Indeed,  the  3-D  ray-tracing  code  shows  that  the  O-Modes  start  with  LCP 
and  arrive  as  RCP.  Likewise,  the  X-Modes  start  with  RCP  and  arrive  as  LCP.  This  same  polarization 
behavior  was  observed  for  the  simulations  in  Table  23.  Therefore,  polarization  propagation  and  its  effect 
on  delays  and  fading  must  always  be  part  of  HF  link  analysis. 
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8.2  OVER  THE  WATER 


This  section  simulates  an  HF-NVIS  channel  operating  over  the  water.  This  example  supports 
assessments  of  stealthy  HF  communications  and  makes  explicit  the  shortfalls  of  using  the  standard  HF 
channel  modeling  versus  the  coherent  modeling  delivered  by  the  3-D  ray-tracing  code.  Figure  161  shows 
the  link  between  the  transmitter  on  Point  Loma  and  a  receiver  located  on  San  Clemente  Island.  Table  24 
compares  the  hop  times  extracted  from  the  nearest  ionogram  of  Figure  162  and  the  simulated  rays  of 
Figure  163. 


Figure  161.  Over-the-water  link  between  Point  Loma  and  San  Clemente  Island. 


Table  24.  Comparing  simulated  ray  path  delays  to  the  delays  estimated  from  the  Point  Arguello,  CA 
ionogram  at  15  December  201 5  at  10:00  PST. 


Ray  paths  at 
5.82  MHz 

Simulation 

delays 

Ionogram 

delays 

1-Hop  O-Mode 

1.77  ms 

1 .73  ms 

1-Hop  X-Mode 

1.98  ms 

1.80  ms 

2-Hop  O-Mode 

3.47  ms 

3.47  ms 

2-Hop  X-Mode 

3.90  ms 

3.63  ms 

3-Hop  O-Mode 

5.18  ms 

— 

3-Hop  X-Mode 

5.82  ms 

— 
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Figure  162.  Point  Arguello,  CA  ionogram  at  15  December  2015  at  10:00  PST.  Open  source  ionograms 
courtesy  of  the  Global  Ionospheric  Radio  Observatory  at  the  University  of  Massachusetts,  Lowell  [64], 


Geographic  Longitude  Degrees) 
Geographic  Latitude  (Degrees) 


Geographic  Latitude  Degrees) 


Geographic  Longitude  D^Erees) 


Geographic  Longitude  •Degree;!  Geographic  Longitude  •Degrees;! 


Figure  163.  1-,  2-,  and  3-Hop  between  Point  Loma  and  San  Clemente  Island,  CA;  15  December  2015  at 
10:00  PST. 


Ill 


The  transmitter  broadcasts  from  a  horizontal  dipole  while  the  receiver  employs  a  vertical  loop  made 
from  copper.  Table  25  details  the  simulated  antenna  parameters.  These  antennas  coherently  couple  into  the 
O-  and  X-Mode  ray  paths. 


Table  25.  Simulated  transmitter  and  receiver  antennas. 


TX:  horizontal  dipole 

RX:  vertical  loop 

Orientation 

12° 

Orientation 

45° 

Length 

25.2984  m 

Loop  Radius 

0.5  m 

Height 

12.192m 

Wire  Radius 

0.00762  m 

Figure  164  reports  on  the  over-the-water  channel  simulation.  The  polarization  switches  on  both  modes. 
The  O-Mode  delays  are  to,  too ,  and  tooo-  Likewise  for  the  X-Mode  delays.  As  a  communications 
channel,  the  1-Hop  ray  paths  show  that  the  coherent  processing  on  the  1-Hop  O-  and  X-Modes  will  boost 
receiver  diversity. 


Point  Loma  to  San  Clemente 


5.82  MHz:  15-DEC-2015  10:00  PST 


m-sec:  Modified  Absorption 

Figure  164.  “Channel”  between  between  Point  Loma  and  San  Clemente  Island;  15  December  2015  at 
10:00  PST. 

Table  26  shows  this  channel  is  governed  by  the  absorption.  The  reflection  off  the  sea  water  is 
negligible  and  the  spreading  loss  only  picks  up  4  to  5  dB  per  hop.  However,  absorption  doubles  in  dB  for 
each  hop.  Therefore,  low  absorption  promises  a  multi-tap  channel;  high  absorption  reduces  this  HF-NVIS 
channel  to  1-Hop  O-  and  X-Modes.  The  absorption  also  affects  how  this  HF-NVIS  link  is  detected  by  a 
more  distant  receiver.  Low  absorption  permits  greater  propagation  of  the  1-Hop  O-Mode 
power — assuming  this  particular  frequency  and  propagation  conditions — compromising  the  “stealthiness” 
of  the  HF-NVIS  link.  Therefore,  this  3-D  ray-tracing  code  coupled  with  the  semi-empirical  absorption 
models  can  bound  the  throughput  versus  probability-of-detection  for  such  HF-NVIS  over-the-water  links. 
The  table  shows  that  absorption  is  the  limiting  factor  in  successive  hops.  Indeed,  absorption  scales  by  the 
number  of  hops  in  dB.  Therefore,  absorption  modeling  is  significant  only  when  absorption  is  small — high 
absorption  shuts  down  multiple  hops.  Using  this  reduced  absorption  model,  this  simulation  shows  stealthy 
HF-NVIS  is  still  viable,  but  interference  from  the  2-Hop  O-Mode  requires  consideration. 
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Table  26.  Simulated  propagation  losses  between  between  Point  Loma  and  San  Clemente  Island;  15 
December  201 5  at  1 0:00  PST. 


5.82  MHz 

Delay 

Total  loss 

Spreading 

Mode  shedding 

Absorption 

Reflection 

1-Hop  O-Mode 

1.77  ms 

108.6  dB 

99.6  dB 

3.1  dB 

5.9  dB 

N/A 

1-Hop  X-Mode 

1.92  ms 

1 14.9  dB 

99.1  dB 

3.1  dB 

12.7  dB 

N/A 

2-Hop  O-Mode 

3.45  ms 

120.2  dB 

105.4  dB 

3.2  dB 

11.5  dB 

0.10 

2-Hop  X-Mode 

3.76  ms 

133.1  dB 

104.8  dB 

3.3  dB 

24.9  dB 

0.10 

3-Hop  O-Mode 

5.15  ms 

129.6  dB 

108.9  dB 

3.3  dB 

17.1  dB 

0.21 

3-Hop  X-Mode 

5.62  ms 

149.2  dB 

108.3  dB 

3.5  dB 

37.3  dB 

0.21 

8.3  NETWORKING  FREQUENCIES 

Frequency  selection  is  a  basic  task  of  any  HF  link  or  HF  network.  The  3-D  ray-tracing  code  can 
estimate  link  quality  at  frequencies  without  having  to  sound  the  channel — assuming  a  moderate  absorption 
model.  For  this  example,  the  antennas  are  located  according  to  the  09:41  PST  scenario  and  sized  as  in 
Table  27.  The  ray-tracing  code  uses  the  hybrid  ionogram  obtained  by  spatially  averaging  the  ionograms  at 
Point  Arguello,  CA,  Boulder,  CO,  and  Austin,  TX.  The  end-to-end  simulated  received  power  is  scaled  to 
the  simulated  power  at  5.82  MHz.  Figure  165  shows  this  relative  received  power  as  a  function  of 
frequency.  Software  tools  such  as  Voice  of  America  Coverage  Analysis  Program  (VOACAP)  can  also 
compute  the  total  channel  loss  including  propagation  and  antenna  patterns,  but  this  3-D  ray-tracing  code  is 
unique  by  including  both  mode  and  polarization  effects. 

Table  27.  Antennas  on  the  09:41  PST  link. 


TX:  horizontal  dipole  RX:  vertical  loop 


Orientation 

12° 

Orientation 

45° 

Length 

25.2984  m 

Loop  Radius 

0.5  m 

Height 

12.192m 

Wire  Radius 

0.00762  m 

Latitude 

32°  41 1  54.97"  N 

Latitude 

32°  41’  21.84"  N 

Longitude 

117°  14'  52.96"  W 

Longitude 

116°  18'  27.36"  W 

Figure  165.  Relative  received  power  over  the  09:41  PST  link. 
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8.4  SUMMARY 


The  report  started  with  Sklar’s  observation  that  radio  design  starts  from  channel  characteristics  and 
system  requirements.  The  measurements  and  analysis  show  that  these  mid-latitude  daytime  HF-NVIS 
channels  are  characterized  by  strong  1-Hop  O-  and  X-Modes  arriving  within  0.2  ms  of  one  another  with 
opposite  polarization.  These  extended  simulations  generalize  this  channel  characterization.  Section  8.1 
demonstrates  that  these  modes  arrive  in  mixed  order.  Section  8.2  points  out  that  the  2-Hop  modes  deliver 
additional  diversity  when  operating  over  the  water.  Section  8.3  shows  that  frequency  selection  is  crucial  to 
link  quality.  These  channels  point  to  a  radio  design  that  exploits  both  the  O-  and  X-Modes  rather  than 
suffer  from  Faraday  fading.  A  single-antenna  wideband  system  requires  both  modes  to  be  separated  in 
time.  Dual-polarization  systems  drop  the  time  separation,  provided  the  modes  are  separated  in 
polarization.  More  robust  designs  require  more  breadth  in  the  HF-NVIS  channels. 
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9.  DEVELOPING  THE  HF-NVIS  PORTFOLIO 

This  report  analyzes  mobile  wideband  HF-NVIS  channel  soundings  at  5.82  MHz  collected  on  10 
February  2016  on  a  west-to-east  route.  Two  other  mobile  wideband  HF-NVIS  experiments  expand  upon 
this  effort: 

•  A  mid-latitude  north-south  mobile  route  on  1  Sept  2016 

•  A  mid-latitude  shore-to-ship  test  during  August  2016 

In  the  west-to-east  experiments,  the  Earth’s  magnetic  field  exerts  maximal  effect  on  the  O-  and 
X-Modes.  The  north-to-south  experiment  provides  a  comparison  when  the  Earth’s  magnetic  field  effect  is 
minimized:  the  arrival  times  of  the  O-  and  X-Modes  should  be  reduced,  the  Faraday  fading  effect  should 
be  more  pronounced,  and  these  two  experiments  bound  the  delay  between  these  modes.  However,  one-day 
experiments  only  capture  part  of  the  ionosphere’s  daily  cycle.  The  week-long  shore-to-ship  experiments 
provide  a  multi-day  evolution  of  the  ionosphere  and  employ  dual-polarization  antennas  to  assess  Faraday 
fading  and  the  possible  SIMO  gains.  However,  these  three  experiments  are  limited  to  the  mid-latitudes,  a 
handful  of  frequencies,  and  at  a  limited  number  of  days  in  a  particular  season.  Recommendations  for 
expanding  this  HF-NVIS  portfolio  and  modeling  capability  conclude  this  report. 

Antenna:  Comparison  of  the  whip  to  the  loop  showed  that  the  upward  pointing  pattern  of  the  loop 
reduced  the  RF  noise  level  by  20  dB  or  more  (Section  6).  This  loop  was  simply  placed  on  the  van  as  shown 
in  Figure  61.  There  were  no  modifications  to  reduce  horizontal  sidelobes  or  get  a  tighter  main  beam.  Both 
modifications  will  reduce  the  RF  noise  or,  equivalently,  boost  the  SNR.  Increased  SNR  has  significant 
implication  for  the  mobile  receiver’s  wideband  HF-NVIS  link  budget  (Section  2.3).  Therefore,  antenna 
pattern  design  boosting  vertical  signal  strength  and  reducing  horizontal  RF  noise  for  a  mobile  platform  is 
an  excellent  HF-NVIS  engineering  task. 

Fading:  Both  the  O-  and  X-Modes  exhibited  non- stationary  envelope  fading  consistent  with  “fast”  and 
“slow”  fading  over  the  58-second  soundings.  Neither  fading  process  is  readily  available  from  ionograms 
sampled  at  15-minute  increments.  The  fading  analysis  simply  fit  the  envelopes  with  the  Rician  distribution 
although  other  distributions  were  tested.  Extracting  a  physics-based  model  of  the  “rough- surface” 
scattering  from  the  ionograms  is  an  excellent  HF-NVIS  topic.  Another  approach  could  use  the  clutter 
spectrum  coupled  with  the  ray-tracing  code  to  establish  fading  on  the  ray  paths  [35].  Fikewise,  the  deep 
fades  in  both  the  O-  and  X-Modes  are  puzzling  given  that  the  phase  is  so  smooth.  One  conjecture  is  that 
ground  effects  may  cause  some  fading  as  the  mobile  travels  though  mountainous  terrains. 

Doppler:  The  Doppler  shifts  in  both  the  O-  and  X-Modes  are  well-modeled  using  the  quadratic  phase 
over  the  58-second  soundings  (Sections  7.1,  7.2).  Automatically  extracting  a  physics-based  model  of  this 
Doppler  shift  from  the  ionograms  is  an  excellent  enhancement  of  the  3-D  ray-tracing  code  [43].  Fikewise, 
adapting  the  Doppler  spread  computations  of  Pickering  [62]  to  the  O-  and  X-Mode  ray  paths  is  excellent 
enhancement  of  the  3-D  ray-tracing  code. 

Correlation:  Although  the  O-  and  the  X-Modes  showed  low  cross  correlation,  this  observation  is 
specific  to  this  particular  day,  this  frequency,  and  the  route  in  the  west-to-east  direction.  Therefore,  a 
physics-based  explanation  is  sought  to  explain  this  correlation.  The  north-to-south  channel  soundings  will 
assess  the  role  the  Earth’s  magnetic  field  has  in  separating  the  modes  on  this  west-to-east  direction. 

Noise:  A  positive  result  of  these  experiments  shows  that  the  HF-NVIS  link  with  a  loop  antenna 
reduces  noise  by  20  dB  in  comparison  with  the  standard  whip.  Therefore,  a  noise  analysis  comparing  all 
three  experiments  incorporating  the  antenna  patterns  (Section  5.1),  the  theoretical  noise  levels 
(Section  2.3),  and  the  measured  noise  levels  (Section  6.2)  will  deliver  wideband  HF-NVIS  noise  models 
for  the  mid-latitudes. 
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Interference:  These  measurements  verified  Bark’s  [5]  observation  that  interference  plagues  the  HF 
link  (Section  6.3).  Indeed,  both  the  delay  spread  and  the  cross  correlation  showed  strong  interference 
sources  along  the  rural  route  and  a  wide  noise  swath  driving  through  El  Centro,  CA.  Currently, 
interference  models  for  the  mobile  wideband  HF-NVIS  channel  do  not  exist.  Therefore,  an  interference 
analysis  comparing  this  experiment  and  the  north-south  experiment  will  deliver  wideband  HF-NVIS 
interference  models  for  this  mobile  link. 

Polarization:  The  evolution  of  the  polarization  along  each  mode  was  predicted  in  1934  by  Booker  [9]. 
One  advance  of  the  3-D  ray-tracing  code  is  the  tracking  the  polarization  along  the  ray  paths  and  the 
interaction  with  the  ground  bounces.  Understanding  the  polarization  state  at  the  receiver  is  critical  if  SIMO 
gains  are  to  be  realized  by  polarization  diversity.  The  dual-polarization  sounding  will  quantify  the  SIMO 
gains  for  this  polarization  diversity.  More  generally,  the  3-D  ray-tracing  code  also  supports  the  propagation 
of  the  Poynting  vector  along  each  mode.  Better  HF  direction  finding  and  vector  sensing  is  enabled  by 
separating  these  Poynting  vectors. 

Absorption:  Although  the  3-D  ray-tracing  code  models  absorption  using  the  Hamiltonian,  the 
parameters  of  the  Hamiltonian  are  estimated  from  the  time-varying  ionosphere.  The  excellent  review  paper 
by  Pederick  and  Carvera  [59]  makes  explicit  the  shortfalls  in  many  of  the  absorption  models.  Incorporating 
their  semi-empirical  model  is  a  current  task  to  upgrade  the  3-D  ray-tracing  code. 

Reciprocity:  Non-reciprocity  implies  that  the  full-duplex  channel  emulators  for  HF  radios  must  use 
different  channel  models  for  each  link.  Comparing  measured  and  simulated  non-reciprocity  effects  is  an 
excellent  task  supporting  HF  networking. 

Finally,  efforts  to  acquire  wideband  HF-NVIS  high-latitude  measurements  are  underway.  Figure  4  and 
the  high-latitude  long-haul  measurements  (Section  4.5)  show  that  this  ionosphere  is  challenging  both  in  its 
static  structures  (e.g.,  polar  walls)  and  the  substantial  time-varying  magnetic  fields.  Maslin’s  observation 
[46,  Section  8.2.3]  that  predictions  can  only  be  as  good  as  the  model  shows  that  multiple  high-latitude 
soundings  are  required  to  encompass  these  channels. 
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A.  CALIBRATION  OF  THE  CHANNEL  SOUNDING  SYSTEM 


A.1  CALIBRATION  FACTOR  FOR  RECEIVED  POWER 

This  calibration  relates  the  levels  of  the  digitized  waveform  on  the  receive  software-defined  radios 
(SDRs)  to  the  actual  incident  power  at  the  receive  antenna  feed.  Two  waveforms  are  used  for  calibration:  a 
pure  tone  and  the  channel  sounding  waveform.  The  pure  tone  calibration  produces  an  adjustment  factor  to 
go  from  digitized  samples  to  incident  power.  The  channel  sounding  waveform  calibration  produces  an 
adjustment  factor  to  go  from  a  correlation  function  peak  to  incident  power.  Because  the  transmit  and 
receive  gain  may  vary  with  center  frequency,  the  calibration  was  performed  at  different  frequencies  over 
2  to  12  MHz. 

A  block  diagram  of  the  calibration  setup  is  shown  in  Figure  A-l.  The  transmitting  radio  is  a  National 
Instruments®  USRP  with  an  LFTX  daughterboard.  Power  amplifiers  that  would  normally  go  after  the 
transmitting  radio  are  bypassed  and  33  dB  of  attenuation  is  instead  added.  The  attenuated  waveforms  go 
into  the  receive  chain.  The  first  receive  chain  blocks  are  a  transient  power  disrupter  and  a  power  limiter. 
Then,  unlike  the  over-the-air  channel  sounding,  the  signal  does  not  encounter  a  180°  hybrid  coupler 
because  the  signal  only  enters  the  receive  chain  from  one  port  instead  of  a  two-port  receive  loop.  The 
hybrid  coupler  was  measured  to  attenuate  each  port  by  about  0.1  dB  over  the  entire  2  to  12  MHz  band. 
Then,  in  this  calibration,  the  signal  resumes  the  same  path  as  with  the  over-the-air  measurements  through  a 
bandpass  filter,  a  low-noise  amplifier  (LNA)  rated  to  43  dB  gain,  and  finally  into  an  RFSpace®  SDR.  This 
receive  SDR  has  internal  attenuation  that  can  be  set  at  either  0,  10,  20,  or  30  dB.  The  rated  receive  SDR 
maximum  input  power  after  internal  attenuation  is  13  dBm.  For  all  calibration  measurements,  the  SDR  had 
30  dB  of  internal  attenuation  to  prevent  exceeding  the  SDR’s  maximum  input  power.  This  30-dB  internal 
attenuation  is  taken  into  account  when  computing  adjustments  to  dBm,  so  all  adjustments  to  dBm  assume 
no  attenuation  inside  the  receive  SDR. 


Measured 

power 


Transient  power 
disruptor 


Figure  A-1 .  Block  diagram  of  the  calibration  using  the  Nl  USRP  transmit  radio  (with  external  power 
amplifiers  bypassed),  external  attenuators,  and  the  receive  chain  and  receive  SDR  (RFSPACE®  SDR-IP). 

Measurements  of  the  transmitted  power  were  made  on  a  Keysight  Fieldfox®  Spectrum  Analyzer  at  the 
point  in  Figure  A-l  denoted  by  “Measured  power.”  The  transmit  radio  power  varies  based  on  the 
bandwidth  of  the  signal,  so  the  powers  of  both  a  pure  tone  and  the  25-kHz  channel  sounding  waveform 
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were  measured.  The  measurement  bandwidths  while  measuring  pure  tones  and  the  sounding  waveform 
were  10  kHz  and  50  kHz,  respectively.  These  bandwidths  were  chosen  so  all  of  the  signal  frequency 
content  was  sure  to  be  included  in  the  power  measurement  without  adding  excessive  internal  noise  in  the 
power  meter.  Table  A-l  gives  the  power  from  pure  tones.  The  adjustment  factors  used  to  go  from  “digital 
power”  10  log10  ( E  [\sr^ [n]\2] )  1°  power  incident  on  the  receive  chain  in  dBm  are  used  as  follows: 

Raw  power  at  receive  chain  (dBm)  =  10  log10  ( E  [n\  |2] )  +  Radio  n  adjustment  (dB), 
where  sr^d[ti\  is  the  digitized  signal  at  the  receive  SDR  and  E[-]  is  the  expectation  operator. 


Table  A-1 .  Pure  tone  power  measurements  and  adjustment  factors  for  each  radio  to  go  from  “digital 
power”  to  incident  power  at  the  receive  chain.  Not  included  in  the  adjustment  factors  is  the  loss  through  a 
hybrid  coupler,  which  was  measured  to  be  approximately  0.1  dB  over  the  entire  HF  band. 


Freq.  (MHz) 

Measured  pwr. 
(dBm) 

Radiol:  101og10  (e  [H«^>[n]|2]) 

Radio  2:  101og10  (e  [|^j>  [n] |2]) 

Radio  1  adjustment  to  dBm 
(dB) 

Radio  2  adjustment  to  dBm 
(dB) 

2.000 

-29.03 

-12.72 

-21.55 

-46.31 

-37.48 

2.387 

-29.03 

-11.35 

-20.45 

47.68 

-38.58 

5.161 

-29.14 

-7.44 

-17.04 

-51.70 

42.10 

5.820 

-29.16 

-7.21 

-16.78 

-51.95 

42.38 

12.000 

-29.86 

-6.15 

-15.56 

-53.71 

44.30 

A  similar  procedure  determines  calibration  factors  to  go  from  the  peaks  in  the  correlation  of  the 
received  signal  with  the  channel  sounding  waveform  and  the  actual  incident  power  of  the  channel  sounding 
waveform  at  the  receiver.  The  channel  sounding  waveform  is  a  BPSK-modulated  pseudo-noise  code  (PN 
code  10)  sequence  with  an  autocorrelation  that  is  nearly  a  delta  function.  The  autocorrelation  is  not  a 
perfect  delta  function  due  to  the  root-raised-cosine  filtering  of  the  transmitted  waveform.  The  channel 
sounding  waveform  was  chosen  so  that  the  separate  propagation  paths  between  the  transmitter  and  receiver 
are  denoted  by  distinct  peaks  after  the  received  signal  is  correlated  with  the  transmitted  waveform.  To 
compute  the  loss  through  the  channel,  the  incident  power  of  the  signal  at  the  receiver  from  these 
propagation  paths  must  be  calculated.  Consequently,  if  the  transmit  power  delivered  to  the  antenna  and  the 
transmit  and  receive  antenna  realized  gains  (Institute  of  Electrical  and  Electronics  Engineers  [IEEE] 
Standard  145-1993)  are  known,  the  attenuation  of  each  peak  or  propagation  path  is  computable. 

The  channel  sounding  waveform  begins  each  second  on  the  second  with  22  PN  10  sequences  sent  at  a 
rate  of  25,000  symbols  per  second.  Because  the  PN  10  sequence  is  1023  bits  long,  the  channel  sounding 
waveform  has  a  duration  of  0.9002  seconds.  The  remainder  of  each  second  is  followed  by  silence.  Because 
the  correlation  to  determine  channel  power  only  is  concerned  with  the  times  the  waveform  is  received  and 
not  the  silent  periods,  the  calibration  power  measurements  do  not  include  a  silent  period  and  rather 
continuously  repeat  the  PN  10  sequence. 

Channel  correlation  is  calculated  using  the  procedure  outlined  in  Section  5.3.  Because  in  this 
calibration  there  is  a  6-ft  cable  connecting  the  transmit  SDR  to  the  external  attenuators  and  receive  chain, 
the  correlation  should  be  a  single  peak  with  a  delay  of  zero  because  the  propagation  time  is  much  faster 
than  the  sampling  rate.  An  example  correlation  peak  from  Radio  1  for  a  2-MHz  center  frequency  is  shown 
in  Figure  A-2  and  is  calculated  by  20  log10  (\Rrc(t)  |),  where  Rrc(t )  is  the  correlation  of  the  received 
signal  with  the  transmitted  signal.  All  correlations  have  this  similar  shape  with  peaks  given  in  Table  A-2. 

The  adjustment  factors  are  used  to  go  from  correlation  peak  20  log10  (\Rrc(t)  |)  to  power  incident  on 
the  receive  chain  in  dBm  as  follows: 

Sounding  waveform  power  at  receive  chain  (dBm)  =  20  log10  (\Rrc(t)\)  +  Radio  n  adjustment  (dB) 
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Figure  A-2.  Channel  correlation  from  calibration  setup  with  a  center  frequency  of  2  MHz. 

Table  A-2.  Channel  correlation  peak  values  and  adjustment  factors  for  each  radio  to  compute  incident 
power  at  the  receive  chain  of  the  channel  sounding  waveform.  Correlations  are  scaled  by  a  normalization 
factor  of  0.025133.  Not  included  in  the  adjustment  factors  is  the  loss  through  a  hybrid  coupler,  which  was 
measured  to  be  approximately  0.1  dB  over  the  entire  HF  band. 


Freq.  (MHz) 

Measured  pwr. 
(dBm) 

Radio  1:201og10(|fli,c(T)|) 

Radio  2:  20 log10( |U*c(t)|) 

Radio  1  adjustment  to  dBm 
(dB) 

Radio  2  adjustment  to  dBm 
(dB) 

2.000 

-38.04 

16.24 

8.78 

-84.28 

-76.82 

2.387 

-38.07 

17.62 

9.84 

-85.69 

-77.91 

5.161 

-38.27 

21.62 

13.34 

-89.89 

-81.61 

5.820 

-38.25 

21.87 

13.94 

-90.12 

-82.19 

12.000 

-38.91 

22.98 

15.13 

-91.89 

-84.04 

The  digitized  power  from  a  tone  was  on  average  7.9  dB  less  at  Radio  2  than  at  Radio  1.  The  digitized 
correlation  peak  was  on  average  9.3  dB  less  at  Radio  2  than  at  Radio  1.  Secondary  peaks  appear  in  the 
channel  correlations  when  processing  the  Radio  2  calibration  data  but  not  the  Radio  1  calibration  data. 
These  peaks  occurred  at  a  time  delay  greater  than  2  ms,  but  the  only  correlation  peak  should  be  at  0  ms 
delay.  These  secondary  peaks  suggest  the  signal  at  Radio  2  may  have  exceeded  the  radio’s  dynamic  range, 
as  will  be  explained  in  Section  A.2.  However,  the  over-voltage  indicator  on  Radio  2  was  negative  during 
the  calibration.  To  be  conservative,  the  received  power  at  Radio  2  may  be  accurate  only  to  within  2  dB. 

The  receive  chain  as  configured  in  Figure  A-l  represents  the  receive  chain  for  channel  soundings  after 
February  2016.  For  the  channel  sounding  done  in  February  2016,  the  receive  chains  did  not  include  power 
limiters  or  bandpass  filters.  However,  the  calibration  was  remeasured  after  removing  the  power  limiters 
and  bandpass  filter  in  Figure  A-l  and  the  correlation  peaks  differed  by  less  than  0.1  dB  at  each  frequency. 

A.2  EFFECT  OF  RECEIVE  SDR  INTERNAL  ATTENUATION 

The  receive  SDR’s  internal  attenuation  can  be  adjusted  by  10-dB  increments  between  0  dB  and  30  dB. 
As  expected,  increasing  the  internal  attenuation  10  dB  scales  the  digitized  signal  by  about  10  dB. 
Therefore,  the  attenuation  level  must  be  taken  into  account  when  computing  the  channel.  The  next  set  of 
figures  show  the  correlation  peaks  for  different  levels  of  SDR  internal  attenuation.  Figure  A-3  shows  the 
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correlation  when  no  internal  attenuation  is  added  to  the  radio.  The  additional  spikes  in  the  correlation  are 
caused  by  the  received  signal  exceeding  the  dynamic  range  of  the  SDR’s  digitizer.  The  peak  value  should 
not  be  taken  as  reliable  due  to  this  error  in  quantization. 


Figure  A-3.  Channel  correlation  with  no  internal  attenuation.  The  additional  spikes  are  due  to  the  received 
signal  constantly  exceeding  the  dynamic  range  of  the  SDR’s  digitizer. 

Figures  A-4,  A-5,  and  A-6  show  the  channel  correlations  when  the  receive  SDR  has  10,  20  and  30  dB 
of  internal  attenuation,  respectively.  The  level  of  the  maximum  correlation  decreases  10  dB  with  less  than 
0.1  dB  of  error  each  time  the  attenuation  is  increased  10  dB.  The  absolute  level  of  these  correlation  peaks 
should  not  be  compared  directly  to  Section  A.l  because  different  attenuation  was  used  between  transmit 
and  receive  chains  when  each  of  these  measurements  was  taken. 


Figure  A-4.  Channel  correlation  with  10-dB  internal  attenuation. 

A  MATLAB®  function  called  “GetlnternalAttenuation.m”  extracts  the  radio’s  internal  attenuator 
setting  in  the  metadata  of  a  .wav  file  recorded  in  Spectra  Vue®.  The  setting  is  encoded  in  the  .wav  file 
using  the  Resource  Interchange  File  Format  (RIFF).  The  RIFF  subchunk  containing  the  attenuation  is 
labeled  “auxi”  and  the  attenuation  is  a  32-bit  integer  after  sixteen  16-bit  integers  and  five  32-bit  integers  in 
this  subchunk.  This  is  documented  on  page  64  of  the  Spectra  Vue®  manual,  available  at 
http : / /www .moetronix . com/ files/ spectravue .pdf. 
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Figure  A-5.  Channel  correlation  with  20-dB  internal  attenuation. 


Figure  A-6.  Channel  correlation  with  30-dB  internal  attenuation. 
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B.  DOPPLER  SHIFT  MEASUREMENTS 


This  section  undertakes  Doppler  shift  estimation  without  vehicle  motion.  The  implicit  assumption  in 
this  approach  is  that  after  subtracting  the  Doppler  shift  due  to  vehicle  motion,  what  is  left  is  the  Doppler 
shift  produced  by  the  ionosphere.  In  this  context,  the  natural  “test  of  this  approach”  requires  comparison 
against  other  Doppler  measurements,  albeit  limited.  This  section  begins  with  an  estimate  of  the  Doppler 
shift  due  to  vehicle  motion,  then  presents  measured  Doppler  shifts  when  subtracting  the  shift  due  to 
vehicle  motion,  and  finally  fits  these  Doppler  shifts  to  statistical  distributions. 

B.1  DOPPLER  SHIFT  DUE  TO  PLATFORM  MOBILITY 

The  transmitter  and  receiver  platforms  are  ground-based  vehicles.  The  vehicles’  motions  relative  to  the 
direction  of  propagation  induces  a  Doppler  shift  of  the  transmitted  signal.  Because  the  path  taken  for  NVIS 
communication  is  not  the  direct  path  between  the  two  vehicles,  the  Doppler  shift  is  not  simply  a  function 
of  the  relative  velocities  of  the  two  vehicles  but  rather  is  dominated  by  the  angle  of  departure  and  arrival  of 
the  transmitted  signal. 

Assume  two  ground  vehicles  travel  with  velocities  v\  and  v2  in  random  directions,  as  shown  in 
Figure  B-l.  The  transmitted  wave  refracts  off  of  the  ionosphere  in  more  of  a  curved  path  due  to  the 
gradient  in  refraction  index  of  the  ionosphere  as  a  function  of  elevation.  For  simplicity,  assume  a  perfect 
reflection  off  of  a  single  layer  ionosphere  at  some  virtual  height  with  an  angle  of  reflection  Or. 


Figure  B-1.  The  direction  of  propagation  of  an  HF-NVIS  signal  between  two  mobile  ground  vehicles  using 
NVIS  antennas  with  velocities  vx  and  v2. 

The  angle  made  between  the  ground  velocity  of  the  vehicle  v%  and  a  ground  velocity  that  follows  the 
propagation  direction  v[  is  given  by  the  following  equation: 

Vi  =  Vi  COS (0i), 

where  i  is  either  1  or  2  corresponding  to  the  transmitter  and  receiver  in  Figure  B-l.  The  angle  between  the 
ground  direction  velocity  vector  following  propagation  and  the  actual  propagation  velocity  vector  v"  is 
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given  by: 


//  /  fnl\  t  •  (@R 

Vi  =vicos(di)  =^sm(  — 


The  Doppler  shift  caused  by  the  velocities  in  the  direction  of  propagation  of  the  transmitter  and 
receiver  is  given  by: 

}D  = 

c  c  Ao 


where  c  is  the  speed  of  light,  /o  is  the  signal  center  frequency,  and  Ao  is  the  wavelength.  Substituting  to 
express  the  Doppler  shift  as  a  function  of  velocities  and  ground  directions  and  the  ionosphere  reflection 
angle  gives: 


,,  V\  cos(6>i)  -  V2  cos (02 )  . 

Sd  = - 7 - sin 

Ao 


(B-l) 


Simulation  and  measurement  of  angles  of  reflection  for  NVIS  propagation  is  found  in  [84].  In  that 
reference,  simulations  showed  that  launch  elevation  angles  from  65°  to  90°  were  sufficient  for  a  coverage 
area  with  a  radius  of  150  km.  Measurements  of  the  received  angle  of  incidence  for  transmitters  9  to 
165  km  away  sending  a  waveform  at  7  MHz  yielded  a  distribution  of  incidence  angles  between  70°  and 
90°.  Assume  the  angle  of  incidence  is  uniformly  distributed  from  65°  to  90°,  which  means  Or  is  uniformly 
distributed  between  0°  and  50°. 

To  compute  typical  values  for  Doppler  shift,  assume  that  the  ground  angles  0\  and  62  are  uniformly- 
distributed  random  variables  from  0  to  2n.  Assume  the  velocity  is  less  than  70  mph  or  about  32  m/s.  The 
wavelength  Ao  for  frequencies  2  to  30  MHz  ranges  from  150  m  to  10  m. 

Doppler  shifts  are  generated  from  one  hundred  thousand  Monte  Carlo  simulations  using  the  previously 
stated  parameters  at  center  frequencies  from  2  to  12  MHz.  The  Doppler  shifts  are  all  distributed 
symmetrically  around  0  Hz.  The  root  mean  squared  Doppler  shift  for  each  frequency  is  shown  in 
Figure  B-2. 


Figure  B-2.  The  root  mean  squared  Doppler  shift  as  a  function  of  frequency. 

Two  histograms  of  Doppler  shifts  at  center  frequencies  of  2  and  12  MHz,  respectively,  are  shown  in 
Figures  B-3  and  B-4.  These  frequencies  are  seen  as  the  lowest  and  highest  feasible  NVIS  propagation 
frequencies.  The  Doppler  shift  is  likely  to  be  the  highest  at  12  MHz,  but  even  at  this  frequency,  the  shift  is 
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not  likely  to  be  more  than  0.5  Hz.  This  Doppler  shift  will  be  compared  with  the  level  of  Doppler  spread 
induced  by  the  ionosphere  itself  caused  by  the  rising  and  falling  of  its  layers. 


Figure  B-3.  A  histogram  of  the  predicted  Doppler  shifts  when  transmitting  at  2  MHz. 


Figure  B-4.  A  histogram  of  the  predicted  Doppler  shifts  when  transmitting  at  12  MHz. 

B.2  DOPPLER  SHIFT  FROM  MEASURED  CHANNEL  SOUNDING  DATA 

The  use  of  frequency-locked  transmit  and  receive  radios  for  channel  sounding  allows  a  direct 
calculation  of  Doppler  shifts.  Doppler  shifts  are  calculated  from  the  change  in  the  phase  of  a  peak  of  the 
delay-spread  function  over  time.  The  delay-spread  function  is  calculated  by  correlating  the  transmitted 
waveform  with  the  received  signal  [8],  and  the  peaks  of  the  absolute  value  of  this  correlation  correspond  to 
propagation  paths. 

The  distribution  of  Doppler  shifts  is  the  Doppler  spread.  This  section  is  concerned  with  the  Doppler 
spread  of  the  first  hops  of  the  O-  and  X-Modes.  From  the  power  delay  profile  plotted  over  the  entire 
sounding  in  Figure  87,  the  first  hop  delays  can  be  seen  by  inspection  to  lie  between  1.5  ms  and  2.1  ms. 

The  process  of  finding  correlation  peaks  is  automated  according  to  the  following  rules.  The  maximum 
absolute  value  of  the  correlation  function  is  considered  a  peak.  Any  additional  peaks  must  be  located  at  a 
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delay  where  the  absolute  value  of  the  correlation  that  is  larger  than  the  two  neighboring  samples  and  must 
be  greater  than  15  dB  down  from  the  global  peak.  The  15  dB  figure  comes  from  the  fact  that  the  sidelobes 
of  the  autocorrelation  of  the  sounding  waveform  are  about  16  dB  down.  Also,  if  there  are  more  than  four 
peaks  that  meet  these  criteria,  the  maximum  peak  value  is  likely  near  the  noise  floor.  In  this  case,  it  is 
difficult  to  determine  secondary  peaks,  so  only  maximum  peak  is  recorded. 

The  phase  in  radians  of  the  correlation  at  any  peak  is  taken  for  just  under  a  minute,  out  of  convenience 
because  the  receiver  recordings  are  made  into  new  files  each  minute.  The  least-squares  fit  to  the  slope  of 
the  phase  over  time  divided  by  2n  is  the  estimate  of  the  Doppler  shift. 

The  Doppler  shifts  corresponding  to  all  correlation  peaks  are  plotted  in  Figure  B-5.  Of  the  656 
calculated  Doppler  shifts,  there  are  three  points  that  exist  outside  —0.5  to  0.5  Hz  at  —1.7  Hz,  —2.5  Hz  (not 
shown)  and  —3.0  Hz  (not  shown).  The  Doppler  shifts  caused  by  vehicle  motion  are  calculated  using 
Equation  B-l.  The  calculated  Doppler  caused  by  vehicle  motion  is  nearly  the  same  for  the  O-  and 
X-Modes  over  the  same  minute,  because  the  difference  in  delay  only  slightly  changes  the  angle  of 
reflection  off  of  a  virtual  height. 

0.5 

0 


o 

Q 


-1.5 

-2 

Figure  B-5.  Doppler  shifts  of  all  correlation  peaks  throughout  the  El  Centro,  CA  channel  sounding  run. 
Doppler  shifts  of  the  first  hop  O-Mode  are  blue  dots,  the  part  of  the  Doppler  shift  attributable  to  vehicle 
motion  for  the  O-Mode  paths  are  red  dots,  Doppler  shifts  of  the  first  hop  X-Mode  are  green  dots,  and  the 
part  of  the  Doppler  shift  attributable  to  vehicle  motion  for  the  X-Mode  paths  are  magenta  dots 

To  determine  the  Doppler  shift  induced  by  the  ionosphere,  the  Doppler  shift  caused  by  vehicle  motion 
is  subtracted  from  the  total  calculated  Doppler  shift.  The  residual  Doppler  shifts  for  the  O-  and  X-Modes 
are  collected  into  histograms  shown  in  Figure  B-6.  Both  modes  present  mean  Doppler  shifts  of  0.02  Hz  in 
a  univariate  distribution  of  approximately  the  same  variance.  The  O-Mode  variance  around  approximately 
0  Hz  is  slightly  greater.  This  greater  variance  is  probably  because  when  there  is  one  peak  very  close  to  the 
noise  floor,  it  is  assigned  to  the  O-Mode,  and  some  assigned  peaks  very  close  to  the  noise  floor  might  be 
noise  themselves. 

B.3  FITTING  A  DISTRIBUTION  TO  THE  MEASURED  DOPPLER  SPREAD 

Like  ray  path  delays  and  attenuation,  Doppler  shifts  for  the  O-  and  X-Modes  used  in  a  channel 
emulator  are  based  on  a  probability  distribution.  Figure  B-6  compares  the  histograms  of  the  two  Doppler 
shifts.  The  O-Mode  Doppler  shifts  tend  to  be  slightly  lower  than  the  X-Mode  Doppler  shifts  and  both  are 
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slightly  asymmetric. 


Figure  B-6.  A  histogram  of  measured  Doppler  shifts  of  the  first  hop  O-  and  X-Modes  with  Doppler  caused 
by  receive  vehicle  motion  subtracted  off. 


Figures  B-7  and  B-8  assess  the  “goodness-of-ht”  between  the  Gaussian  distribution  and  the  “t 
Location-Scale”  distribution  on  the  O-  and  X-mode  Doppler  shifts.  Table  B-l  compares  the  mean 
(/i  «  average  Doppler  shift),  standard  deviation  (2 a  «  Doppler  spread),  and  the  log-likelihood  [42].  The 
higher  the  log  likelihood  function,  the  better  the  distribution  fits  the  measured  samples.  The  asymmetries 
in  the  measured  data  tend  to  bias  the  Gaussian  estimate  more  than  the  “t  Location-Scale”  distribution. 
These  estimates  agree  with  the  near-zero  Doppler  shifts  estimated  from  the  ionogram  in  Section  4.3,  and 
Doppler  shift  in  the  selected  stops  in  Section  7  that  include  the  vehicle  speed.  If  the  vehicle  speed  only 
shifts  the  Doppler,  the  associated  Doppler  spread  is  unchanged.  The  Doppler  spreads  in  Table  B-l, 
estimated  across  the  entire  experiment,  are  comparable  to  the  Doppler  spreads  reported  in  Table  22 
estimated  only  over  the  non-TID  times.  The  upshot  is  that  these  mobile  HF-NVIS  mid-latitude  rays  show 
near-zero  Doppler  during  the  daytime  hours — the  high  angle  and  relatively  low  mobile  speed  do  not  shift 
the  Doppler  induced  by  the  ionosphere. 


Doppler  shift  (Hz)  Doppler  shift  (Hz) 


Figure  B-7.  Histogram  of  measured  O-Mode 
Doppler  shifts  with  best  fit  Gaussian  and 
t  Location-Scale  distributions. 
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Figure  B-8.  Histogram  of  measured  X-Mode 
Doppler  shifts  with  best  fit  Gaussian  and 
t  Location-Scale  distributions. 


Table  B-1 .  Modeling  the  Doppler  shift  induced  by  the  ionosphere  without  vehicle  speed. 


Ray  path 

a 

Distribution 

Likelihood 

O-Mode 

0.0187  Hz 

0.234  Hz 

Gaussian 

10 

X-Mode 

0.0208  Hz 

0.227  Hz 

Gaussian 

16 

O-Mode 

0.0335  Hz 

0.183  Hz 

t  Location-Scale 

107 

X-Mode 

0.0331  Hz 

0.133  Hz 

t  Location-Scale 

151 
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C.  DOPPLER  SHIFTS  AND  THE  TIME-VARYING  TRANSFER  FUNCTION 


This  section  reviews  the  simple  1-Hop  model  for  delay  and  Doppler  shifts  in  the  HF-NVIS  channel. 
The  model  explains  some  features  observed  in  the  spectrograms  of  the  measured  channel  soundings.  The 
canonical  HF-NVIS  channel  model  consisting  of  only  one  strong  O-  and  X-Mode  has  delay-spread 
function 

h(t,  t)  —  hoe+^27r^ot5(r  —  to)  +  hxe+^27r^xtS(r  —  rx). 

Bello’s  time-varying  transfer  function  links  the  model  to  the  spectrogram.  The  time-varying  transfer 
function  is  Fourier  transform  in  delay  of  the  delay-spread  function  h(t,  r)  [8]: 

/OO 

e-j2wfTh(t,  r)dr. 

-CO 

Assume  the  transmitted  signal  st  G  L1(M)  f]  L2(M)  has  the  spectral  representation 

/oo 

e+^sT(f)dt. 

-oo 

The  time- varying  channel  acts  on  the  transmitted  signal  as 

/OO  /‘OO 

h(t,  t  -  r)sT(r)dr  =  /  e+j2lTftT(t,  f)sT(f)df 
-OO  J  — OO 

and  captures  the  O-  and  X-Modes  as  follows: 

T(t,f)  =  hoe+j27Tfote-j2nfTo  +  hxe+j27rfxte~j2nf T 


if  the  spectrum  of  the  transmitted  signal  is  relatively  flat,  one  estimate  of  the  received  signal’s  spectrogram 
is  |  T(£,  f)s'r(f)  |2-  This  spectrogram  reveals  the  constructive  and  destructive  interference  of  the  O-  and 
X-Modes  as 


\H(t,f)\  =  \h0\ 


1  +  pe 


-i 


where  p  =  hx/ho  and  the  phase  is 


4>(t,  f)  ■=  2 7t(/x  -  fo)t  -  2i r/(rx  -  r0). 


The  canonical  example  of  the  O-  and  X-Mode  interaction  assumes  both  modes  have  equal  power. 


Example  13  (Equal  Modes)  Assume  ho  =  hx  =  1-  In  the  O-Mode,  assume  zero  Doppler  shift  fo  —  0 
Hz  and  time  delay  to  =  2.16  ms.  In  the  X-Mode,  assume  fo  =  0.22-Hz  Doppler  shift  and  time  delay 
rx  —  2.34  ms.  Figure  C-l  plots  iT(t,  /)  over  a  40-kHz  band  for  a  3-second  increment.  For  fixed  time, 
spectrogram  sliced  infrequency  will  have  minima  spaced  at  multiples  of 

A  fN  =  7 - 1 - 7. 

rx  -  to  I 


For  fixed  frequency,  spectrogram  sliced  in  time  will  minima  spaced  at  multiples  of 


A  tN 


1 

l/x  -  fo\ 


The  minima  are  lines  parallel  to 


t  =  -f 


TX  -  TQ 

fx  -  fo 
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seconds:  r  =2.16;  r  =2.34  see 
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kHz:  /o=0;/x=0.22  Hz 

Figure  C-1 .  Time-varying  filter. 


C-2 


D.  3-D  IONOSPHERIC  RAY  TRACING  USING  DRION 


Drlon  is  a  three-dimensional  ionospheric  ray-tracing  program  used  to  analyze  the  channel  performance 
for  HF  communications.  This  program  is  being  developed  in  MATLAB®  based  on  the  Jones-Stephenson 
formulation  of  Hamilton’s  ray-tracing  equations  [36].  Electromagnetic  waves  propagating  through  the 
ionosphere  will  split  into  two  possible  waves.  These  are  called  the  ordinary  (O-Mode)  and  extraordinary 
(X-Mode)  waves.  Drlon  computes  both  of  these  magneto-ionic  waves  for  single  and  multiple  hops 
between  the  transmitter  and  receiver.  Hamilton’s  ray-tracing  equations  written  in  spherical  polar 
coordinates  are  given  by 


dr 
~dP' 
dd 
~dP’ 
d(f> 
dP 7 
dkr 
dP 1 
dke 
dP' 
dk<j) 
dP 7 


1  dH/dkr 
~c  DH/Ou 
_  1  dH/dke 
rc  dH/duj 
1  dH/dka, 
rc  sin  6  dH/duj 

1  dH/dr  ,  d9  ,  .  „  # 

+k8T5i  +  k<t>  smeT5i 
c  oH /  ouj  dP  dP 

l/l  dH/ dd  7dr  ,  „  dr/) 

~  -~auTr> - +  k4>rcoseT5 7 

r  \c  oH / ouj  dP  dP' 

1  (ldH/dt  7  .  ndr  l  nd0 

- ^  517T/1T —  k<t>smd^si  ~  k^rcosd^rsj 

r  sm9  \c  oH  /duo  dPf  dPr 


where  P'  =  ct  is  the  group  path  length  and  H  is  the  user-selected  Hamiltonian.  These  ray-tracing 
equations  are  solved  simultaneously  using  an  ordinary  differential  equation  (ODE)  solver.  The  solution  of 
these  equations  gives  the  spherical  polar  coordinates  of  each  point  along  the  ray  path  (r,  0,  </>)  and  the 
propagation  vector  in  the  wave  normal  direction  (£;r,  k$,  k^).  Other  ray  path  quantities  are  found  once  the 
ray  path  is  computed. 

The  infinitesimal  Doppler  shift  resulting  from  a  time- varying  electron  density  profile  is 

djfp)  1  c)H/dt 

dP'  27 rc  dH/duo 

The  total  Doppler  shift  is  computed  by  integrating  along  the  ray  path. 

The  phase  path  length  is  the  phase  shift  along  the  ray  path  divided  by  the  free  space  wavenumber.  The 
phase  path  length  is  computed  by  integrating  [36,  Eq.  16]: 

]T1  -<  Zp  DH_  _i_  k  dH  ,  k  M. 

dP  _  1  Kr  dkr  +  dke  ^  ^  dk 0 

dP'  ~  dH/duo 

along  the  ray  path.  The  phase  path  length  quantity  is  critical  to  determine  the  wave  polarization  and  the 
polarization  fading  at  the  receiver. 

Power  absorbed  in  the  ionosphere  along  an  infinitesimal  segment  of  the  ray  path  is  [36,  Eq.  17]: 

dA  _  10  Q(^n2)  dP 

dP'  ~~  ln(10)  k°  ¥  ~dPr 

The  total  amount  of  power  (in  dB)  lost  due  electron  collisions  with  neutral  particles  is  the  integral  along 
the  ray  path. 
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The  geometrical  path  length  is  computed  by  integrating  [36,  Eq.  18] 


dP'  c  |  dH/du 


along  the  ray  path. 


Both  polarization 

_  .  Y?  ± 

P  ~  J  2Yl  (U  —  X) 


and  longitudinal  polarization 

pL  =  J  (U  -  X)  ^  ~~  n  ) 

of  the  O-  and  X-Modes  are  computed  along  the  ray  path  [11].  These  polarization  compute  the  magnitudes 
of  the  field  components  along  the  ray  path,  the  wave  polarization  along  the  ray  path,  and  finally,  the  correct 
voltage  induced  on  the  receive  antenna. 


Hamiltonian:  Drlon  includes  several  options  for  the  Hamiltonian.  The  simplest  Hamiltonian  is 


where  k 2  is  the  square  of  the  propagation  vector  calculated  from  the  solution  of  the  ray-tracing  equations 
and  n2  is  the  square  of  the  refractive  index  of  the  ionosphere.  This  Hamiltonian  can  be  used  with  either 
the  Appleton-Hartree  formula  [11]  or  the  Sen-Wyller  formula  [68]  for  the  index  of  refraction  and  works 
well  for  most  situations.  However,  the  ray-tracing  solution  breaks  down  when  the  ray  path  goes  through  a 
spitze.  A  spitze  occurs  when  a  ray  path  never  becomes  horizontal,  but  has  a  cusp  at  the  height  where  X  — 
1  [11].  In  this  case,  a  Hamiltonian  based  on  the  Booker  Quartic  formula  is  required.  The  solution  of  the 
Booker  Quartic  equation  is  the  Appleton-Hartree  formula  for  the  index  of  refraction.  The  Hamiltonian 
based  on  the  Booker  Quartic, 


H  =  sft  {  [(17  -  X)  U2  -  Y2U]  cAkA  +  x(k-Yj2  cAk2+ 


-2 U  (U  -  X)2  +  Y2  (2 U  -  X) 
( U-X)2-Y 2]  ( U-X)uA } 


c2k2uj2  -  X 


2  2  i 

C  OJ  + 


is  valid  in  all  ray-tracing  cases,  including  ray  paths  near  a  spitze.  A  similar  Hamiltonian  included  in  Drlon 
uses  a  quartic  equation  whose  solution  is  the  Sen-Wyller  mono-energetic  index  of  refraction. 


Improvements:  The  Hamiltonian  encoded  in  Drlon  supports  substantial  improvements.  First, 
improvements  and  validation  of  the  absorption  calculation  needs  to  be  performed.  The  absorption  is 
difficult  to  compute  since  it  depends  on  accurate  knowledge  of  the  electron  density,  neutral  particle  density, 
and  electron  temperature.  These  quantities  are  used  to  compute  the  electron  collision  frequency.  The 
literature  has  many  conflicting  definitions  of  the  electron  collision  frequency.  These  differences  in 
definition  need  to  be  reconciled  and  validated  with  propagation  measurements.  Second,  improved 
ionospheric  models  need  to  be  added.  This  includes  the  automated  assimilation  of  ionosonde  data.  There 
also  needs  to  include  a  method  to  fit  the  ionosonde  electron  density  profile  to  an  analytical  model  that  can 
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be  effectively  used  by  the  ODE  solver  to  find  the  ray  paths.  In  addition  to  ionosonde  data,  ways  to  include 
traveling  ionospheric  disturbances  and  scintillation  need  to  be  determined.  Finally,  with  the  inclusion  of 
the  types  of  non- stratified  ionospheric  phenomena,  the  ray-homing  algorithm  needs  to  be  improved.  It  is 
necessary  to  generalize  this  algorithm  so  that  accurate  ray  paths  between  fixed  transmitter  and  receiver 
locations  can  be  determined  in  a  fast  and  robust  manner. 
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